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REPORT  SUMMARY 

During  the  period  1  October,  1971  -  30  September,  1972,  the  ongoing 
program  concentrated  on  a  number  of  research  areas  that  were  felt  to  be 
of  significant  and  basic  importance  to  the  successful  development  of 
chemical  and  molecular  laser  technology.  The  principal  areas  of  effort 
were:  measurement  of  V  +  V  and  V  R,T  rates  among  hydrogen  halide 
molecules  and  the  temperature  dependence  of  the  energy  transfer  rates 
previously  determined  for  the  HF,  DF,  HF/CO^,  DF/CC^  systems;  studies  in 
rapid  subsonic  mixing  with  emphasis  on  the  use  of  turbulence  to  enhance 
rapid  mixing;  vibrational  relaxation  studies  of  CC^COO0!)  in  the  pressure 
range  1-100  atm.;  TEA  HBr  chemical  laser  studies;  EPR  study  of  reaction 
rates  in  the  MO  +  F^  dissociative  reaction;  chemical  laser  studies  of 
the  reaction  of  atomic  oxygen  and  acetylene  including  the  determination 
of  V-V  energy  transfer  probabilities  between  CO  molecules;  studies  of 
the  partition  of  exothermicity  in  the  chemical  laser  system; 

experimental  study  of  the  relaxation  of  CO^iOM))  and  HF(v  -  1)  by 
collisions  with  H  and  F  atoms. 

In  all  cases  the  primary  method  of  study  was  the  execution  of 
laboratory  experiments  designed  to  provide  quantitative  results  followed 
by  appropriate  integration  with  theory.  A  detailed  statement  of  results 
and  state  of  progress  will  be  found  in  Sections  1  through  5  of  this  report. 

This  research  is  motivated  by  the  mutual  interest  of  the  sponsor  and 
contractor  in  chemical  lasers  as  efficient  sources  of  coherent  optical 
emission.  Emphasis  is  placed  upon  the  development  of  new  and  efficient 
lasers  operating  at  wavelengths  not  previously  available,  and  upon  the 
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basic  processes  that  underlie  all  molecular  lasers .  The  processes 
quantitatively  studied  included  vibrational  energy  exchange  and  relaxa¬ 
tion  processes,  rates  of  chemical  reactions,  vibrational  distributions 
resulting  from  important  chemical  laser  reactions  and  rapid  gas  mixing 


processes. 


I.  "Vibrational  Energy  Transfer  and  De-excitation  in  the  HF-HC1,  HF-HBr, 
HF-HI  and  HF-DF  Systems"  (T.  A.  Cool) 


Several  reaction  schemes  employed  in  present  cw  HF  chemical  lasers  lead 
to  the  formation  of  vibrationally  excited  HF  in  the  presence  of  HC1,  HBr,  HI 
and  DF.  We  have  measured  several  key  rates  for  energy  transfer  and  deactiva¬ 
tion  for  such  chemical  laser  systems  at  temperatures  of  300  and  350°K.  The 
large  rates  for  these  processes  can  be  attributed  to  energy  transfer  to  rota¬ 
tion  under  the  influence  of  strong  attractive  (hydrogen  bonding)  intermolecular 
forces.  The  energy  transfer  processes  studied  during  the  past  year  are  given 

below  and  the  rate  constants  are  summarized  in  Table  I. 
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It  is  believed  that  all  of  the  above  processes  involve  large  changes  in 

rotational  quantum  number.  This  is  suggested  by  the  large  values  of  the  rate 

constants  despite  the  generally  large  vibrational  energy  defects  for  the  above 

processes.  Moreover,  it  is  apparent  from  the  data  that  the  k'  values  are  too 

e 

large  to  be  explained  unless  energy  transfer  to  rotation  is  important  in  mini¬ 
mizing  the  magnitude  of  the  overall  internal  energy  defect.  Thus,  the  above 
rates  are  actually  an  aggregation  of  vibration  rotation  transitions,  e.g., 

HX(v  =  1,  +  HY (v  =  0,  J2)  ->■  HX(v  =  0,  +  HY(v  =  0,  J2)  (9) 

»  > 

where,  in  general,  Jx  ?  ±  1  and  J2  f  Jj  +  1. 

A  paper  concerning  the  above  work  will  be  presented  at  a  forthcoming  ACS 
meeting  in  January;  additional  work  completed  during  the  current  contract  year 
has  been  published  [J.  Chem.  Phys.,  56,  5863-5878  (1972)]  and  is  included  in 
the  Appendix  to  this  report. 
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TABI.E  I:  MEASURED  RATES 


System 


HF-HC1 


Temperature 

(°K) 

300°K 


Rate  Constants 
(10^  sec*^  torr  H 
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e  lz  - 


ke  +  kn  =  1.55  +  0.1 
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The  laser  excited  fluorescence  method  has  been  employed  to  determine  the  key  rate  constants  for  energy 
transfer  and  deactivation  processes  in  the  HF,  DF,  HF-COi,  and  DF-CO,  chemical  laser  systems  at  a 
temperature  of  350*K.  The  self-deactivation  rates  for  HF(»«1)  and  DF(v-l)  molecules  by  ground 
state  molecules  were  found  to  be  *Hr-nr-5  25±0.30X10'  see*'-torr1  and  kDP-Dr-2.0±0.2X10*  sec*1- 
tore-1,  respectively.  The  measured  rates  of  V—V  transfer  from  HF(s- 1)  and  DF(»-1)  to  the  COj(OOM) 
state  were  4iif-coj-*3.7±0.3X10*  sec'1- tore'1  and  kDr-coi”lI-5=fc2.5X10*  sec'1- torr-1.  The  respective 
deactivation  rates  of  COi(00*l)  by  ground  state  HF  and  DF  were  determined  to  be  kcot-ar~3.6ckO.3X 
10*  sec'*  •  tor  r'*  and  kcot-Dr  “  1 .9^0.4X10*  sec*1-  torr-1.  The  large  rates  for  these  processes  can  be  attributed 
to  energy  transfer  to  rotation  under  the  influence  of  a  sizable  attractive  (hydrogen  bonded)  intermolecular 
potential  well  and  enhanced  repulsion  at  close  range. 


I.  INTRODUCTION 

Much  progress  has  been  achieved  recently  in  the 
development  of  continuous-wave  (cw)  chemical  lasers 
despite  a  serious  lack  of  understanding  of  the  mechan¬ 
isms  of  those  vibrational  energy  transfer  processes 
which  play  major  roles  in  chemical  laser  operation.1-1 
Present  HF  and  DF  cw  chemical  lasers  have  operated 
under  conditions  of  partial  inversion  because  the  rate  of 
redistribution  of  vibrational  energy  exceeds  the  com¬ 
bined  rate  at  which  gas  mixing  and  chemical  reaction 
proceed.  Under  such  conditions  the  antecedent  excita¬ 
tion  process  to  stimulated  emisv'  i  fr  n  a  given  mole¬ 
cule  is  more  likely  to  be  vibrational  energy  transfer 
than  direct  specific  excitation  by  chemical  reaction. 
In  the  DF-CO*  and  HF-COj  “transfer”'  chemical 
lasers,  vibrational  energy  transfer  provides  a  means  for 
the  nonthermal  energy  release  of  chemical  reaction  to  be 
effectively  utilized  as  laser  output  from  totally  in¬ 
verted  COj.7 

Prior  to  the  completion  of  the  present  work,  reliable 
data  were  available  for  the  self-deactivation  of  HF  as 
reported  by  Airey  and  Fried,'  and  more  recently  by 
Solomon  el  al*  and  Bott  and  Cohen.10  No  data  were 
available,  however,  concerning  the  rate  of  vibrational 
deactivation  of  DF  or  the  rates  for  the  important 
processes  by  which  vibrational  energy  is  transferred 
from  HF  and  DF  to  the  COi(00°l)  upper  laser  level. 
Moreover,  the  rates  for  the  vibrational  deactivation  of 
COj(00°l)  molecules  by  collisions  with  HF  and  DF 
molecules  were  unknown  although  the  results  of  con¬ 
current  rate  measurements  are  now  available.11 

We  have  used  the  laser  excited  fluorescence  method" 
to  determine  the  rate  constants  for  these  key  processes 
at  350°K,  a  typical  temperature  for  chemical  laser 
operation.  The  experimental  results  for  HF  and  DF 
vibrational  deactivation  exhibit  good  qualitative  agree¬ 
ment  with  the  vibration-rotation  theory  of  Shin,  when 
the  effects  of  strong  attractive  interactions  are  con- 
sidered.w*-,,b 


The  rates  of  vibrational  energy  transfer  from  HF 
and  r  to  COj  are  too  large  to  be  accounted  for  with 
presen;  approximate  theories.  Interpretation  of  the 
results  presented  here  suggests  that  the  magnitude  of 
these  rates  may  be  explained  in  terms  of  a  hydrogen 
bonded  attractive  intermolecular  potential  well  for 
HF-CQi  and  DF-COj  collision  pairs  which  leads  to 
enhanced  repulsion  at  close  ranges.  A  successful  theory 
must  necessarily  allow  a  major  portion  of  the  difference 
in  vibrational  energy  between  initial  and  final  states  to 
be  taken  up  as  rotational  motion  of  the  collision 
partners. 

The  previously  measured  rate  of  vibration  to  vibra¬ 
tion  (F-»K)  transfer  between  D,(r=  1)  and  CO,(00°1) 
states"  is  shown  to  be  satisfactorily  described  when  the 
vibrational  energy  resonance  defect  is  converted  into 
rotational  motion  of  Dt  under  the  interaction  of  the 
quadrupole  moment  of  deuterium  and  the  asymmetric- 
Stietching  transition  dipole  moment  of  CO*. 

Deactivation  rates  for  the  COj(QO°l)  vibrational 
state  by  HF  and  DF  are  fast  and  appear  to  exhibit  the 
effects  of  a  strong  a!*ractive  interaction  and  vibra¬ 
tion-rotation  energy  tran  sfer. 

II.  EXPERIMENTAL  APPARATUS  AND 
TECHNIQUES 

A.  Laser  Source  and  Absorption  Cell 

The  small  scale  transverse  flow  cw  chemical  laser  and 
the  heated  absorption  cell  used  for  the  present  measure¬ 
ments  have  been  described  elsewhere.’*  Laser  operation 
was  based  upon  the  continuous  mixing  of  F  atoms  and 
Fj  with  Hj  and  D*.u  The  chopper  modulated  output  of 
this  device  consisted  of  triangular  pulses  of  4  psec 
FWHM  with  a  repetition  frequency  that  was  variable 
between  400  and  2500  Hz.  The  total  pulse  energy  was 
typically  4  jrJ;  50%  of  the  pulse  energy  was  present  in 
the  1— >0  band  transitions.  The  standard  deviation  in 
pulse  height  was  less  than  5%;  the  cw  laser  output  was 
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stable  to  1  part  in  20  in  the  2  min  intervals  during 
which  fluorescence  data  were  recorded. 

A  fixed  cell  temperature  of  350±2°K  was  chosen  for 
all  of  the  work  reported  here.  At  the  partial  pressures 
of  IIF  and  DF  used  (0.05  to  1.0  torr)  the  dimer  con¬ 
centrations  in  HF  and  DF  arc  less  than  5  ppm  in 
equilibrium  at  this  temperature.1*  The  cell  temperature 
was  monitored  with  a  stainless  steel  jacketed  thermo¬ 
couple  placed  in  direct  contact  with  the  gas  within  the 
absorption  cell.  The  thermocouple  junction  was  located 
near  the  center  of  the  detector  field  of  view  about  2.5 
mm  outside  of  the  path  of  the  laser  beam  and  about 
1  mm  away  from  the  nearest  wall  surface. 

Cell  pressures  were  measured  with  a  0-10  torr  MKS 
capacitance  manometer  which  was  accurate  to  within 
0.1%.  The  calibration  of  the  manometer  electronics 
package  supplied  by  the  manufacturer  was  confirmed 
with  a  digital  voltmeter  capable  of  an  accuracy  of  1 
part  in  10s.  The  total  cell  pressure  was  always  greater 
than  10  torr  to  prevent  excessive  diffusion  of  vibra- 
tionally  excited  species  to  the  cell  walls. 

The  gas  sample  was  mixed  and  preheated  before 
entering  the  heated  cell.  This  was  accomplished  by 
injecting  room  temperature  HF'  (DF)  and  CO»  through 
needle  valves  into  a  heated  argon  flow.  The  combined 
flow  passed  along  aim  length  of  heated  stainless  steel 
tube  of  4  mm  bore  with  a  mean  residence  time  exceeding 
1.5  sec  to  provide  adequate  mixing  and  equilibration 
before  entering  the  cell.  The  heated  cell  was  designed 
to  provide  a  laminar  flow  without  recirculation  or 
separation  in  the  cell  region  within  the  detector  field 
of  view.  Typical  gas  flow  rates  through  the  cell  were 
0.020  mM/scc.  A  mean  residence  time  within  the  cell 
was  0.2  sec.  Mean  cell  flow  velocities  were  20  cm/sec. 


These  flow  rates  were  sufficient  to  minimize  the  effects 
of  surface  reactions  in  causing  adsorption  or  contamin¬ 
ation  of  the  HF  or  DF  gas  samples  at  the  cell  pressures 
and  temperatures  employed  in  these  studies.  The  flow 
velocities  were  low  enough,  however,  to  avoid  con¬ 
vective  effects  which  would  lead  to  an  apparent  shorten¬ 
ing  of  the  fluorescence  lifetime  because  of  the  sweeping 
of  excited  molecules  outside  of  the  detector  field  of  view. 

The  field  of  view  of  the  detector  included  a  sample 
region  extending  approximately  2  cm  along  the  laser 
optical  axis  from  the  beam  entrance  window  (most  of 
the  signal  originated  from  the  region  nearest  the 
entrance  window).  Narrow  band  interference  filters 
were  placed  between  the  detector  and  cell  to  eliminate 
the  scattered  P-branch  HF  or  DF'  laser  radiation. 
F'ilters  for  DF  with  a  peak  transmission  of  60%  for  a 
pass  band  located  between  the  1%  transmission 
limits  of  3.228  and  3.423  microns  were  used;  for  HF, 
filters  with  a  peak  transmission  of  75%  at  2.47  y  and 
1%  transmission  limits  at  2.27  and  2.64  microns  were 
used.  These  filters  passed  most  of  the  fluorescence  on 
the  R  branches  for  the  (t>=  1— *0)  HF'  and  DF  bands. 
Fluorescence  from  the  CO2(n0nl)  level  was  sampled 
through  filters  with  75%  peak  transmission  at  4.3  y 
and  a  F'VVHM  of  0.4  y. 

The  absorbed  laser  pulse  energies  were  about  three 
orders  of  magnitude  below  the  values  reported  by  other 
investigators  in  single  pulse  experiments.8-17  The  frac¬ 
tion  of  HF  (DF')  molecules  excited  did  not  exceed 
0.01  and  thus  no  special  precautions  were  necessary  to 
prevent  excessive  sample  heating  or  excitation  of  the 
second  vibrational  level  of  HF'  (DF')  by  the  laser 
pulses. 

B.  Gas  Handling  and  Impurities 

The  present  gas  handling  system  shown  schematically 
in  Fig.  1  was  similar  to  that  employed  by  Airev  and 
Fried.8  Those  portions  of  the  system  exposed  to  HF 
were  of  stainless  steel  (except  the  aluminum  absorption 
cell)  and  were  carefully  passivated  to  insure  repro¬ 
ducible  results.  The  absorption  cell  and  HF  handling 
system  were  exposed  to  100  torr  of  HF  for  24  h.  The 
distillation  cylinders  (Ci  and  C»)  were  additionally 
exposed  to  900  torr  of  HF  for  one  week.  The  entire 
system  (including  argon  and  COj  handling  portions) 
was  pumped  down  to  well  below  10~8  torr  and  baked 
out  at  90°C  for  several  clays.  After  this  the  leak  rate  at 
room  temperature  for  the  entire  system  of  1.5  liter 
volume  was  not  detectable  (less  than  10"*  torr/h) . 
A  leak  rate  of  5X10"*  torr/h  was  measured  for  the 
0.1  liter  volume  absorption  cell  and  connecting  lines. 

The  absorption  cell  and  connecting  lines  were  then 
brought  up  to  the  operating  temperature  of  350°K  and 
a  flow  of  10  torr  argon  and  1  torr  of  HF  was  maintained 
for  several  hours.  This  resulted  in  an  increased  out- 
gassing  rate  which  stabilized  at  a  value  of  about  0.01 
torr/min.  No  further  significant  change  from  this  out- 
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gassing  rate  was  observed  in  the  entire  series  of  experi- 
ments  on  HF.  When  all  1IF  data  were  recorded,  the 
system  was  prepared  again  according  to  the  above 
procedures  before  any  DF  data  were  taken;  similar 
leak/outgassing  rates  were  measured  for  DF. 

The  major  impurities  (in  parts  per  million)  of  the 
gases  used  in  this  study  were  as  follows:  argon:  Mathe- 
son  UHP  (02<2,  Nj<4,  C02<1,  HjO<3);  C02: 
Air  Products  Research  Grade  (O2+N2=10,  r0<10, 
CH|=1A  HjO<1.5)  ;  HF:  Air  Products  CP  Grade, 
typical  analysis  (H;SiFe=30,  SO2=40,  H2=300); 
DF:  Matheson  CP  grade,  typical  analysis  (98%  min¬ 
imum  purity,  HF<2%).  The  HF  and  DF  as  actually 
supplied  contained  about  20%  H2  and  D2,  respectively. 

Additional  purification  of  gases  was  accomplished  as 
follows.  A  dry-ice  cold  trap  (see  Fig.  1)  was  employed 
in  the  argon  flow  to  minimize  the  water  vapor  impurity. 
The  HF  and  DF  were  purified  to  remove  H2  (D2) 
and  H:0  (D^)  according  to  the  method  of  Airey  and 
Fried.8  With  this  procedure  a  5  liter  sample  of  HF(DF) 
gas  was  admitted  to  the  distillation  cylinder  (Ci  of 
Fig.  1)  and  was  solidified  at  77°K.  Gaseous  H2  (D2) 
was  then  removed  by  pumping  the  distillation  cylinder 
to  below  0.010  torr.  The  cylinder  was  then  allowed  to 
warm  to  195°K  and  a  trap  to  trap  distillation  between 
cylinder  C'i  and  a  second  cylinder  C2  cooled  to  77°K 
was  performed  at  least  twice  before  data  were  taken. 
In  each  stage  of  the  distillation  the  last  5%  of  the 
original  sample  was  discarded;  the  empty  cylinder 
was  baked  out  under  vacuum  at  90°C  for  24  h  before 
the  next  distillation.  No  detectable  difference  was  ob¬ 
served  between  fluorescence  trices  taken  after  six 
dist  'lations  and  those  taken  after  two  distillations. 

T  ie  procedure  used  for  establishing  a  given  absorp¬ 
tion  cell  gas  composition  was  as  follows.  First  the  argon 
flow  rate  was  set  to  give  the  desired  argon  cell,  pressure 
(usually  10  torr)  as  indicated  by  the  capacitance 
manometer  M.  Under  these  conditions  choked  flow 
existed  through  both  valves  V2  and  V3;  the  pressure 
ratio  across  these  valves  always  exceeded  10: 1  and 
15:1,  respectively.  The  manometer  was  then  balanced 
by  opening  valve  V)  to  admit  a  reference  air  flow  and 
by  adjustment  of  Vi  until  the  null  position  for  the 
manometer  diaphragm  was  established.  This  balance 
was  stable  to  within  ±10/i  during  the  3-4  min  period  of 
a  tvpical  run.  Absolute  argon  pressures  were  obtained  to 
within  1%  with  pressure  gauge  G  (Wallace  and 
Tiernan,  0-  20,  or  Q-100  torr)  after  corrections  for  any 
capacitance  manometer  imbalance  were  made. 

A  given  partial  pressure  of  C02  was  then  established 
within  the  cell  when  V«  was  opened  and  adjusted  to 
give  the  desired  pressure  increase  as  indicated  by 
the  manometer  M.  In  a  similar  fashion  a  desired 
pressure  of  HF  (DF)  was  finally  admitted  by  means 
of  valve  V$.  The  flows  through  valves  Vs  and  V4  were 
always  choked.  The  partial  pressure  of  C02  could  be 
maintained  to  within  ±10/i  at  the  lower  range  of 
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Fjc.  2.  Vibrational  energy  level  diagram  which  illustrates  the 
laser  excited  fluorescenre  transitions  in  HF,  DE,  and  COj. 


pressures  employed  or  ±2%  at  the  higher  pressures 
used.  For  HF(DF)  the  corresponding  accuracies  of 
partial  pressure  settings  were  ±20m  or  ±5%. 

After  the  run,  valves  V4  and  Vs  were  turned  off 
separately  and  the  pressure  changes  indicated  by  the 
manometer  were  recorded.  These  values  were  in  agree¬ 
ment  with  the  initially  recorded  pressures  to  the  ac¬ 
curacies  stated  above;  an  average  of  the  initial  and  final 
values  was  usually  used  in  data  reduction.  The  sum  of 
HF(DF)  and  C02  partial  pressures  was  monitored  with 
the  capacitance  manometer  throughout  the  run. 
Occasionally  a  drift  in  pressure  would  occur  during  a 
run  which  was  clearly  outside  the  above  stated  error 
limits;  such  runs  were  terminated  and  the  fluorescence 
data  discarded. 

C.  Instrumentation  for  Fluorescence  Detection 

The  electronics  instrumentation  employed  for  this 
work  has  already  been  described.14  The  ouput  of  the 
indium  antimonide  PV  detector  was  amplified  and  used 
to  drive  either  a  P.A.R.  Model  CW-1  boxcar  integrator, 
or  Model  TDH-9  waveform  eductor.  The  time  bases 
for  both  instruments  were  determined  to  within  0.5?,' 
wilh  a  100  kHz  crystal  calibrator.  The  scanning  circuilb 
of  the  boxcar  integrator/waveform  eductor  were 
triggered  by  a  phototransistor  synchronized  to  the 
laser  pulse.  The  time  origin  for  the  fluorescence  traces 
coincided  with  tiie  peak  of  the  baser  pulse  and  could  be 
determined  to  within  ±0.5  /isec. 

Each  fluorescence  trace  represented  the  average  of  as 
many  as  300  000  separate  event?.  Before  signal  averag¬ 
ing,  the  signal  to  noise  (S/N)  ratios  ranged  from  a  high 
of  3/1  for  C02  fluorescence  from  the  DF-C02  system 
to  about  1/30  for  DF  fluorescence  from  the  DF  and 
DF-C02  systems.  All  of  the  HF  data  and  most  of  the 
DF  data  were  taken  with  the  boxcar  integrator,  which 
provided  an  improvement  in  the  S/N  ratio  by  factors 
ranging  from  30  to  80.  The  weaker  DF  fluorescence 
from  the  DF  and  DF-C02  sys  .-nis  (owing  to  the  rela¬ 
tively  long  spontaneous  emission  lifetime  of  DF  com¬ 
pared  with  both  HF  and  C02)18  usually  required  use  of 
the  waveform  eductor  which  gave  an  improvement  in 
S/N  ratio  by  a  factor  exceeding  350.  After  signal 
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averaging,  the  S/N  ratio  of  these  data  generally  ranged 
r.  om  10/1  to  100/1. 

During  the  2  min  time  period  of  a  typical  run,  the 
base  line  of  the  fluorescence  trace  was  usually  constant 
to  within  3%;  all  data  which  exhibited  base  line  shifts 
of  more  than  5%  were  discarded.  The  observation  that 
data  taken  with  the  waveform  eductor  were  in  excel¬ 
lent  agreement  with  those  data  taken  with  the  boxcar 
integrator  confirmed  that  the  stability  of  the  laser 
output  during  a  given  run  was  adequate  for  the  present 
measurements. 

III.  EXPERIMENTAL  RESULTS 

A.  Vibrational  Energy  Transfer  Processes 

The  processes  of  importance  in  describing  the  present 
measurements  of  vibrational  relaxation  in  DF  and 
mixtures  of  DF  and  C02  are  the  following  (refer  to  the 
energy  level  diagram  of  Fig.  2) : 

vibration  to  rotation  (V— >R)  energy  transfer  in  DF: 

DF(B»i)y1')+DF(#=olyI')-^DF(t-o,yI) 

+DF(t'=0, 72)+A/*!=  2907  cm-1  (1) 

(where  A E  is  the  vibrational  energy  difference  between 
initial  and  final  states) ; 

vibration  to  vibration  (V— *V)  transfer  between 
DF  and  COj: 

k. 

DF(t=  1,-70  +  C02(00°0)  *=>  DF(»=  0,  J)  +COJ(00°1) 

k,' 

+  A£=558  cm-1;  (2) 

de-excitation  of  the  upper  COj  laser  level  by  DF 
through  intramolecular  vibrational  energy  transfer 
in  COj  and  energy  transfer  to  rotation  in  DF: 

7'io 

COj(00Dl)-fDF(ti=0,  J')—*  CO2(wm'0) 

-f-DF(»=0,  J);  (3) 
intramolecular  energy  transfer  in  C02: 


and  fluorescence  from  the  C02[(0f)°l)— *(00<'0)]4.3/i 
band 


Trj 

CO2(00"l)  —  CO2(00('0) +/;■/,,  A/=±l.  (7) 

In  Processes  (l)-(3)  only  certain  assignments  of  the 
rotational  quantum  numbers  will  not  require  a  large 
energy  transfer  to  translational  motion.  For  present 
purposes  we  will  consider  the  above  rate  constants  as 
total  rates  summed  over  all  rotational  quantum 
number  combinations  which  provide  an  appreciable 
contiibution.  It  is  assumed  that  equilibrium  exists  for 
rotational  states  within  a  given  vibrational  level. 

It  has  been  proposed7  that  the  coupling  of  vibrational 
energy  between  HF(y=l)  and  CO2(00°l)  levels  may 
involve  the  participation  of  an  intermediate  COj 
state  in  contrast  to  the  direct  coupling  of  Process  (2) 
indicated  for  the  DF-C02  system  by  the  present  ex¬ 
perimental  results.19  Therefore,  in  the  HF-C02  case 
one  must  consider  both  the  direct  coupling  processes 


k. 

HF(r=  1,  y')  +  CO,((XH))  — ♦  HF(r=0,y) 

+  CO2(00°l)-f  AE=  1609  cm-1  (2) 
and  the  alternate  processes 


HF(t‘=  1)  +  C02(00’0)  HF(i)=0)-+-COj(10°l) 

ky, 

+AE  =  242  cm-1  (2') 
or 

kbt 

HF(r=  1)  +  C0j(00°0)  <=>  HF(t,=0)-fCOj(02°l) 

kfb 

+AE=349  cm-1  (2') 
and 

C02  ( 10° . )  +  C02  ( OCW)  CO,  ( 10°0)  -f  COj  (00»1 ) 


+AE=-21cm-1  (2") 
or 

ktl 

C0j(02nl)  +  C0j(00n0)  -+ C0j(02°0)-|-C0j(00nl) 


COj(00°l)-)-M  -4  CO2(»tw'0),+M  (4) 

(where  the  species  M  is  either  C02  or  argon,  with  the 
respective  rate  constants  kloc  and  ku>a) ; 

vibration  to  translation  ( V—*T)  deactivation  of  DF 
by  argon  ( V-*T  deactivation  of  DF  and  HF  by  C02 
is  assumed  negligible  under  the  conditions  of  the 
present  experiments.) : 

kbaa 

DF(»=  1)+Ar  -»  DF(t>=0)+Ar;  (5) 

fluorescence  from  the  DF(n=l— >0)  fundamental 
band 

Tn 

DF(t>=  1)  —>  DF(r=0)-f-Ai/1(  AJ  —  ±1,  (6) 


+A£—  —26  cm-1.  (2") 

Processes  analogous  to  those  of  Eqs.  (1)  and  (3)— (7) 
also  apply  to  the  HF  and  HF-CO?  systems.  For  the 
HF-COj  system  we  must  consider  the  additional 
fluorescences  of  the  “hot”  bands” 

CO2[(10"l),  (02°1)-4(10°0),  (02"0) 

A/=±l.  (8) 

Processes  (8)  present  a  major  experimental  difficulty 
in  the  determination  of  which  of  the  coupling  mechan¬ 
isms  (2)  or  (2')  and  (2")  gives  rise  to  the  observed 
f  orescence  in  the  4.3m  region.  Recent  high  resolution 
absorption  measurements”  indicate  that  the  matrix 
elements  for  these  bands  ate  at  least  as  large  as  that 
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for  the  C03[(00°I)->(00°0)]  band  of  Process  (7). 
Thus,  since  fluorescence  from  the  bands  associated 
with  both  Processes  (7)  and  (8)  would  lie  within  the 
pass  band  of  the  4.3p  interference  filter  used  here,  we 
have  no  means  in  the  present  experiments  to  distinguish 
between  the  alternate  coupling  mechanisms. 

B.  Self -Deactivation  of  HF  and  DF 

1.  Vibrational  Deactivation  of  DF(v=l)  by  DF( v=0) 

Fluorescence  from  the  fundamental  v=  1— >0  band  of 
DF  was  observed  from  mixtures  of  DF  with  argon.  The 
intensity  of  this  fluorescence  was  proportional  to  the 


T ime,  (,u  sec) 


Time,  (/»iec) 

Fro.  3.  A:  Decay  of  DF(r=  I)  — *DF (v=0)  fluorescence  from 
the  mixture:  DF- 0.305  torr,  argon  =  9,97  torr,  at  349°K.  B: 
Semilogarithmic  graph  of  a  smooth  curve  fitted  to  the  data  of 
Fig.  3A. 

DF(v=  1)  concentration  given  by  the  expression 

[DF(#=  l)]/CDF]o=exp(  — f/r),  (9) 

where  [DF]0  is  the  initial  DF(v=*  1)  population  pro¬ 
duced  by  absorption  from  the  DF  laser  pulse,  and 

{pr)~'=  k  DF-nfX  dH-^df-a  *X Ar+  ( prn)~l .  (10) 

Here  p  is  the  total  pressure  of  the  gas  mixture;  r  is  the 
relaxation  time,  r„  is  the  total  radiative  lifetime  cf  the 
v=  1  state  of  DF1S  and  f'fnr  and  X&,  are  the  respective 
mole  fractions  of  DF  and  argon. 

Figure  3A  shows  a  typical  DF  fluorescence  trace; 
such  DF  relaxation  data  always  exhibited  a  single 
exponential  time  constant  as  is  illustrated  for  the  data 
of  Fig.  3A  Iv  the  semilogarithmic  plot  of  Fig.  3B. 
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Fig.  4.  Observed  DF  fluorescence  decay  rates  for  Dl'-argon 
mixtures  at  350°K  as  a  function  of  DF  mole  fraction. 

Figure  4  contains  all  of  the  fluorescence  decay  data 
taken  for  mixtures  of  DF  and  argon.  Here  values  of 
obtained  from  the  slopes  of  semilogarithmic 
plots  of  the  data  for  each  run  are  plotted  against  the 
mole  fraction  of  DF  in  the  conventional  manner.  Most 
of  the  data  were  taken  for  an  argon  partial  pressure  of 
10  torr;  however,  a  few  points  were  taken  at  an  argon 
pressure  of  100  torr  to  establish  the  location  of  the 
intercept  of  the  data  with  the  ordinate.  Since  the 
intercept  of  Fig.  4  is  at  the  origin  within  experimental 
uncertainty,21  the  neglect  of  the  last  two  terms  of  Flq. 
(10)  is  justifiable  in  agreement  with  analogous  findings 
from  other  studies  of  hydrogen  halides.8  "  21  The  DF 
partial  pressures  of  Fig.  4  range  from  0.05  torr  to  1.05 
torr.  This  range  was  selected  because  fluorescence 
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Fic.  5.  A:  Decay  of  HF(v—  1)— »HF(r  =  0)  fluorescence  from 
the  mixture:  HF- 0.465  torr,  argon-9.75  torr,  at  351°K.  It: 
Semilogarithmic  graph  of  a  smooth  curve  fitted  to  the  data  of 
Fig.  5.\. 
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ElC.  6.  Observed  HE  fluorescence  decay  rates  for  H E-argon 

mixtures  at  350°K  as  a  function  of  ill'  mole  fraction. 

signals  were  significantly  weaker  at  higher  and  lower 
pressures.  In  addition,  at  somewhat  higher  pressures  the 
accuracy  of  measurement  would  be  limited  by  the  time 
response  of  our  apparatus;  at  lower  pressures  the 
partial  pressure  determinations  would  be  less  certain. 
The  data  of  Fig.  4  give  ltba^kn f  df-  2.O±O.2X104 
sec"  ‘•tori--1  at  550° K. 

2.  Vibrational  Deactivation  of  11  F(v-\)  by  IlF{v-0) 

Fluorescence  from  the  2.7/x  v=  1— *0  band  of  HF  was 
observed  and  analyzed  for  mixtures  of  HF'  with  argon 
in  a  manner  similar  to  that  discussed  for  DF.  Single 
exponenlial  decays  such  as  shown  in  Figs.  5A  and  5B 
were  observed  for  all  data.  Figure  6  shows  the  relaxation 
times  for  all  of  the  HF"  fluorescence  data  as  a  function 
of  the  mole  fraction  of  HF.  The  smaller  scatter  of  these 
data  compared  to  those  of  Fig.  4  probably  reflects  the 
shorter  radiative  lifetime  of  IIF18  which  affords  a  better 
fluorescence  S/N  ratio  than  was  achieved  for  DF 
(compare  Figs.  3A  and  5A).  Most  of  the  data  were 
taken  at  a  10  torr  argon  pressure  with  HF  partial 
pressures  ranging  from  0.05  to  1.06  torr;  a  few  ad¬ 
ditional  traces  were  recorded  for  100  torr  of  argon  and 
HF  pressures  from  0.10  to  0.35  torr  to  confirm  the 
negligible  effect  of  Process  (5)  compared  with  Process 
(1)  as  was  demonstrated  for  DF  with  Fig.  4.21  The 
slope  of  the  linear  fit  to  the  data  of  Fig.  6  gives  the  self- 
deactivation  rate  of  HF  as  ^im=^hf-hf=5.25± 
.30X10*  sec_l-torr-1  at  350°K.  Our  rate  constant  is 
about  44%  below  the  deactivation  rate  for  HF’  previ¬ 
ously  measured  by  Airey  and  Fried.8  Since  the  gas 
handling  and  purification  procedures  followed  by  the 
two  groups  were  virtually  identical,  it  is  difficult  to 
account  for  this  discrepancy  which  is  somewhat  outside 
our  assigned  error  limits. 

C.  Measurements  of  Vibrational  Energy  T.ansfer  in 
the  HF-COj  and  DF-COj  Systems 

/.  Energy  Transfer  in  DF-COi  Mixtures 


from  the  4.3/t  (00°1— >00°0)  band  of  C02.  It  was  thus 
possible  to  monitor  the  concentrations  of  both  the 
DF(b=1)  and  CO2(00°l)  molecules  following  the  DF' 
laser  pulse.  The  time  dependence  of  these  concentra¬ 
tions  is  given  by  the  solutions  of  the  rate  equations 
corresponding  to  the  processes  of  Eqs.  (l)-(7).  The 
solutions  are  given  by  tile  sum  of  two  exponential 
terms12: 

[DF(„=l)]/CDF]«=C(A1-fll)/(X.-Xl)] 

X  exp ( - X2f)  +  [ (fl! - X2)  / (X, -  X2)  ]  exp ( -X,0 ,  (11) 

[C02 (00°1 ) ]/(]D F]0  =  [ (Xi — a, )  ( a, — X2)  'A,(X,-X2)] 

X[exp(— XjO  —  e.xp(— Xff)],  (12) 

where  ai/p~  kt,aXm:-{-krXcoi\-kuaX/u-\-(pT,l)~>t  and 

bi/p  —  ke'Xar, 

and 


(hx)/P~  J {<irl  ^  Ft(<ri— th)*  f-4«26IJ1'2) ,  (13) 


'hi  p  —  k,Xcnt 


bt/ p  —  (£,'  +  £lo)  A,DF+£lOC-XrcOi+Fio°ATAr+  (PT't)~'- 


In  Eq.  (13)  the  square  root  term  is  always  positive  and 
X[  corresponds  to  the  +  sign.  It  can  readily  be  seen 
that  Xi  and  X2  are  positive;  Xi>0]  and  X2<a,.12 

The  rate  constant  k,  for  the  transfer  of  vibrational 
energy  from  DF  to  C02  by  Process  (2)  can  be  de¬ 
termined  most  easily  from  observations  of  DF(v--l) 
fluorescence  under  conditions  of  small  concentrations  of 
DF.  In  the  limiting  case  when  the  concentration  of  DF' 
is  sufficiently  small,  Eq.  (11)  reduces  to  the  simp! 
form  of  Eq.  (9) ,  where  now  the  relaxation  time  is  given 
by  the  expression 


{pr)~l  —  kbaX-DV+k  cXcOi+kbfX  K,-\-  (pTri)~' .  (  14) 

Only  the  first  two  terms  on  the  rhs  of  (14)  are  needed 
to  interpret  the  data  since,  as  before,  the  last  two  terms 
are  of  negligible  magnitude  (thus,  in  the  discussion  of 
the  DF-COj  and  HF-C02  systems,  the  reference 


Fluorescence  was  obseived  from  mixtures  of  DF  and 
COi  in  argon  on  both  the  (»=1— »0)  band  of  DF  and 


Eig.  7.  Observed  DE  fluorescence  decay  rates  for  DE-CO>- 
argon  mixtures  at  350°K  at  a  function  of  the  reduced  mole  fraction 
of  COj;  .Ycoi,**-X’coi/(-X'Dr+i'fcoi)- 
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pressure  p  is  the  sum  of  only  the  C02  and  DF  (HF) 
partial  pressures).  Figure  7  shows  relaxation  rate  data 
(solid  points)  for  various  values  of  *W—  (Xco,)/ 
(Xdf+Xco,)  ranging  from  0.46  to  0.84  (the  open  data 
points  and  the  dashed  lines  of  Fig.  7  are  discussed  at  the 
end  of  this  section) .  These  data  exhibited  a  well  de¬ 
fined  single  exponential  behavior  and  {pf)  1  values 
were  obtained  from  scmilogarithmic  plots  in  a  manner 
similar  to  that  discussed  in  connection  with  Figs.  3A 
and  3B.  Extrapolation  of  these  data  with  the  solid  line 
of  Fig.  7  to  the  value  Xco/^l  yields  a  value  for  the 
rate  constant  for  V-+V  transfer  from  DF(r=l)  to 
CO,(00°1)  of  k\  =1.75  ±0.25  X105  sec'1 -torr"1.  [The 
solid  line  of  Fig.  7  is  not  simply  a  reasonable  fit  to  the 
data  shown,  but  a’so  has  a  slope  constrained  by  the 
values  of  k,'  and  km  as  will  become  apparent  after  some 
further  discussion.] 

Equation  (12)  describes  the  rise  and  subsequent 
decay  of  CO2(00°l)  concentration.  The  rate  constant 
kt0  for  the  de-excitation  of  CO2((J0°l)  by  DF(»=0) 
was  measured  from  examination  of  CO2(00°l)  fluor¬ 
escence  for  times  sufficiently  long  aft'-,  the  initiating 
laser  pulse  that  the  decay  of  CP^00°1)  population 
could  be  accurately  described  by  the  single  exponential 
term  exp(  — X2/)  appearing  in  Eq.  (12).  That  is,  when 
exp(-X5/)»exp(-Xi/),  then  (£r)“‘  =  X2//>.  When  X2 


Time,  (»see) 


Fig.  8.  A:  Decay  of  CO2[(00°l)-»(W0) ]  fluorescence  during 
ihc  latter  stages  of  relaxation  for  Ihe  mixture:  DI —0.340  torr, 
C  0,-0.720  torr,  argon-9.91  torr,  at  351'K.  B:  Semilogatithmic 
graph  of  a  smooth  curve  fitted  to  the  data  of  lig.  8A. 


1. 1  ,  i .x. — • - J 

Q1  r*  10 

Fic.  9.  Observed  CO,  fluorescence  decay  rates  for  the  latter 
stages  of  relaxation  for  DF-CO«-argon  mixtures  at  350  K  as  a 
function  0?  the  reduced  mole  fraction  of  DF;XDF,“ADF/(3tD»+ 
Xco,). 


was  computed  with  the  Eqs.  (13),  the  terms  and 
(pTr,)~l  were  again  dropped;  the  term  kwm  has  a  small 
value53  of  50  sec-1 -torr-1  and  was  also  neglected;  the 
value  of  r,,  is  long  enough1*  that  the  last  term  appearing 
in  the  expression  for  h  can  usually  be  neglected  with¬ 
out  appreciable  error.54 

A  representative  CO2(00°l)  fluorescence  trace  (Fig. 
8A)  was  used  to  construct  Fig.  8B,  which  illustrates  the 
linear  time  dependence  of  the  logarithm  of  fluorescence 
intensity  observed  for  all  DF-CO,  mixtures  during  the 
latter  stages  of  relaxation.  Figure  9  shows  the  observed 
dependence  of  (pr)“'  values  determined  from  such 
plots  upon  the  ratio  Xm',  where  Xdf'^(Xdf)/ 
(XDf+Xco,),  for  all  of  the  C02  fluorescence  data 
The  cell  compositions  were  systematically  chosen  to 
cover  the  partial  pressure  ranges  for  DF  from  0.05  to 
1.0  torr,  for  CO,  from  0.05  to  10.0  torr,  and  for  argon 
from  10  to  20  torr. 

The  value  for  the  rate  constant  tu  may  be  deduced 
from  the  data  of  Fig.  9.  In  Eqs.  (13)  the  rate  constant 
kioc  is  known  to  be  335  sec-1 -torr-1  at  350  K53;  the 
values  of  and  k,  have  been  given  above.  Dctadcd 
balancing  provides  the  relationship  kc'=k,  exp(  — 55S 
hc/kT) ;  where,  as  has  been  stated,  k,  and  k,'  are  over-all 
rate  constants  (summed  on  J  and  J')  and  rotational 
equilibrium  at  the  temperature  T  is  assumed.  Use  of 
these  rate  constants  and  the  relationship,  ( pr)-l  =  \}/p , 
for  three  different  assumed  values  for  km  leads  to  the 
curves  shown  with  the  data  of  Fig.  9.  The  solid  curve  has 
been  drawn  for  the  value  km—  1.9X  10*  sec  **torr  ‘.The 
upper  and  lower  dashed  curves  are  for  £10 =  2.3X  101 
sec"1 -torr1  and  1.5X10*  sec"1 -torr1,  respectively; 
these  curves  represent  reasonable  upper  and  lower 
bounds  on  km  when  the  data  of  Fig.  9  and  the  error 
limits  on  k„  k,',  k^,  and  kmc  are  considered.  Thus,  we 
find  the  value  for  the  rate  of  de-excitation  of  the  upper 
laser  level  of  C02  by  DF  to  be  ^10=  1.9±0.4X  10* 
sec-1  •  torr-1  at  35G°K. 

The  foregoing  determination  of  the  rate  constants 
kt„,  k„  and  for  the  processes  of  Eqs.  (l)-(3)  has 
been  facilitated  by  the  examination  of  certain  ex- 
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exponential  decay  terms.  The  slow  decay  time  (r  =  108  ^sec)  is 
determined  from  a  straight  line  fit  to  the  data  of  the  latter  stages 
of  relaxation.  The  difference  between  this  line  and  the  DF  fl  lorcs- 
cence  observed  early  in  the  relaxation  yields  a  second  si  aight 
line  with  a  decay  time  r-  28.8  *isec. 

perimental  conditions  for  which  the  fluorescence 
signals  may  be  represented  as  single  exponential  decays. 
It  is  appropriate  to  now  verify  that  the  use  of  the  rate 
constant  values  thus  determined  docs  in  fact  pernit  an 
accurate  description  of  the  observed  fluorescence  traces 
throughout  the  entire  accessible  range  of  experimental 
parameters.  Figure  10  shows  a  semilogarithmic  plot  of 
DF(»=  1)  fluorescence  data  from  a  trace  taken  for  the 
composition:  DF:  0.51  torr,  CO2:  0.11  torr,  and  argon: 
9.72  torr.  Here  a  double  exponential  behavior  is  clearly 
resolvable.  From  Fig.  10  we  measure  directly  the 
values:  Xi/^= 5.6X10*  sec-1  •torr-1  and  h/p-  1.5X 
104  sec-1,  torr-1.  These  values  may  be  compared  with 
calculated  values  based  upon  the  measured  rate 
constants  and  Eq.  (13)  of  Xi/^*=  6.2X10*  sec-1*torr-1 
and  X* /p<=  1.6X10*  sec-1 -torr-1. 

A  further  illustration  of  the  good  agreement  between 
experimental  fluorescence  traces  and  calculated  fluor¬ 
escence  behavior  is  given  in  Figs.  11A,  11B.  and  11C. 
These  traces  were  all  taken  for  essentially  the  same  gas 
composition.  In  order  to  adequately  resolve  both  early 
and  late  time  behavior  of  the  CO*  fluorescence,  Figs. 
11 A  and  1  IB  were  taken  for  identical  conditions  but 
were  recorded  with  differing  time  bases.  Figure  1 1C 
gives  the  corresponding  DF  fluorescence  decay.  The 
curves  shown  with  the  data  were  calculated  from 
Eqs.  ( 1 1)— ( 13)  with  the  rate  constants  quoted  above 
and  were  normalised  to  the  experimental  traces.  Not 
all  of  the  fluorescence  traces  showed  the  excellent 
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agreement  with  calculations  illustrated  in  Figs.  1 1  A— 
1 1C;  however,  in  every  case  the  agreement  was  within 
the  uncertainty  associated  with  the  data  point  scatter 
of  Figs.  4,  7,  and  9. 

In  concluding  this  section,  we  must  confirm  the 
validity  of  the  extrapolation  made  in  Fig.  7  in  the  de¬ 
termination  of  the  rate  constant  k,  from  DF  fluor¬ 
escence  data.  The  extrapolation  was  based  on  the 
fact  that  single  exponential  behavior  was  observed  for 
the  solid  points  of  Fig.  7,  and  upon  the  observation 
that  as  Xcot  tends  toward  unity  (Ydf'— *0)  the 
rhs  of  Eq.  (11)  approaches  the  single  term:  exp(— ajf). 
We  are  now  in  a  position  to  show  that  the  solid  line  of 
Fig.  7  is  correctly  drawn  within  narrow  limits.  Since 
the  solid  data  points  of  Fig.  7  represent  single  ex¬ 
ponential  decays,  these  (pr)~'  values  must  correspond 
to  either  Xi /p  or  \?/p  of  Eqs.  (11)  and  (12).  The 
upper  and  lower  dashed  lines  of  Fig.  7  show  the  re¬ 
spective  variations  with  of  the  values  of  \\/p 

and  X*/ p  calculated  from  our  rate  constant  values  and 
Eq.  (13).  The  right-hand  portion  of  the  upper  dashed 
line  coincides  with  the  solid  curve  drawn  to  determine 
k,.  A  direct  comparison  w’th  experiment  is  provided 


Fig.  11.  Observed  and  calculated  DF  and  CO*  fluorescence 
from  a  DF-COr-*rgon  mixture  at  350*k.  A:  CO*  fluorescence 
on  a  long  time  scale  for  the  mixture:  DF-0.290  torr,  COi- 
0.105  torr,  argon  - 19.7  torr.  B:  COj  fluorescence  on  a  short  time 
scale  for  the  same  conditions  as  in  A).  C:  DF  fluorescence  for 
the  mixture:  DF-U.305  toir,  CO*— 0.105  torr,  argon-20.0 
torr. 
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from  two  DF  fluorescence  traces  which  exhibited 
clear'y  resolvable  double  exponential  behavior.  The 
two  pairs  of  open  data  points  of  Fig.  7  give  the  \Jp 
and  A  i/p  values  obtained  from  these  traces  in  the 
manner  illustrated  in  Fig.  10  (which  was  drawn  from 
one  of  the  two  DF  traces) .  The  upper  two  data  points 
are  in  agreement  with  both  the  calculated  \\/p  curve 
and  the  other  DF  fluorescence  data  (solid  points) ; 
the  lower  two  points  show  good  agreement  with  the 
calculated  A  ,/p  curve  and  hence  also  with  the  CO, 
fluorescence  data  of  Fig.  9  (the  solid  curve  of  Fig.  9 
is  the  A  i/p  curve  of  Fig.  7) .  We  note  further  that  the 
open  data  points  serve  to  confirm  the  Afco,'  =  0  in¬ 
tercepts  of  the  calculated  dashed  curves  rather  well. 
That  is,  in  the  limit  as  Afco,'  =  0  the  value-  of  Ai /p 
and  A t/p  are  the  quantities  k,'+kia  and  ku,,  respectively; 
here  again  the  data  are  consistent  with  the  rate  con¬ 
stant  values. 

2.  Energy  Transfer  in  HF-COi  Mixtures 

Rate  equations  for  the  processes  of  F.qs.  (l)-(8) 
were  solved38  to  provide  expressions  for  the  temporal 
variations  in  HF(r=l)  and  CO2(00°l)  populations,  fn 
contrast  to  the  DF-CO,  case  where  only  the  direct 
coupling  model  (DCM)  was  considered,  an  inter¬ 
mediate  state  model  (ISM)  with  the  alternate  coupling 
scheme  of  Eqs.  (2')  and  (2")  was  considered  in  ad¬ 
dition  to  the  DCM  for  the  IIF-COi  case.  Though  the 
detailed  expressions  are  too  lengthy  to  be  usefully  in¬ 
cluded  here,  one  may  apprehend  that  because  the 
Processes  (2")  are  quite  fast  (iess  than  50  collisions)’7 
compared  to  the  measured  rate  of  coupling  from 
HF(r=l)  to  CO2(00,'l)  states  (about  200  collisions), 
the  temporal  variations  of  the  HF(v=  1)  and  CO,(00°1) 
populations  differ  little  from  thos  predicted  by  the 
DCM  expressions  (11)  and  (12)  of  the  last  section. 
More  precisely,  this  is  true  if  the  inequalities  ky,- f 
F2o»^'m  and  kt\  h 7J(/»/'M  are  satisfied,  where  kK, 
ky,,  and  kn  have  been  defined  in  Eqs.  (2')  and  (2") 
and  kta  and  kyf  are  total  rate  constants  for  deactivation 
of  the  intermediate  CC>2(  10°1)  or  CO2(02°l)  states  by 
ground  state  HF  and  CO-,  respectively.  The  second 
of  these  inequalities  probably  holds,  as  stated  above; 
however,  it  is  not  certain  that  the  first  inequality  is 
satisfied  since  k2n  is  unknown  (kw  cannot  be  much 
larger  than  kta  since  the  HF-C02  laser  i:  known  to  be 
capable  of  reasonably  efficient  performance38) . 

Over  a  very  small  range  of  experimental  parameters, 
i.t.,  when  AVo2«ArnK,  there  are  differences  in  the 
initial  increases  of  CXMOO’l)  population  predicted 
by  the  two  models,  as  Fig.  12  illustrates.  Two  theo¬ 
retical  curves  are  given  in  Fig.  12  which  are  solutions 
to  the  rate  equations  for  the  DCM  and  ISM  assump¬ 
tions,  respectively;  the  curves  have  been  calculated 
with  the  measured  rate  constants  of  this  section.  A 
small  correction  for  the  finite  width  of  the  laser  pulse 
has  been  accounted  for  by  the  introduction  of  a  source 


Eio.  12.  Calculated  initial  rise  in  the  COifOtEl)  population 
(at  low  CO,  concentration)  for  two  different  HF-CO,  coupling 
mechanisms  (see  text)  compared  with  observed  CO,  fluorescence 
at  350°K  from  the  mixture:  HF  » 0.630  torr,  CO, - 0.050  torr, 
and  argon  “10.0  torr.  Calculated  curves  have  been  normalized 
to  the  peak  fluorescence;  the  temporal  profile  of  the  laser  excita¬ 
tion  pulse  is  also  shown. 

term  into  the  rate  equations  having  a  time  profile 
similar  to  that  shown  in  Fig.  12.  In  the  absence  of 
better  information  we  have  assumed  ffejo^io  and 
kioe=kio‘.  The  rate  constant  for  Process  (2")  was 
taken  to  be  ku=  10*  sec-1* torr-1  which  corresponds  to 
about  10  collisions.  The  rise  in  fluorescence  intensity 
predicted  with  the  ISM  exhibits  a  significant  lag  with 
respect  to  both  the  variation  predicted  with  the  DCM 
and  the  experimental  fluorescence  tra;e.  The  fact  that 
the  experimental  data  and  the  prediction  based  on  the 
DCM  are  in  excellent  agreement  does  not  imply  that 
the  ISM  is  invalid  since  we  must  account  for  the 
potential  existence  of  fluorescence  from  the  CO2(10°l) 
or  C02(02°l)  states  superimposed  upon  that  detected 
from  the  CO2(00°l)  state.  Because  the  strengths  of  the 
“hot”  bands  at  4.3  u  [Eq.  (8)  3  are  comparable  to  that 
of  the  fundamental  CO^QOO0!)— >(00P0)J  band,  no 
measurable  differences  in  the  temporal  variations  of  the 
aggregate  4.3  y  fluorescences  are  predicted  by  the  two 
models  under  the  present  experimental  conditions. 
Thus  until  an  experiment  is  performed  with  enough 
fluorescence  intensity  to  permit  the  resolution  of  rota¬ 
tional  components  of  the  CO;[(00#1)— »(0(P0)]  band 
and  the  identification  of  the  fluorescence  (if  any)  from 
intermediate  states,  an  unfortunate  ambiguity  must 
remain  in  the  V— *V  coupling  mechanism  between 
HF  and  C02.  (See  “Second  Note  Added  in  Proof”  at 
end  of  paper.) 

The  fluorescence  data  for  the  HF-CO,  system  are 
adequately  descr.oed  by  the  equations: 

[HF(r=  l)]/[HF3o=exp(— (Jjf),  (11') 

[COi*]/[H  FJo = [<ij/  (<ii  -  W)  ][exp( — b-J)  -  exp(ii|  ()  ], 

(12') 

where 

ih/p  =  k*X  co„ 

<ll/p=  fo»A'|iF  +  A*A'co,+£t»‘-YAr+  (^T,,)-1, 
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Fic.  13.  Observed  HF  fluorescence  decay  rates  for  HI'-COj-Ar 
mixtures  at  350°K  as  a  function  of  the  reduced  mole  fraction  of 
CO-;  Xcoi'“ATco2/(Ahf+Xcoj). 


and 

VP-*ioA,iir»^io,.Vco,+*io,.VAr-l-(#r,)'‘.  (Id') 

Here  k*  is  interpreted  as  the  rate  constant  It.  with  the 
DCM,  or  as  the  rate  constant  kbi  with  the  ISM  as¬ 
sumption.  These  equations  are  consistent  with  the 
neglect  of  k,'  compared  with  k,  lor  the  DCM,  and  also 
apply  for  the  ISM  as  a  consequence  of  the  foregoing 
considerations.  COj*  denotes  those  states  which  give 
rise  to  fluorescence  within  the  pass  band  of  the  COj 
filters.  Since  it  is  expected  that  *h»£m,  then  sufficiently 
late  in  the  fluorescence  decay  period  COj*  will  represent 
only  the  COa(00*l)  state  with  either  the  ISM  or  DCM 
assumptions.  This  circumstance  permits  the  determin¬ 
ation  of  the  rate  constant  k\a. 

The  procedure  for  rate  constant  determination  from 
the  fluorescence  data  for  the  HF-COj  system  was 
similar  to  that  already  discussed  for  the  DF-COi 
case.  The  HF(ti=l)  fluorescence  from  HF-COi 
mixtures  always  exhibited  a  single  exponential  decay. 
Relaxation  times  for  these  data  are  given  in  Fig.  13  as 
a  function  of  the  quantity  Xcc,'*=  (A'co,)/(A'hf4- 
Xcot).  Gas  compositions  consisted  of  10  torr  argon 
with  selected  partial  pressures  for  HF  and  CO*  in  the 
respective  ranges  from  0.04  to  0.41  torr  and  0,04  to 
1.50  torr.  The  solid  line  fit  to  these  data  (Fig.  13) 
gives  a  value  for  the  ATcoi'^l  intercept  of  ( pr)~ ’= 
3.7±0.3X10*  sec_,«torr_1;  according  to  Eq.  (IT)  this 
value  is  the  rate  constant  k*.  We  note  that  the  left- 
hand  intercept  of  the  solid  line  of  Fig.  13  at  Xcoj'“0 
is  (pr)_l=  5.25*0.4X10*  sec-1 -torr-1  which  is  in 
excellent  agreement  with  the  value  of  ku  for  the  self- 
deactivaiion  rate  of  HF  already  discussed.  This  agree¬ 
ment  and  the  linearity  of  the  solid  curve  of  Fig.  13 
is  exactly  the  behavior  predicted  by  Eq.  (IT). 

The  rate  constant  k\0  for  deactivation  of  CO»(00pl) 
by  HF  was  determined  by  examination  of  the  late  time 
behavior  of  the  4 ,3n  fluorescence  which  exhibited 
a  single  exponential  decay.  Relaxation  times  for  the 
COi  fluorescence  data  are  given  in  Figs.  14A  and  14B 
for  mixtures  of  HF  and  COi  and  HF,  COi,  and  argon, 


respectively.  Taken  together  these  data  cover  the 
entire  range  of  -Yhf'  =  (Yhf)/(-Yhf+.Yco,)  from  zero 
to  unity.  The  left-hand  intercept  (Yhf'  =  0)  of  Fig. 
14A  gives  a  value  of  (/>t)-,  =  350±50  sec-1 -torr-1  in 
good  agreement  with  the  known  value15  for  the  rate 
constant  fcio*.  The  straight  solid  lines  of  Figs.  14A  and 
14B  represent  a  reasonable  fit  to  all  of  the  data,  drawn 
to  give  less  weight  to  the  relatively  less  reliable  (/>r)-1 
values  obtained  for  1Yhf'>0.8.  A  value  of  (^r)-l= 
S.fiiO.SXlO4  sec-1  •  torr-1  =  £io  is  obtained  for  Yhf'~  1. 

Figures  ISA,  15B,  and  15C,  taken  under  nearly 
identical  cell  conditions,  illustrate  the  close  correspond¬ 
ence  between  the  fall  of  the  HF  fluorescence  and  the  rise 
of  the  COj  fluorescence  in  agreement  with  the  predic¬ 
tions  of  Eqs.  (IT)  and  (12'). 

The  measured  rate  constants  and  the  corresponding 
values  of  fr  are  summarized  in  Table  I  for  all  of  the 
processes  studied  here. 

IV.  DISCUSSION 

A.  Vibrational  Deactivation  of  HF  and  DF 

Vibrational  relaxation  data  for  HF  at  high  temper¬ 
atures  have  recently  become  available.  Bott  and 
Cohen10  and  Solomon  tl al*  have  performed  independent 


Fio.  14.  Observed  COi  fluorescence  relaxation  times  for  tbe 
latter  stages  of  relaxation  from  HF-COi  mixtures.  A:  No  argon 
present ;  COi  at  either  10  or  100  torr;  and  HF  pressures  from  0.05 
to  1.08  torr.  B:  Argon  pressure- 10  ton:  COi  pressures  from 
0.05  to  2.0  torr;  and  HF  pressures  from  0.05  to  0.90  torr. 
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TAHt.r.  I.  Summary  of  rate  measurements  at  350°K. 


Rate  constant 

pr 

Process 

10*  sec*1 -torr') 

(»i  sec-aim) 

DF(r=l)+DF(e=0)— *2DF(t-0) 

*-2.04:0.2 

0.066 

HF(e-l)  +HF(t>-0)— »2HF(»-0) 

*-5.254:0  30 

0.0  Vj 

DF(r-l)  +C0,{00»0)—  DF(»-0)  +CO,(00»1) 

*-17.25  :2.5 

0.0075 

HF(c—  1)  +C0i(00"0)  — »HF(e-0)+COj(00*l) 

*-3.74:0,3 

0.036 

CO,(OOM )  +DF(t>  -  0)  — CO,(imK))  +DF(r  -  0) 

*-1.94:0.4 

0.069 

CO, (00*1)  +HF(e-0)—  CO,(riw*0)  +HF(»-0) 

*-3.64:0  3 

0.037 

studies  of  the  vibrational  relaxation  of  HF  behind 
incident  shock  waves  in  the  temperature  range  1400- 
4000°K.  These  data  and  the  value  for  HF  deactivation 
at  350°K  given  by  Airey  and  Fried*  are  shown  in  Fig.  16, 
along  with  the  HF  and  DF  deactivation  rates  of  the 
present  work,  plotted  in  the  customary  Landau- 
Teller  fashion  with  pr  as  a  function  of  T~'n.  The  shock 
tube  data  agree  closely  and  both  sets  give  a  linear 
variation  in  pr  with  T~m  for  temperatures  above 
1400°K.  The  values  of  pr  at  350°K  arc  much  lower  than 
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Fig.  15.  Observed  and  calculated  HF  and  COt  fluorescence 
from  an  HF-COj-argon  mixture  at  350°K.  The  calculated  curves 
have  been  normalized  to  the  fluorescence  traces.  A:  COi  fluores¬ 
cent  on  a  long  lime  scale  for  the  mixture:  HF-0.195  torr,  CO,- 
0.360  torr,  argon -9.90  torr.  B:  COj  fluorescence  on  a  short  time 
Kale  for  the  same  conditions  as  in  A.  C:  HF  fluoreKence  for  the 
mixture:  HF— 0.200  torr,  COi-0.360  torr,  argon—  10.0  torr. 


the  high  temperatuic  Jala,  ind;  itive  of  the  effects  of  a 
strong  attractive  HF-HF  n  .action  which  has  a 
predominant  influence  on  ene  gy  transfer  for  kT <£JZ 
where  A£  is  the  attractive  potential  well  depth.  The 
magnitude  of  the  low  temper ’  lure  rates,  and  the  fact 
that  the  DF-DF  deactivation  rate  is  less  than  half  that 
for  HF-HF  deactivation,  is  conclusive  evidence  that 
nearly  all  of  the  initial  vibrational  energy  is  converted 
into  rotational  rather  than  translational  energy  of  the 
collision  pair  in  agreement  with  earlier  results  for  HC1 
and  DC1  presented  by  Chen  and  Moore.,T 

The  qualitative  features  of  the  HF  and  DF  relaxation 
data  can  be  understood  in  terms  of  an  approximate 
classical  calculation  of  vibrational  deactivation  offered 
by  Shin.u*  In  recognition  of  the  importance  of  energy 
transfer  from  vibration  to  rotation  for  molecules 
having  small  moments  of  inertia  so  that  rotational 
velocities  can  exceed  translational  velocities,"  Shin  has 
modeled  HX-HX  (X  =  F,  Cl,  Br,  I)  collisions  in  terms 
of  the  interaction  between  a  rotation-averaged  vibrator 
molecule  (breathing  sphere)  and  a  rigid  rotator 
inolcculi.,5b  The  strong  hydrogen  bonding  interaction 
which  is  expected  to  be  a  predominant  feature  in  the 
encounters  between  HF  (DF)  molecules  and  to  a 
lesser  extent  for  1IC1  and  IIBr  collision  pairs  is  ac¬ 
counted  for  in  an  approximate  manner  by  considering 
the  attractive  interaction  between  two  ideal  dipoles. 
Despite  the  questionable  validity  of  the  model  as  a 
representation  of  the  actual  interaction  at  close  range, 
and  the  several  approximations  necessary  for  a  mathe- 


Temperolure.l’K) 


Fig.  16.  f.andau-Teller  plots  o(  the  sel/deactivation  times  for 
HF  and  DF  at  1  atm  pressure. 
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matically  tractable  analysis,  the  Shin  theory  provides 
results  that  agree  very  well  with  experimental  data  for 
the  II F,  HC1,  DC1,  HBr,  and  III  molecules.13* ,3b 
It  appears  that  the  major  effect  of  the  attractive  forces 
is  confined  to  the  acceleration  produced  by  the  at¬ 
tractive  potential  which  results  in  a  significantly 
enhanced  repulsion  at  close  range.131’  For  this  reason,  a 
rather  inexact  modeling  of  the  attractive  interaction  is 
apparently  sufficient  to  adequately  represent  existing 
experimental  observations. 14 

The  solid  curves  shown  with  the  experimental  data 
of  Fig.  16  are  the  calculated  results  given  by  Shin  in 
Ref.  13(a).  The  calculated  curves  exhibit  the  correct 
qualitative  variation  of  pr  with  temperature,  and  the 
predicted  values  of  fr  are  certainly  of  the  correct 
magnitude.  The  quantitative  differences  between  ex¬ 
periment  and  theory  should  not  be  regarded  too 
seriously  al  present  in  view  of  the  approximate  nature 
of  the  theory  and  the  uinntainty  in  some  of  tne 
physical  parameters  that  en'er  the  calculation,  llott 
and  Cohen  have  indicated  that  a  good  fit  with  the  high 
temperature  shock  tube  data  can  be  obtained  with  a 
somewhat  smaller  value  of  the  Lennard-Jones  param¬ 
eter  than  the  value  <r  =  3  A  used  by  Shin.10 

The  value  for  the  ratio  of  relaxation  times,  ror/tiiF, 
obtained  in  the  present  work  is  2.6  at  350°K;  this  is 
perhaps  significantly  different  from  the  ratio  rDF/riiF= 
7.6  calculated  by  Shin.13*  Unfortunately  no  other  data 
are  presently  available  concerning  vibrational  re¬ 
laxation  in  DF;  high  temperature  shock  tube  data  arc 
needed  for.  a  more  precise  comparison  with  calculated 
values. 

B.  Vibration-Vibration  Energy  Transfer  from 
HF  and  DF  to  COj 

Perhaps  the  most  interesting  observation  of  the 
present  study  was  the  extremely  rapid  rate  of  coupling 
of  vibrational  energy  from  DF  to  COj;  the  measured 
rate,  *=  17.5±2.5X10*  sec_,*torr“l,  corresponds  to 
only  about  40  collisions,  one  of  the  fastest  rates  of  in- 
termolecular  V-*V  transfer  yet  measured.10  The 
magnitude  of  this  rate  is  inexplicable  in  terms  of 
existing  theories  of  V— *V  transfer  because  of  the  large 
over-all  vibrational  resonance  defect  (558  cm-1)  for 
Process  (2).  Because  the  observed  rise  in  CO^Wl) 
fluorescence  exactly  matched  the  decay  in  DF(»=  1) 
fluorescence  as  discussed  in  the  last  section,  and  since 
coupling  by  means  of  an  intermediate  state  in  COj 
is  not  plausible  in  this  case,"  there  can  be  little  doubt 
that  the  direct  coupling  process  of  Eq.  (2)  is  correct. 

Moore  and  co-workers31  have  measured  both  the 
forward  and  backward  rates  for  the  analogous  V—*V 
coupling  between  HC1  and  COj  and  verified  that  the 
direct  coupling  process 

HCl(t-  1,  y')+CO,(00>0)— *HCl(t=0,  J) 

+C0i(00#l)+A£“537  cm-1,  (15) 


is  correct  and  only  requires  about  80  collisions  (k~ 
9.4X104  src_,‘torr-1)  at  300JK.  They  also  determined 
that  the  corresponding  process  foi  V— *V  coupling 
between  HI  and  CO},  with  AZs  =  — 119  era-1,  requires 
about  60  collisions;  &=  12.6X104  sec'Morr"1,  as  de¬ 
termined  from  the  reverse  rate  at  300°K. 

As  mentioned  in  the  last  section,  a  clear  choice  be¬ 
tween  the  alternate  mechanisms  for  V— *V  coupling 
between  HF  and  CO}  [Processes  (2)  and  (2')]  could 
not  be  made  on  the  basis  of  the  present  experimental 
results.  Nevertheless,  the  rate  of  K— *V  transfer  from 
I1F  to  CO2  was  determined  to  be  k~ 3.7X10*  sec-1, 
torr-1  at  350°K  which  corresponds  to  about  200 
collisions.  No  data  [.Vo/e  added  in  proof:  Data  on  these 
systems  have  recently  been  reported;  see  H.  L.  Chen, 
J.  Chem.  Phys.  55,  5557  (1971).]  appear  to  be  avail¬ 
able  for  the  rates  of  coupling  to  CO  *  from  HBr,  DBr, 
DC1,  and  DT;  however,  it  is  of  some  interest  to  consider 
1  lie  process 

Ds(1  =  1,/')tCO}(000)— D}(e=  ),7)+ CChfOO0!) 

+  AE=641  cnr1.  (16) 

The  reported  rate  for  this  proce.s,"  k=  104  sec-,*iorr-1 
at  300°K,  corresponds  to  about  1500  collisions;  the 
relatively  slow  rate  reflects  a  different  mechanism 
for  V— *V  transfer  in  this  case  as  will  be  discussed  in  the 
following. 

For  comparison  with  experiment,  calculations  of 
V—*V  transfer  rates  for  the  HF-CO},  DF-COt,  HC1- 
CO},  HI-COj,  and  Dj-COi  systems  have  been  per¬ 
formed"  within  the  framework  of  the  theory  of  V-*V 
transfer  by  multipole  moment  interactions  discussed 
by  Sharma  and  Brau.33  The  primary  contribution  to 
the  over-all  rate  of  transfer  comes  from  the  first  non¬ 
vanishing  term  in  the  multipole  moment  expansion  in 
each  case;  i.e.,  dipole-dipole  transition  moment 
interactions  in  all  except  the  D*-CO»  case  which  has  a 
rate  determined  by  the  interaction  of  quadrupole  and 
dipole  transition  moments  of  D»  and  CO»,  respectively. 

A  description  of  these  calculations  is  given  in  the 
Appendix.  The  calculated  V—*V  transfer  probabilities 
are  given  in  Column  8  of  Table  II  for  comparison  with 
the  experimental  probabilities  of  Column  9.  The  cal¬ 
culated  energy  transfer  probabilities  are  several  orders 
of  magnitude  too  small  in  every  instance  except  for  the 
D*-COi  case,  which  is  in  excellent  agreement  with  ex¬ 
periment. 

Considering  first  the  HF-CQ»,  DF-COi,  HCl-COt, 
and  HI-COi  systems,  we  note  that  only  certain  pairs 
of  vibration-rotation  transitions  correspond  to  values 
of  u  less  than  15  cm-1,  where  u  is  the  amount  of  energy 
that  must  go  into  translational  motion.  A  given  pair  of 
allowed  vibration-rotation  transitions  can  give  a  large 
contribution  to  the  over-all  probability  of  transfer 
only  when  «<15  cm-1  for  the  dipole-dipole  case, 
provided  that  the  initial  state  populations  are  not  too 
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Table  II.  Comparison  of  experimental  V->V  transfer  probabilities  with  values  calculated  with  multipole  moment  interaction  theory. 


Squared  vibrational  matrix  elements* 


Molecular 

system 

1  (0|<?.(,,l  1)1*  I  (»/!<?,<*>  |0)  1* 

(10-*»esu*-cm*) 

AEb 

(cm'1) 

<*LJ  * 

(A) 

T 

(°K) 

P 

P UiMnt  * 

P *«Ptl 

HF-CO, ' 

HF(1— 4)):97« 

CO  ([00°0)— »(1091)]:  13h 

242 

3.45 

350 

6. 7X10** 

1.7X10"* 

4.9X10"* 

DF-CO, 

DF(l-4):70* 

CO  [(00°0)-*(00°1)]:  950h 

558 

3.45 

350 

0.30 

1.0x10-* 

2.4X10-* 

HC1-CO* 

HCI(1— *0):  45’ 

C0[(00°0)-(00»l)]:950h 

537 

3.60 

300 

0.26 

1.1X10-* 

1.3X10-* 

Hl-COj 

HI(1— ^)) :  0. 161* 

C0.[(00°0)-*(00»1)]:  950“ 

-119 

4.01 

300 

1. 23X10-* 

2  0X1O-* 

1. 76X10-* 

1  <olO»(,,|i)l* 

1  (1  IQ®|0)|* 

(10-*1  esu’-cm*) 

(10-“esu’-cm‘) 

D,-CO, 

D,;  0,(6):  2. 76* 

C0,[(00°0)-*(00°1)]:  950h 

641 

3.67 

300 

... 

4.7X10-4 

6.7X10-*- 

*  The  notation,  of  Ref.  33  is  adopted  here. 

b  Vibrational  energy  difference  between  band  centers. 

"  <a»=*  (<r,+<r,)/2  is  the  Lennard-Jones  collision  diameter,  calculated  from  data  in  J.  O.  Herschfelder,  C.  F.  Curtis,  and  R  B  Bird 
Molecular  Theory  of  Cases  and  Liquids  (Wiley,  New  York,  1954),  pp.  597, 1200.  For  HF  and  DF  a  value  of  *,-3.0  A  has  been  adopted 
[see  Ref.  13(a)], 

d  P,  denotes  the  theoretical  probability  for  exact  resonance  (see  Ref.  34). 

•  Calculations  of  Ref.  32  (see  Appendix). 

' The  Pr°“ss  HF<*=  *)  +CO,(00°0)— *HF(ii  =  0) -f  COi(02°l)  t-AE=342  cm-*  can  be  neglected  by  comparison  owing  to  the  larger 
resonance  defect  and  the  smaller  CO,  matrix  element  for  this  process. 

«  R.  J.  Lovell  and  W.  F.  Herget,  J.  Opt.  Soc.  Am.  52,  1374  (1962). 

11  H.  Statz,  C.  L.  Tang,  and  G.  F.  Roster,  J.  Appl.  Phys.  37,  4278  (1966). 

■  See  Ref.  18. 

i  W.  S.  Benedict,  R.  Herman,  and  G.  E.  Moore,  J.  Chem.  Phys.  26,  1671  (1957) ;  R.  A.  Toth,  R.  H.  Hunt,  and  E.  K.  Plyler,  I.  Mol. 
Spectry.  35,  110  (1970). 

k  W.  Benesch,  J.  Chem.  Phys.  39,  1048  ( 1963). 

1  Calculation  is  based  on  matrix  elements  for  the  0,(5),  0,(6),  and  0,(7)  transitions  in  deuterium  obtained  by  extrapolation  from 
the  results  for  hydrogen  given  by  G.  Karl  and  J.  D.  Poll,  J.  Chem.  Phys.  46,  2944  (1967)  and  T.  C.  James,  J.  Mol.  Spectry.  32,  512 
( 1969) ;  the  value  given  is  for  the  O,  (6)  transition. 

“  Reference  12. 


small.  In  the  HF-COj  calculation,  for  example,  the 
P(6)  transition  in  HF  combines  with  a  large  fraction 
of  the  allowed  C02(OO°O)— »C02(10°1)  vibration-rota¬ 
tion  transitions  ,o  give  values  of  oj<15  cm-1,  However, 
the  population  of  the  J'—  5  rotational  state  is  only 
about  1/13  of  the  total  HF(t=l)  population  at 
350°K;  unfortunately  because  of  the  large  rotational 
spacing  of  HF,  other  possible  P  and  R  branch  transi¬ 
tions  correspond  to  larger  values  of  u  and  thus  the  total 
calculated  probability  of  transfer  is  small.  The  large 
effect  attributable  to  the  fact  that  u  values  greatly 
exceed  15  cm-1  for  the  DF-C02  and  HC1-CO*  systems 
can  be  ascertained  by  comparison  of  Column  8  with 
the  probabilities  of  Column  7  which  have  been  cal¬ 
culated  with  the  totally  unrealistic  assumption  of 
exact  resonance11;  i.e.  w  =  0  for  all  possible  transitions. 

It  is  clear  from  the  results  of  Table  II  that  the  long 
range  multipole  interaction  theory  cannot  explain 
V-*V  transfer  in  the  HF-C02,  DF-CC>2,  HCl-COz, 
and  HI-CO2  systems.  There  are  other  examples  of 
V-*V  transfer  processes  that  do  not  agree  with  existing 
theory,  as  discussed  by  Stephenson  and  Moore.25 
Just  as  in  the  cases  cited  by  these  authors,  conventional 


theories  of  V-*V  transfer  of  the  “SSH  type”35  based 
upon  transitions  caused  primarily  by  short  range  re¬ 
pulsive  forces  are  also  quite  inadequate  in  the  present 
case.  As  is  well  known,  these  theories  predict  transfei 
probabilities  many  orders  of  magnitude  too  low  at  the 
present  temperatures  for  processes  with  resonance 
defects  as  large  as  are  considered  here. 

One  source  of  the  difficulty  in  reconciling  the  present 
experimental  data  with  theoretical  analysis  seems  clear. 
A  successful  theory  must  necessarily  allow  the  major 
portion  of  the  vibrational  energy  defect  between  initial  and 
final  stales  to  be  taken  up  as  rotational  motion  of  the 
collision  partners.  The  rotational  selection  rules  which 
are  a  feature  of  the  long  range  multipole  interaction 
calculation  do  not  permit  this  in  general  when  the  rota¬ 
tional  energy  level  spacings  are  large.  Even  though 
consideration  of  interactions  between  higher  moments 
and  extension  of  the  calculations  to  higher  orders  of 
perturbation  theory  should  provide  some  increase  in 
the  theoretical  values,  it  is  unreasonable  to  expect 
them  to  even  remotely  approach  the  large  experimental 
rates.  We  are  thus  led  to  abandon  the  long  range 
interaction  small-perturbation  approach  as  an  ex- 
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planation  for  the  HF-COi,  DF-COj,  HCl-COi,  and 
HI-COi  data.  Happily,  however,  the  Dj-COj  case  is 
another  story. 

Sltaima  lias  explained  deactivation  of  C()j(01n0)  by 
Hi  and  Di  in  terms  of  a  dipole-quadrupole  interaction.3* 
Kolker37  has  offered  a  similar  calculation  for  the  de¬ 
activation  of  COi(00°l)  by  Hi.  Calculations33  sum¬ 
marized  in  the  Appendix  indicate  that  the  observed 
probability  for  Process  (16)  (see  Table  II)  can  be 
also  satisfactorily  explained  in  terms  of  an  interaction 
between  the  dipole  transition  moment  for  the  asym¬ 
metric  stretch  of  COi  and  a  quadrupole  transition 
moment  of  Di  which  permits  energy  transfer  by  means 
of  the  following  processes: 

D,  (r  =  1 ,  J  =  3)  +  COi  (00"0)  — >D,(r = 0,  J  =  5) 

+CO,(O0ll)-f-AF=  101  cm-1,  (16') 

D:(e=l,7=4)  +  COa(0(yO)— ♦D,(r=0,/=6) 

+  CO,(000l)+A£=-20cm-1,  (16") 

D,(r=  1,7=5)  +  C0i(0(y0)— *D,(r=0,/=  7) 

+CO,(000l)+A£=-141cm-1.  (16'") 

The  temperature  dependence  of  the  calculated  D*-COt 
V-*V  transfer  probability  is  given  in  Fig.  17.  The  lower 
curve  indicates  the  contribution  to  the  total  probability 
of  transfer  made  by  Process  (16").  The  probability  for 
transfer  decreases  markedly  for  temperatures  below 
300°K  because  of  the  reduced  population  in  the  7  =  4 
rotational  level  of  Dj.  At  higher  temperatures  the 
probability  decreases  with  increasing  temperature  as  is 
characteristic  of  weak  attractive  interactions.  The  good 
agreement  with  experiment  at  300°K  is  encouraging 
since  there  are  no  “adjustable  parameters”  in  the 
theory. 

It  should  be  noted  that  the  mechanism  proposed  by 
Moore  et  al.u  for  the  Dj  to  COi  V—*V  process,  i.e., 
Dt(«“l)+C0a(00n0)-^D,(V=0)+C01(01,l)  followed 
by  C01(01'l)+CO,(00D0)-*CO,(01'0)+C0,(00lll)  is 


unlikely  within  the  framework  of  the  multipole  moment 
interaction  theory  because  of  the  strongly  forbidden 
nature  of  the  (.‘()..(00"0)— rCOjfOl'l)  transition.3* 

An  important  conclusion  to  be  drawn  from  the 
foregoing  discussion  can  be  best  illustrated  by  com¬ 
parison  of  the  DF-COj  and  D2-CO2  cases.  The  dif¬ 
ferences  in  the  rotational  velocities  and  energy  level 
spacings  are  not  large  enough  to  account  for  the  sig¬ 
nificantly  greater  V—*V  transfer  rate  for  the  DF-COj 
system  with  simply  the  recognition  that  most  of  the 
vibrational  energy  defect  must  be  converted  to  rota¬ 
tional  motion;  this  conversion  is  clearly  accomplished  by 
different  mechanisms  in  the  two  cases.  The  apparent 
success  of  the  weak  long  range  interaction  hypothesis 
for  the  Dj-CO:  system  and  the  inapplicability  of  this 
assumption  for  the  DF-CO2  system  suggests  that  in  the 
latter  case  a  strong  attractive  interaction  is  involved. 
This  interaction  causes  a  strongly  enhanced  repulsion 
at  short  range  which  results  in  a  violation  of  the  rota¬ 
tional  selection  rules. 

A  search  of  t lie  literature  reveals  several  references  to 
possible  complex  formation  in  HF-COj  and  IiCl-COi 
gas  mixtures.40-*3  Measurements  of  microwave  ab¬ 
sorption  near  the  far  infrared  rotational  lines  of  HC1 
and  HBr  give  evidence  of  sufficiently  strong  attractive 
interactions  between  HCl-COi  and  HBr-COj  collision 
pairs  to  cause  a  violation  of  the  rotational  selection 
rules.40  Similar  data  exist  for  HF-COj  mixtures.41 
Mixtures  of  HC1  and  HBr  with  nonpolar  gases  such  as 
He  or  Hj  do  not  show  similar  effects.40  Values  for  the 
apparent  binding  energies  of  HF-COj  and  IICl-COj 
complexes  have  been  given  as  3.5  and  2.5  kcal/molc, 
respectively.43 

In  view  of  the  above,  and  since  it  is  known  that  the 
COj  molecule  possesses  an  electron  charge  distribution 
strongly  localized  in  the  vicinity  of  the  oxygen  atoms,44 
it  is  probable  that  the  hydrogen  bonded  configurations 

F— H  •  •  •  0=C=0,  F— D  •  •  •  0=0=0 

Cl — H**  *0=0=0,  and  Br — H***O=C=0 

are  influential  in  vibrational  energy  transfer  processes. 
It  is  doubtful,  however,  that  the  hydrogen  bonding 
energies  exceed  2-4  kcal/mole  in  contrast  with  the 


Table  III.  Probability  for  deactivation  of  C0»(00°t). 


Molecular  system 

T 

(’*) 

^ UDtl 

COj-HK 

350 

4.75X10"' 

COj-DF 

350 

2.55X10-' 

COi-Ht 

300 

2.41X10-** 

COrO: 

300 

2. 88X10"*  • 

COrHCl 

300 

5.6X10-**- 

*  Reference  38. 
h  Reference  31. 
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relatively  large  (6.8  kcal/moi  )  fending  in  HF  pairs. 
Since  the  dipole  moments  bar  .<uc  progressively  smaller 
in  the  sequence  1.83,  1  .OS,  0.80  ».»d  0.42  D  for  DF, 
HCi,  HBr,  and  HI,  respectively  the  case  for  the 
existence  of  a  hydrogen  bonded  BI-CO*  configuration 
is  not  quite  as  convincing  in  the  absence  of  experi¬ 
mental  data.  Notice,  however,  that  since  the  vibra¬ 
tional  resonance  defects  for  energy  transfer  from  each 
of  these  molecules  to  CO*  decrease  in  the  sequent 
558,  537,  210,  and  —119  cm-1,  a  compensating  eflfev... 
exists  to  explain  the  large  rate  of  V-*V  transfer  from 
HI  to  CO*. 

C.  Deactivation  of  COi(00°l)  by  HF  and  DF 

Here  again  because  of  the  possibility  of  relatively 
strong  attractive  interactions  we  expect  relatively  large 
probabilities  for  deactivation  in  comparison  with  H* 
and  D*,  as  Table  HI  illustrates.  Because  the  rates  for 
deactivation  for  HF  and  DF  are  both  large  and  since 
the  probability  for  deactivation  by  HF  exceeds  that 
for  DF,  vibration-rotation  transfer  is  indicated; 
however,  conclusions  are  difficult  to  draw  when  the 
COi(wm'O)  states  are  not  defined. 

Preliminary  results  indicate  an  inverse  temperature 
dependence  for  the  deactivation  probability  of 
CCjfOfPl)  by  HF4*;  if  this  result  is  confirmed  with  more 
careful  measurements,  it  would  4end  to  support  the 
hypothesis  that  attractive  influences  are  important  in 
the  deactivation  of  CO*  by  HF  and  DF.  The  deactiva¬ 
tion  rate  of  COi(00°l)  in  collisions  with  HSl  D*,  and 
HD  increases  with  increasing  temperatures1*;  therefore 
a  strong  attractive  interaction  is  not  a  feature  of  the 
intermolecular  potentials  in  these  cases. 

Until  more  data  concerning  deactivation  of  the 
COjfOO9!)  level  by  other  hydrogen  and  deuterium 
halides  become  available,  it  is  difficult  to  generalize, 
However,  the  relatively  small  value  of  the  probability 
for  deactivation  by  HCI  compared  with  that  by  DF 
(see  Table  III)  is  quite  likely  the  effect  of  the  weaker 
intermolecular  attractive  potential  for  the  HCI-CO* 
case  since  the  similarity  of  the  energy  levels  in  HCI  and 
DF  would  argue  against  a  change  in  mechanism  between 
the  two  cases. 

V.  CONCLUSIONS 

(1)  All  three  types  of  vibrational  energy  transfer 
processes  studied  here  have  large  rates  near  room 
temperature  which  can  be  attributed  to  energy  transfer 
to  rotation  under  the  influence  of  a  sizable  attractive 
intermolecular  potential  well  and  enhanced  repulsion 
at  close  range. 

(2)  The  approximations  inherent  in  the  first  order 
perturbation  theory  of  multipole  moment  interactions 
are  invalid  for  treatment  of  the  energy  transfers  to 
COt(00°l)  from  HX(r=l)  and  DX(r=l)  molecules. 

(3)  The  vibrational  energy  transfer  from  D*(r=  1) 


to  CO*(0(P1)  can  be  explained  in  terms  of  a  quad- 
rupole-dipole  interaction  at  long  range  which  causes 
a  significant  transfer  to  rotation  of  orthodeuterium  at 
room  temperature. 

(4)  A  hydrogen  bonding  interaction  between  HF  or 
DF  and  CO*  is  suggested  here  to  account  for  the  rapid 
transfer  rates  for  the  HF-CO*  and  DF-CO*  systems. 
A  similar  mechanism  is  believed  to  occur  in  HX-HX 
deactivation.1111  Since  CO*  has  no  permanent  dipole 
moment,  the  approximation  made  by  Shinu*  in  treat¬ 
ing  the  HX-HX  interaction  as  one  between  ideal 
dipoles  will  not  account  for  the  HX-CO*  results  and 
may  be  questionable  in  the  HX-HX  case  as  well. 
Nevertheless,  Shin’s  results  agree  well  with  experiment 
for  HX  and  DX  deactivation. 

Second  Mote  Added  in  Proof:  \V.  H.  Green  and  J.  K. 
Hancock  have  recently  employed  a  filter  cell  to  deter¬ 
mine  that  the  contribution  of  the  hot  bands  to  the 
total  fluorescence  is  negligible  for  the  conditions  of 
Fig.  12.  We  may  therefore  conclude  from  the  data  of 
Fig.  12  that  the  coupling  of  vibrational  energy'  from 
HF  to  CO*  occurs  by  the  direct  process  of  Eq.  (2). 
We  thank  them  for  their  timely  interest  in  the  resolu¬ 
tion  of  this  important  question. 
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APPENDIX 

The  probability  for  vibrational  energy  transfer  from 
a  molecule  1  to  a  molecule  2  under  the  interaction  of 
the  respective  inultipoles  lt  and  /*  can  be  written 

P(*u,  T)=  Z  fdi)f(ji) [ <Pi,i:(ffu, «,  T)  >„ 

>i  .ir.ii' .it'  I 

+  — ,  r  >«**],  (ad 

°u'  J.  j  J 

where  f(jt)  and  /(/*)  are  the  rotational  distribution 
functions  for  molecules  I  and  2  with  initial  rotational 
quantum  numbers  j\  and  jt.  The  summation  includes 
all  of  the  respective  allowed  choices  of  initial  and  final 
rotational  quantum  numbers  and  jr~*jY  for 

molecules  1  and  2.  The  velocity  averaged  probability 
(Piih(b,  «,  T))*  is  defined  in  Eq.  (13)  of  Ref.  33. 

The  integration  over  impact  parameter  b  of  Eq. 
(Al)  follows  the  approximation  of  J.  Fiutak  and  J.  Van 
Kranendonk  [Can.  J.  Phys.  41,  21  (1963)]  which 
those  authors  have  referred  to  as  "Case  II.”  This  ap¬ 
proximation  is  necessary  because  the  constant  velocity 
trajectory  assumption  becomes  invalid  for  close 
collisions,  (b<ou);  Expression  (Al)  leads  to  a  some¬ 
what  overestimated  value  for  the  over-all  probability 
of  transfer.**-44 
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The  parameters  which  enter  the  computer  calcula¬ 
tion  of  u,  r))„  are  given  in  Table  II;  the 

notation  used  here  is  the  same  as  that  of  Ref.  33. 


•  Supported  by  the  Advanced  Research  Projects  Agency  under 
ONR  contract  N00014-67-A-0077-0006,  by  the  U.S,  Air  Force 
Office  of  Scientific  Research  under  Grant  AFOSR-70-1952,  and 
by  the  National  Aeronautics  and  Space  Administration  under 
Grant  NGL  33-010-064. 

f  NSF  Trainee,  present  address:  Hughes  Research  Laboratories, 
Malibu,  Calif. 

•  D.  J.  Spencer,  H.  Mirels,  T.  A.  Jacobs,  and  R.  W.  F.  Gross, 
Appl.  Phys.  Letters  16,  235  (1970). 

•  T.  A.  Cool,  R.  R.  Stephens,  and  J.  A.  Shirley,  J.  Appl.  Phys. 
41,  4038  (1970). 

C.  Wittig,  J.  C.  Hassler,  and  P.  D.  Coten.  n,  Nature  226, 
645  (1970). 

4  J.  A.  Glaze,  J.  Finzi,  and  W.  F.  Kruplte,  Appl.  Phys.  Letters 
18,  173  (1971). 

:  V..  A.  Tteinzer,  Intern.  J.  Chem.  Kinetics  2,  335  (1970). 

'The  te..n  “chemical  laser”  refers  to  the  fact  that  the  laser 
pumping  requires  conversion  of  the  energy  release  of  chemical 
reactions  into  the  specific  (nonthermal)  vibrational  er  Station 
of  reaction  products  [K..  E.  Shuler,  Appl.  Opt.  Suppl.  2, 1  *65)  ]. 

The  term  “transfer  ctiemical  laser"  distinguishes  those  terns 
for  which  a  vibrational  energy  transfer  process  from  a  product 
molecule  to  another  molecule  (e.g.,  CO,)  is  requisite  to  laser 
operation. 

1 T.  A.  Cool,  T.  J.  Falk,  and  R.  R.  Stephens,  Appl.  Phys.  Letters 
15,  318  (1969). 

•  J.  R.  Airey  and  S.  F.  Fried,  Chem.  Phys.  Letters  8  ,  23  ( 1971) . 

•  W.  C.  Solomon,  J.  A.  Blauer,  F.  C.  Jaye,  and  J.  G.  Hnat, 
Intern.  J.  Chem.  Kinetics  3,  215  (1971). 

>»  J.  F.  Bott  and  N.  Cohen,  J.  Chem.  Phys.  55,  5124  (1971). 

11  R.  S.  Chang,  R.  A.  McFarlane,  and  G.  J.  Wolga,  J.  Chem. 
Phys.  6,  667  (1972). 

"C.  B.  Moore,  Fluorescence,  edited  by  C.  G.  Guilbault  (Dekker, 
New  York,  1967),  pp.  133-198. 

"  (a)  H.  K.  Shin,  Chem.  Phys.  Letters  10,  81  (1971)  ;  (b) 
If.  K.  Shin,  J.  Phys.  Chem.  75,  1079  (1971). 

HR.  R.  Stephens  and  T.  A.  Cooi,  Rev.  Sci.  Instr.  42,  1489 
(1971). 

“  An  rf  discharge  was  used  for  partial  dissociation  of  F,  (instead 
of  SFi  as  described  in  Ref.  14) ;  further  details  are  given  in  Ref.  26. 

44  G.  T.  Armstrong  and  R.  S.  Jessup,  J.  Res.  Natl.  Bur.  Std. 
64 A,  49  (1960). 

11  H.  L.  Chen  and  C.  B.  Moore,  J.  Chem.  Phys.  54,  4072 
(1971). 

"  The  radiative  lifetime  of  the  00,(00*1)  state  is  approximately 
2.4  mser  [n.  Weber,  R.  J.  Holm,  and  S.  S.  IVnncr,  J.  Chem. 
Phys.  20,  1820  (1552)  J.  R.  J.  Lovell  a  ml  W.  F.  Ilvrg.t  Q.  Opt. 
Soc.  Am.  52,  1374  (1962)]  have  given  a  value  for  the  dipole 
derivative  of  HFof  3ft-  1.5X10~"esu  and  a  value  for  the  square 
of  the  transition  dipole  moment  of  |  R,Y -0.97X10-"  esu'-cm1 
for  HF.  The  square  of  the  transition  dipole  moment  for  DF 
may  be  estimated  from  the  relationship  |  {Mir,)lB,/u„ 

if  we  assume  that  Mi  is  the  same  for  both  HF  and  DF;  this 
leads  to  the  estimate  |  Ro,|,"7X10-!’  esu'-cm1  for  DF.  The 
radiative  lifetimes  corresponding  to  the  above  matrix  elements 
are  about  5.5  and  19  msec,  respectively. 

"The  measured  rate  of  V—>V  transfer  from  DF  to  COi  cor¬ 


responds  to  only  about  40  collisions;  the  creation  of  an  inter¬ 
mediate  state,  e.g.,  DF(s“l)+CO,(00°0)-»DF(»-0)-|- 

COi(01M)-t-AE——  73  cm'1,  followed  by  the  process 

C0,(0l'l)  +  C0,(00“0)-<0,(01'0)+C0,(00°l)  seems  unlikely 
tiecause  of  the  smallness  of  the  matrix  element  associated  with 
the  CO, (00*0)—  CO,(01M)  transition. 

»  R.  Oberley,  K.  N.  Rao,  Y.  H.  Hahn,  and  T.  K.  McCubbin, 
J.  Mol.  Spectry.  25,  138  (1963). 

11  An  estimate  of  the  highest  possible  intercept  of  the  ordinate 
of  Fig.  4  for  Xdi— *0  consistent  with  5%  error  estimates  on  the 
100  torr  argon  data  points  leads  to  the  bound:  £».*<100  sec-1* 
torr-1.  The  same  bound  on  the  HF-argon  deactivation  rate, 
4*.* <100  sec-1-torr_1,  was  obtained  from  the  HF  deta  of  Fig  6. 

"  C.  T.  Bowman  and  D.  J.  Seery,  J.  Chem.  Phys.  50,  1504 
(1969). 

"J.  T.  Yardley  and  C.  B.  Moore,  J.  Chem.  Phys.  46,  4491 
(1967). 

“  If  the  observed  relaxation  times  r  exceeded  50  *»sec,  then 
corrections  were  made  to  the  decay  times  for  the  effect  of  the 
2.4  msec  spontaneous  emission  lifetime. 

“  J,  C.  Stephenson  and  C.  B.  Moore,  J.  Chem.  Phys,  52,  2333 
(1970). 

"  R.  R.  Stephens,  Ph.D.  thesis,  Cornell  University. 

n  Near  resonant  processes  of  this  type  in  CO,  are  known  to  be 
fast.  Examples  are  the  processes  CO,(10“0)-l-CO,(00*0)—» 
2CO,(01'0)+A/i--51  cm-1,  and  CO,(02*0)+CO,(00»0)-» 
2CO,  (Ol'O)-f  A/2-  —  49  cm"1  which  require  about  50  collisions 
[C.  K.  Rhodes,  M.  J.  Kelly,  and  A.  Javan,  J.  Chem.  Phys.  48, 
5730  (1968)]. 

"J.  A.  Shirley,  R.  N.  Sileo,  R.  R.  Stephens,  and  T.  A.  Cool, 
“Purely  chemical  laser  operation  in  the  HF,  DF,  HF-CO,  and 
DF-CO,  systems,”  AIAA  Paper  No.  71-27,  AIAA  9th  Aerospace 
Sciences  Meeting,  New’  York,  25-27  January  1971. 

”  C.  B.  Moore,  J.  Chem.  Phys.  43,  2979  (1965). 

"Note  that  the  measured  rate  is  over  4  times  faster  than  the 
value  estimated  indirectly  by  N.  G.  Basov  el  al.,  Appl.  Opt.  10, 
814(1971). 

'*  H.  L.  Chen,  J.  C.  Stephenson,  and  C.  B.  Moore,  Chem.  Phys. 
Letters  2,  593  (1968). 

**T.  A.  Cool  and  J.  C.  Kershenstein,  "Vibrational  Energy 
Transfer  between  D,  and  CO,”  (unpublished) . 

w  R.  D.  Sharma  and  C.  A.  Brau,  J.  Chem.  Phys.  50,  924 
(1969). 

14  J.  C.  Stephenson,  R.  E.  Wood,  and  C.  B.  Moore,  J.  Chem. 
Phys.  48,  4760  (1968). 

u  K.  F.  Herzfeld  and  T.  A.  Litovitz,  Absorption  and  Dispersion 
o)  Ultrasonic  Waves  (Academic,  New  York,  1959),  Chap.  7. 

••  R.  D.  Sharma,  T.  Chem.  Phys.  50,  919  (1969). 

»’H.  J.  Koltcr,  Chem.  Phys.  Letters  10,  498  (1967). 

"J.  Stephenson,  R.  E.  Wood,  and  C.  B.  Moore,  J.  Chem.  Phys. 
54,3097  (1971). 

"S.  S.  Pcnner,  Quantitative  Molecular  Spectroscope  and  Gas 
Iiiniss'vitirs  (Addison -Wesley,  Reading,  Mass.,  1959),  p.  152. 

44  A.  A.  Maryott  anil  G.  Biinbaum,  J.  Chem.  Phys.  47,  3200 
(1967). 

41 M.  R.  Atwood,  H.  Vu,  and  B.  Vodar,  Spectrochim.  Acta  A23, 
553  (1967). 

41 S.  Bratoz  and  M.  L.  Martin,  J.  Chem.  Phys.  42, 1051  (1965). 

41 T.  G.  Burke  and  D.  F.  Smith,  J.  Mol.  Spectry.  3,  381  (1967). 

44  J.  E.  Lennard-Jones,  J.  Chem.  Phys.  20,  1024  (1952). 

41  Unpublished  data  (this  laboratory). 

44  C.  J.  Gray  and  J.  Van  Kranendonk,  Can.  J.  Phys.  44,  2411 
(1966). 


II •  Rapid  Gas  Mixing  Studies 


(D.  Turcotte) 


In  moat  chemical  lasers  the  reactants  are  mixed  and  then  react 
spontaneously.  The  mixing  process  introduces  a  mixing  time  which  may 
or  may  not  be  shorter  than  the  characteristic  tine  for  the  chemical 
reaction.  In  order  to  get  the  expected  performance  from  a  chemical 
laser  the  mixing  time  should  be  of  the  same  order  as  the  chemical  reaction 
time.  The  failure  of  several  chemical  lasers  to  perform  up  to  the  expected 
power  levels  has  been  blamed  in  part  on  Inadequate  mixing. 

Rapid  mixing  may  be  achieved  by  decreasing  the  mixing  scale.  However 
the  design  of  very  small  injectors  is  expensive  and  may  in  fact  be 
impossible  for  large  devices.  For  small  injectors  mixing  is  expected  to  be 
laminar  and  the  laminar  theory  of  mixing  is  well  developed.  One  method  of 
providing  more  rapid  mixing  is  to  introduce  turbulence.  It  is  known  that 
turbulent  mixing  may  be  orders  of  magnitude  more  rapid  than  laminar  mixing. 
Turbulence  may  be  produced  simply  by  using  injectors  which  are  sufficiently 
large  to  produce  turbulent  eddies  or  the  turbulence  may  be  introduced 
artificially  by  using  grids  upstream  of  the  mixing  region.  No  adequate 
theory  exists  for  turbulent  mixing  and  very  few  experiments  have  been  carried 
out  on  turbulent  mixing  at  low  pressures. 

In  order  to  study  mixing  at  low  pressures  a  mixing  wind  tunnel  has 
been  designed  and  built.  The  mixing  tunnel  is  illustrated  in  Figure  1. 

The  plenum  chamber  of  the  wind  tunnel  is  divided  into  two  parts  by  a  splitter 
plate.  The  two  gases  that  are  to  be  mixed  are  metered  into  the  two  plenum 
chambers.  All  experiments  used  helium  and  argon.  The  splitter  plate  ends 


-20- 


/ 


-21- 


with  a  sharp  edge  at  the  beginning  of  the  constant  area  mixing  section.  The 
mixing  section  has  a  3/8"  by  3/4"  cross  section  and  is  4"  long.  The  mixing 
section  exhausts  into  a  plenum  which  is  maintained  at  a  low  pressure  by  a 
vacuum  pump.  Just  upstream  of  the  mixing  region  grids  of  various  size  wires 
could  be  placed  to  induce  turbulence . 

Flow  velocities  were  obtained  using  a  pitot  probe.  The  pitot  probe 
was  also  used  to  obtain  the  variation  in  static  pressure  across  the 
mixing  zone.  Concentration  profiles  were  obtained  using  a  sampling  probe. 
Gas  samples  extracted  from  the  mixing  region  were  analyzed  using  a  thermal 
conductivity  cell.  This  cell  is  illustrated  in  Figure  2.  The  cell  contains 
four  heated,  thin  resistance  wires.  The  sampled  gas  passes  over  two  of 
the  wires  while  a  reference  gas  (in  this  case  argon)  passes  over  the  two 
remaining  wires.  The  four  wires  comprise  a  bridge  circuit  which  measures 
the  difference  in  the  voltage  drop  between  a  wire  in  the  reference  gas  and 
a  wire  in  the  sampled  gas.  This  voltage  difference  AV  i$  related  to  the 
thermal  conductivities  k  of  the  two  gases  by 

4V-cfc-fe) 

where  C  is  a  constant  that  is  evaluated  by  celibration.  With  the  probe 
sampling  pure  helium  the  two  thermal  conductivities  are  known,  AV  is 
measured,  and  C  is  determined.  Using  this  value  of  C,  AV  is  measured 
across  the  mixing  zone  and  the  thermal  conductivity  profile  is  determined. 

The  dependence  of  thermal  conductivity  on  the  relative  ratios  of  helium  and 
argon  is  known1  so  that  the  concentration  profile  is  then  determined.  A 
complete  schematic  diagram  is  shown  in  Figure  3. 

S.  Chapman  and  T.  G.  Cowling.  The  Mathematical  Theory  of  Non-Uniform 

Gases  (Cambridge,  1960). 
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Lamina  r  Measurements 


For  the  constant  velocity  (u) ,  laminar  mixing  of  helium  and  argon 

it  is  expected  that  the  molal  concentration  of  helium,  c„  ,  will  satisfy 

H6 

the  equation 
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(2) 


where  x  is  measured  along  the  channel  and  y  is  measured  across  it  with 
°HeA  tbe  binary  diffusion  coefficient  for  the  mixing  of  helium  and  argon. 
Since  mixing  occurs  only  in  the  central  part  of  the  channel  it  is  appropriate 
to  apply  the  boundary  conditions  c  ->-0  as  /-*■-<*>  an<j  cu  -*•  1  as 
y  +  00 •  The  solution  of  eq.  (1)  that  satisfies  these  boundary  conditions  is 


cHe  ‘  I  [1  +  erf(2  '  5 


HeA 


x 


)  )] 


(3) 


Three  measured  concentration  profiles  are  given  in  Figure  4.  These  were 

obtained  at  x  ■  5.2  cm.  Velocities  of  100  and  200  m/sec  and  pressures  of 

20  and  30  torr  were  considered.  The  measurements  are  compared  with  the 

profile  given  in  eq.  (2)  using  the  accepted  value^  for  D  ., 

HeA 

2 

■  17.8  cm  /sec  for  p  ■  30  torr.  The  agreement  between  theory  and 
experiment  is  reasonably  good. 

The  most  significant  result  of  the  laminar  measurements  is  not 
illustrated  in  the  above  analysis.  The  observed  center  line  for  mixing  is 
given  in  Table  1  for  the  three  cases  considered.  With  equal  velocities 
of  helium  and  argon  about  60%  of  the  channel  is  filled  with  argon.  If 
equal  amounts  of  argon  and  helium  were  metered  into  the  channel,  it  was 
found  that  the  velocity  of  the  helium  would  be  nearly 
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Table  1 

pressure 

torr 

velocity 

m/s 

fraction  of  channel 
filled  with  argon 

PA  '  P 

1  pHe  ' 

30 

100 

0.605 

1.35 

30 

200 

0.575 

1.27 

20 

200 

0.622 

1.22 

twice  as  large  as  the  velocity  of  argon  and  that  about  two-thirds  of  the 
channel  would  be  filled  with  argon.  The  thickness  of  the  mixing  region  in 
this  situation  would  be  about  one-half  that  observed  for  constant  velocity 
mixing. 

In  order  to  investigate  this  peculiar  behavior  the  static  pressure 

was  measured  across  the  mixing  zone.  If  the  variation  in  static  pressure 

across  the  zone  is  a  significant  fraction  of  i  p„  u2  then  an  influence 

He 

of  mixing  on  flow  is  to  be  expected.  In  Table  1  the  measured  values  of 
<PA  “  pHe^*  PHe  u  are  Siven  for  the  three  cases.  It  is  seen  that  this 
ratio  is  in  fact  of  order  unity  so  that  it  is  not  surprising  that  the  mixing 
Induced  pressure  gradient  has  a  significant  effect  on  the  flow. 

This  lateral  pressure  gradient  is  clearly  a  deviation  from  ideal  mixing. 
In  ideal  mixing  the  interaction  forces  between  the  two  gases  are  equal  and 
opposite  as  they  mix  so  that  a  symmetric  mixing  profile  is  obtained  as 
shown  in  Figure  4.  However  the  argon  being  heavier  than  helium  has  a 
greater  inertia.  In  the  rapid  mixing  studies  considered  here  this  inertia 
is  significant  and  the  slowness  sf  the  argon  to  move  into  the  helium 
results  in  the  higher  pressure  in  the  argon  which  is  measured  here.  T»> 
date  a  quantitative  theory  has  not  been  developed. 
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Turbulent  Measurements 

In  order  to  study  the  influence  of  turbulence  on  mixing,  screens 
were  placed  across  the  flow  at  the  entrance  to  the  mixing  section.  The 

results  for  two  screens  are  shown  in  Figure  5.  These  measurements  were 

carried  out  with  p  -  30  torr,  u  -  100  m/s,  and  x  -  5.2  sec.  The  screens 
had  wire  diameters  of  0.012  in.  and  0.018  in.  For  the  finer  screen  the 
concentration  in  helium-rich  mixtures  was  unaffected  but  a  turbulent 
enhancement  of  mixing  is  noted  for  the  argon-rich  mixtures.  This  difference 
can  be  explained  by  the  different  densities  for  argon  and  helium  and  the 
resultant  differences  in  Reynolds  number  based  on  screen  wire  diameter. 

The  Reynolds  number  for  the  screen  in  argon  is  a  factor  of  ten  larger  than 

the  Reynolds  number  for  the  screen  in  helium.  For  the  coarser  screen  no 

additional  mixing  is  observed  in  the  argon  but  the  mixing  in  the  helium 
is  now  as  effective  as  in  the  argon.  Apparently  the  transition  Reynolds 
number  for  the  screen  in  helium  has  now  been  exceeded  and  turbulence 
is  being  produced. 

It  should  be  noted  that  the  pressure  drop  across  the  screen  in  the 
argon  is  much  greater  than  the  pressure  drop  in  the  helium.  This  is 
because  the  pressure  drop  is  proportional  to  pu2  .  Because  of  this 
effect  the  helium  may  be  accelerated  through  the  screen  relative  to  the 
argon  and  a  "jet"  effect  obtained,  thereby  inhibiting  mixirg.  It  is 
necessary  that  turbulence  be  introduced  before  mixing  takes  place  to 
prevent  this  effect. 

In  Figure  6  the  development  of  the  turbulent  mixing  is  shown. 

Mixing  profiles  at  the  three  measuring  stations  are  compared  with  the 
laminar  measurement  at  x  ■  5.2  cm  and  the  laminar  theory.  These 
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results  are  for  the  0.012  in.  screen.  It  should  be  noted  that  the  increase 
in  mixing  observed  in  this  figure  is  in  addition  to  the  spread  of  the  laminar 
mixing  profile. 

In  conclusion  two  important  results  have  been  obtained: 

(1)  Turbulence  can  be  introduced  to  enhance  mixing  at  low  pressures. 

(2)  Transverse  pressure  gradients  heretofore  neglecter  can  strongly 
influence  relative  velocities  of  the  mixing  gas  streams  and 
thereby  inhibit  mixing.  Care  must  be  taken  to  prevent  this 
eir'ect.  If  a  simple  jet  of  a  low  molncular  weight  gas  is 
introduced  into  a  high  molecular  weight  gas  the  transverse 
pressure  effect  will  accelerate  the  low  molecular  weight  jet 
thereby  reducing  the  mixing  predicted  by  the  usual  laminar 
theory  by  a  factor  of  two  to  three. 
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Thermoconductivity  Cell  Circuit 
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III.  Vibrational  Relaxation 


A.  High  Pressure  Relaxation  in  CO2(00rl)  (Wolga) 

1.  HBr  Chemical  Laser  Studies 

The  object  of  this  program  is  to  determine  the  relaxation  time  of 
C02(00#1)  due  to  collisions  with  C02,  N2>  He,  etc.,  in  the  pressure  range 
1-100  atm-,  and  at  temperatures  up  to  6008K.  The  experimental  plan  is  to 
use  laser  induced  fluorescence  from  the  C02(00°1)  state  with  laser  pumping 
by  v  ■  2  -*■  1  transitions  in  the  HBr  chemical  laser.  In  the  process  of 
developing  an  HBr  laser  with  adequate  performance  as  a  pump  source  for 
this  exj  .iraent  extensive  studies  were  carried  out  to  study  beneficial 
additives  if  any,  and  possibly  improved  H  atom  sources.  The  latter  study 
was  motivated  by  the  fact  that  the  chain  propagating  step  in  the 
Hj  +  Br2  chemical  reaction,  Br  +  H2  -*■  HBr  +  H,  is  endothermic  by 
16.7  K  cal/mole,  is  very  slow,  and  does  not  therefore  contribute  to  the 
reaction  within  the  short  time  of  the  chemical  laser  pulse.  A  detailed 
discussion  of  our  studies  of  H  atom  sources  for  pulsed,  electrical  discharge 
initiated  hydrogen  halide  chemical  lasers,  and  of  our  design  and  construction 
of  a  powerful  HBr  chemical  laser  with  substantial  emission  at  4.2  y  was 
included  in  our  Semi  Annual  Report  for  1972  and  will  not  be  repeated  here. 

2.  The  remainder  of  the  past  contract  year  was  devoted  to  completing 
preparations  for  the  measurement  of  the  high  pressure  V-T  relaxation  of  C02(00#1). 
The  HBr  laser  yields  pulses  of  several  microsecond  duration.  Since  the  C02 
fluorescence  will  be  monitored  in  this  experiment  it  fs  necessary  that  the 
exciting  laser  pulse  be  short  compared  to  the  V-T  relaxation  time  of  COj. 

We  plan  to  employ  an  electrooptic  shutter  to  select  a  portion  of  the  HBr  laser 
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pulse,  external  to  the  laser  cavity  thereby  achieving  the  required  short, 
variable  pumping  pulse.  An  electrooptic  shutter  has  been  designed  and 
constructed.  It  uses  an  (110)  oriented,  40  mm  long  single  crystal  slab 
of  GaAs  that  is  electroded  transverse  to  the  optical  path  through  the 
crystal  on  110  faces.  An  appropriate  triggerable,  high  voltage  source 
was  purchased  to  drive  the  GaAs  crystal.  The  GaAs  crystal  operates  in 
conjunction  with  a  Ge  slab  analyzer  to  divert  a  selected  portion  of  the 
output  laser  pulse. 

An  optical  cell  was  designed  and  constructed  for  carrying  out  the 
laser  Induced  fluorescence  measurement.  The  cell  is  capable  of  operating 
from  good  vacuum  conditions  to  an  overpressure  of  100  atmospheres.  Its 
design  is  compatible  with  Installation  in  a  dewar  for  experiments  at 
low  temperatures  and  it  can  be  externally  heated  as  well.  The  path 
length  within  the  cell  is  made  variable  to  accommodate  the  very  short 
absorption  lengths  to  be  expected  at  high  pressures.  A  schematic  drawing 
of  the  cell  is  shown  in  the  figure. 

At  the  time  this  report  was  in  preparation  experiments  were  in  progress. 
B.  Vibrational  Relaxation  and  Chemical  Reactions  (Bauer) 

During  the  past  summer  we  undertook  a  concerted  effort  to  assemble  and 
test  several  computer  programs  for  modeling  chemically  reacting  systems 
which  V-T  and  V-V  processes  are  specifically  included.  While  such  programs 
are  available  in  many  other  laboratories  we  have  not  been  able  to  obtain 
sufficiently  detailed  statements  to  permit  their  direct  use  on  our  IBM  360-65 
computer.  Another  drawback  (besides  the  problem  of  adaptation)  is  the 
length  and  consequent  excessive  cost  of  running  these  computer  programs  in 
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in  our  installation.  During  July-September  a  modified  Runga-Kutta 
integration  scheme  for  multiple,  concurrent,  first  order  non-linear 
differential  equations  was  set  up  and  tested  on  two  reactions: 

[F2  +  NO]  and  [CS2  +  02  (10*  dissociated)].  This  program  incorporates 

(a)  the  Treanor  modifications  for  rapid  convergence  (b)  reverse  steps 
for  all  reactions,  (c)  V-T  relaxation  processes,  and  (d)  the  option 
of  running  the  system  either  isothermally  or  adiabatically,  at  constant 
pressure,  to  simulate  rapid  flow  or  shock  conditions.  In  addition,  a 
program  developed  by  Professor  E.  R.  Fisher  (Wayne  State  University) 
has  been  adopted  for  our  computer.  This  includes  subroutines  for  calculating 
V-V  transition  probabilities,  and  their  incorporation  into  the  kinetic 
scheme.  More  important,  to  permit  integration  cf  "stiff"  equations 
(when  fast  rates  of  destruction  approach  fast  rates  of  production,  such 
that  the  system  enters  a  steady  state  condition)  the  program  switches 
from  integration  of  differential  equations  to  the  solution  of  algebraic 
equations  (Keneshea  technique) .  This  program  is  currently  being  tested 
on  an  artificial  model,  the  dissociation  of  hydrogen  at  room  temperature 
initiated  by  exciting  1*  of  the  molecules  to  the  v  -  1  state. 


IV.  Chemical  Laser  Studies  and  Rate  Determination 


A.  EPR  Study  of  the  NO  +  F,  Reaction  (G.  J.  Wolga) 


The  chemical  decomposition  of  F2  by  reaction  with  NO  has  been 
extensively  employed  in  the  purely  chemical,  DF/HF  -  C02,  transfer 
studies  by  Cool,  Falk  and  others.  Until  now,  the  rates  for  the 


primary  reaction  NO  ^*NOF  +  F  and  the  secondary  reaction 
F  +  NO  A-NOF  have  been  taken  from  the  work  of  Rapp  &  Johnson^ 
who  used  visible  emission  studies  of  light  emitted  from  a  dilute 
diffusion  flame  together  with  extensive  analysis  to  deduce  these 


rates.  These  authors  stated  that  the  interpretation  of  their  results 
was  somewhat  ambiguous.  Since  both  NO  and  F  are  paramagnetic,  their 
concentrations  may  be  followed  quantitatively  in  a  flow  system  using  EPR 
absorption  measurements.  We  have  done  this  by  fitting  a  fast-flow  mixing 
system  to  an  x-band  EPR  spectrometer  and  thereby  measuring  the  growth 
of  F  concentration  and  the  decay  of  NO  concentration  as  a  function  of 
reaction  time.  The  analysis  of  our  results  is  entirely  straightforward  since 
the  experiment  is  a  simple  one.  A  calibrated  flow  of  NO  is  established 
upstream  of  the  microwave  cavity.  A  calibrated  flow  of  F2  is  Injected 
into  the  NO  stream  and  the  F  and  NO  concentrations  are  monitored 
as  a  function  of  position  downstream  of  the  injection  position.  With  a 
constant  flow  speed  a  conversion  of  downstream  position  to  time  of 
reaction  is  easily  made. 

The  rates  we  have  determined  in  this  study  are: 


(1) 


Donald  Rapp  and  Harold  S.  Johnson,  J.  Chem.  Phys.  33,  No.  3,  695 
(Sept.,  1960).  — 
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o  3 

k^  ■  2.6  x  10'  cm  /mole  sec 
kj  m  1 • 1  x  10  cm  /mole^  sec 

The  details  concerning  the  experiment  and  the  rate  determinations  are 
discussed  more  fully  In  the  accompanying  manuscript.  It  Is  clear  from  our 
work  that  mathematical  modeling  studies  of  the  DF  -  C02  laser  utilizing 
the  NO  +  F2  ■>  NOF  +  F  reaction  as  a  source  of  F  atoms  must  use  the 
rates  determined  in  this  work  rather  than  those  of  Rapp  &  Johnson. 
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ABSTRACT 


The  reactions  of  nitric  oxide  with  fluorine  molecules  and  fluorine 


atoms  were  studied  using  EPR  detection.  The  two  rate  constants  determined 

^1  9  3 

are:  the  reaction  NO  +  F«  — FNO  +  F,  k,  *  2.6  x  10  cm  /mole  sec;  and 

c  k 

for  the  reaction  NO  +  F  +  M  — ^  FNO  +  M,  kg  a  1.1  x  10^7  cm®/mole^  sec. 
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REACTION  KINETICS  OF  NO+F2  +  FNO+F  AND 
F+NO+M  -*•  FNO+M  USING  EPR  DETECTION 
by  • 

E.  S.  Mooberry  and  G.  J.  Wolga 
Laboratory  of  Plasma  Studies 
and 

School  of  Electrical  Engineering 
Cornell  University,  Ithaca,  New  York 


Because  of  the  interest  in  the  reaction 
^1 

F2  +  NO  V  FNO  +  F  (1 ) 

as  a  fluorine  atom  source  for  chemical  lasers1  the  rate  constants  for  reac¬ 
tion  (1)  and  also  the  succeeding  step 

k2 

F  +  NO+M  FNO  +  M  (2) 

were  measured.  The  reactions  above  have  been  studied  previously  by  Rapp  and 
2 

Johnston  using  the  dilute  diffusion  flame  method.  At  298°K  the  rate  con- 

in 

stant  for  reaction  (1)  was  determined  to  b*  4.9  x  10  cm  /mole  sec.  In  an 

effort  to  check  this  value  a  study  by  EPR  detection  was  undertaken. 

The  experimental  setup  was  similar  in  design  to  that  of  Westenberg  and 

Deh'aas'  utilizing  a  fixed  EPR  cavity  and  movable  injector.  Reaction  (1) 

was  studied  by  measuring  the  decay  of  [NO]  with  time  under  the  condition 

[NO] <  [Fg]  < [He] . Good  logarithmic  decay  plots  of  [NO]  were  found.  See 

Fig.  1.  Table  I  gives  the  data  for  several  determinations  of  k-j .  The 

fluorine  atom  concentration  was  also  monitored  by  EPR  detection  and  found  to 

Increase  linearly  with  time  for  reaction  distances  up  to  ^0  cm.  The  [F] 

was  found  to  correspond  to  the  loss  of  [NO]  down  the  reaction  tube  so  that 

-1- 

i 
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for  small  distances  the  effect  of  reaction  (2)  was  small,  see  Fig.  2. 
Reaction  (2)  was  studied  by  discharging  [Fg]  in  a  sidearm  and  measuring  the 
F  atom  decay  as  a  function  of  distance  down  the  reaction  tube  under  the 
flow  conditions:  [Fg] ,  [F]  < [NO] <  [He].  Table  II  gives  the  operating  condi¬ 
tions  and  the  values  for  kg  that  were  calculated  from  the  logarithmic  decay 
of  F  atoms  as  plotted  in  Fig.  3. 

As  reaction  (2)  is  55.4  kcal .  exothermic  it  is  assumed  to  be  a  three 

body  reaction.  The  rate  constant  for  the  three  body  reaction  cr  Cl  atoms 

16  6  2  4 

with  NO  has  a  similar  rate  in  Helium  of  3  x  10  cm  /mole  sec.  The 
reactions  of  0  and  H  atoms  with  NO  are  three  body  and  have  rates  of  the  ' 
order  of  1016  also5,6. 

A  third  reaction  must  be  considered  and  that  is  the  three  body  recombin¬ 
ation  of  fluorine  atoms: 

F  +  F  +  M+Fg  +  M  (3) 

where  M  is  He  or  the  wall.  This  reaction  was  found  to  be  negligible  under 
the  present  experimental  conditions.  The  microwave  discharge 
cavity  was  placed  at  several  positions  along  the  reaction  tube  while 
flowing  a  mixture  of  Fg  with  excess  helium  at  a  pressure  of  ca_.  1  torr.  No 
measureable  difference  in  the  EPR  fluorine  atom  signal  height  was  found  at 
different  distances  from  the  cavity.  This  contrasts  with  the  fast  recom¬ 
bination  found  for  chlorine  atoms  with  uncoated  quartz  tubing. 

The  rate  constants  for  the  reaction  of  fluorine  with  nitric  oxide  have 

g 

been  measured  in  the  present  work.  The  value  of  k-|  was  found  to  be  2.6  x  10 
cc/mole  sec  which  is  more  than  an  order  of  magnitude  lower  than  the  pre¬ 
vious  result.  The  present  value  of  k^  was  determined  by  loss  of  NO  in  the 
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reaction  and  by  the  production  of  F  atoms  and  therefore  should  be  a 
reliable  Indication  of  the  reaction  rate.  On  the  other  hand  the  present 
determination  of  kg  is  somewhat  more  uncertain.  The  measured  value 
determined  from  the  direct  reaction  of  F  with  NO  when  used  in  the  computer 
analysis  of  the  Fg  +  NO  reaction  data  leads  to  a  prediction  of  about  half 
the  fluorine  atom  production  actually  observed.  The  value  of  8  x  109  cc/mole 
sec  seems  reasonable  on  comparison  with  the  other  NO  reactions  mentioned 
previously.  Using  these  considerations  the  value  of  8  x  109  cc/mole  sec 
should  probably  be  considered  only  as  an  upper  limit.  Further  study  will 
be  undertaken  to  clarify  the  second  reaction  rate. 
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TABLE  I 

Summary  of  Measurements  of  the  Rate  Coefficient 
k1  for  F2  +  NO  FNO  +  F  298K 


Pressure 

mm 

Velocity 

cm/sec 

molei/cm^ 

cm'Vrliole  sec 

0.60 

1190 

5.9  x  10"9 

'  3.0  x  109 

1.00 

1770 

6.5 

2.6 

1.00 

1840 

8.2 

3.0 

0.83 

1460 

10.6 

2.0 

ave  =  2.6  +  0.5  x  109 
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TABLE  II 

Summary  of  Measurements  of  the  Rate  Coefficient 
for  F  +  NO  +  M  FNO  +  M  298K 


Pressure 

mm 

Velocity 

cm/sec 

[NO]  , 

moles/cm 

3  kJM] 
cm  /mole  sec 

CM]  3 

moles/cm 

cm8/Sol e^sec 

1.0 

1950 

4.4  x  10"9 

8.1  x'109 

5.9  x  10‘8 

1.4  x  1017 

1.3 

2520 

5.9 

8.4 

7.6  x  10"8 

1.1  x  1017 

1.3 

2600 

8.2 

7.9 

7.6  x  10"8 

1.0  x  1017 

ave  ■  1.1  x  10 
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FIGURE  CAPTIONS 

Figure  1.  Plots  of  1n[N0]/[N0]Q  vs.  time.  (1)  Velocity  =  1770  cm/sec, 

P  ■  1.0  torr,  [N0]Q  =  3.2  x  10"9  mole/cm3,  [F£3  =  6.5  x  10-9 
mole/cm3.  (2)  Velocity  *  1840  cm/sec,  P  =  1.0  torr,  [N0]Q  * 
3.2  x  10“9  mole/cm3,  [Fg]  =  8.2  x  10"9  mole/cm3. 

Figure  2.  Plot  of  [F]  vs.  time.  Velocity  ■  1930  cm/sec,  P  =  .92  torr, 
[Fg]  =  5.9  x  10~9  mole/cm3,  [N0]Q  *  3.0  x  10”9  mole/cm3. 
Absolute  F  atom  concentration  determined  by  comparison  of  line 
"d"  of  atomic  fluorine  with  line  "E"  of  molecular  oxygen^. 

Figure  3.  Plot  of  [F]  vs.  time  for  the  reaction  F  +  NO  -*•  FN0.  Velocity  » 
1950  cm/sec,  P  *  1.0  torr,  [NO]  c  4.4  x  10~9  mole/cm3,  [F]  = 
1.35  x  10  mole/cm  assuming  complete  dissociation  of  F?  by 


microwave  discharge. 


u-  0.3 


0.2 


B.  Vibrational  Excitation  In  the  CS^-O^  Laser  Reaction  (S.  H.  Bauer) 

The  major  portion  of  our  study  of  a  pulse  initiated  chemical  CO^ 
laser,  based  on  the  CS2  +  reaction,  has  been  completed  and  the  results 
summarized  in  the  appended  manuscript.  This  has  been  somewhat  revised  and 
submitted  for  publication  in  the  Journal  of  Physical  Chemistry.  The  two 
salient  features  derived  in  this  investigation  are:  (a)  Deduced  the  nascent 
CQ(v)  vibrational  distribution,  as  described  below,  and  concluded  that  in 
a  pulse  initiated  system,  such  as  developed  here,  there  appear  to  be  two 
mechanisms  which  generate  C0^V\  one  identified  with  the  reaction 
CS  +  0  -»•  C0^v^  +  S,  and  another  as  yet  undetermined  process;  (b)  Established 
the  fact  that  the  initiating  electrical  discharge  produces  vibrationally 
excited  CS^  [w  -  0,1,2, 3,4] ,  which  are  long  lived  species. 

In  our  investigations  lasing  was  initiated  by  a  fractional  ys  pulse 

discharge  in  flowing  CS2  +  O2  +  He  mixtures ,  maintained  at  about 
3  torr.  Two  types  of  cavities  were  used;  one  with  two  spherical  mirrors 
and  the  other  with  a  grating  for  wavelength  selection.  In  the  first 
configuration  we  determined  the  optimum  discharge  voltage  to  maximize  laser 
power.  When  the  impressed  voltage  was  higher  than  the  optimum,  CS  and  0 
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were  directly  involved  in  the  lasing  reaction.  When  the  discharge  voltage 
was  lower,  a  chain  reaction  was  the  principal  mechanism  for  lasing.  In 
the  tuned  cavity,  the  measured  t^s  (intervals  between  termination  of  the 
discharge  and  initiation  of  lasing)  followed  a  characteristic  sequence; 
the  minimum  delay  (=  10  us)  occurred  for  the  12  -+■  11  transition,  with 
longer  times  for  higher  and  low  vibrational  states.  The  recorded  delay 
sequence  was  explained  on  the  basis  of  the  initially  generated  distribution 
of  vibrationally  excited  CO.  The  relative  populations  were  also  estimated 
from  chemiluminescence  intensities  recorded  under  conditions  identical  to 
laser  operation  but  in  the  absence  of  mirrors;  these  are:  ^ ^ / . . . . » 

0.13/0.30/0.49/0.72/1.00/1.03/0.90/0.80/0.71/0.62/0.65/0.61/0.57/0.60/0.91/ 

1.20/2.10,  in  general  agreement  with  the  laser  delay  measurements.  While 
this  distribution  is  close  to  that  reported  from  other  laboratories  for 
v  >  8,  there  is  a  significant  differenct  for  the  v  £  7  range,  in  that  our 
data  show  that  substantial  fractions  are  also  generated  in  the  low  v  states. 

To  date  we  have  demonstrated  that  the  rates  of  production  of  total  CO 
and  that  of  destruction  of  total  CS  are  approximately  exponential  and  have 
about  the  same  rate  constant  (Fig.  a).  Our  current  efforts  are  directed 
at  obtaining  more  data  on  the  CS£  +  system  via  time  dependent  UV  emission 
and  UV  absorption  spectra,  to  ascertain  the  fate  of  CS^  species,  and  to 
exploit  the  computer  modeling  programs,  described  in  IIIB,  for  checking  the 
proposed  mechanisms. 

In  a  preliminary  communication  from  this  laboratory  [M.  C.  Lin  and  S.  H. 
Bauer,  Chemical  Physics  Letters,  T_,  223  (1970)]  we  stated  that  mixtures  of 
C3°2  +  °2  Produce  inverted  CO^  populations,  and  do  lase  when  subjected  to  a 
pulse  discharge.  We  have  now  initiated  a  detailed  study  of  this  system, 


/ 


analogous  to  that  concluded  for  the  CS2  +  02  pair.  In  the  first  stage  of 
the  program  shall  to  develop  a  synthetic  procedure  for  generating  abundant 
quantities  of  in  pure  form  and  at  acceptable  yields.  Concurrently 
we  shall  test  for  lasing  under  cw  conditions  by  utilizing  Professor  McFarlane' 
1  Kw  microwave  generator  for  producing  copious  levels  of  dissociated  oxygen. 
Our  primary  objective  is  to  measure  the  nascent  distribution  of  CO^ 
produced  in  +  02  +  He  mixtures  when  subjected  to  a  pulse  discharge. 
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LASING  ACTION  AND  THE 

RELATIVE  POPULATIONS  OF  VEBRATIONALLY  EXCITED  CO  PRODUCED 
IN  PULSE-DISCHARGED  CS  +  O.  +  He  MIXTURES 

tJ  6 

$ 

S.  Tsuchiya  ,  N.  Nielsen  and  S.  H.  Bauer 
Department  of  Chemistry,  Cornell  University,  Ithaca,  New  York  14850 

ABSTRACT 

Lasing  was  initiated  by  a  fractional  its  pulse  discharge  (**15  Kv)  in  flowing 
CS  +  O  +  He  mixtures,  maintained  at  about  3  torr.  Two  types  of  cavities  were 

A  Z 

tested;  one  with  two  spherical  mirrors  and  the  other  with  a  grating  for  wavelength 
selection.  In  the  first  configuration  we  determined  the  optimum  discharge  voltage  to 
maximize  laser  power.  When  the  impressed  voltage  was  higher  than  the  optimum,  CS 
and  O  generated  by  fragmentation  of  the  CS  and  O  were  directly  involved  in  the  lasing 

A  2 

reaction.  When  the  discharge  voltage  was  lower,  a  chain  reaction  was  the  principal 

mechanism  for  lasing.  In  the  tuned  cavity,  the  measured  t.'s  (intervals  between 

a 

termination  of  the  discharge  and  initiation  of  lasing)  followed  a  characteristic  sequence; 

the  minimum  delay  («10  /xa)  occurred  for  the  12  —  11  transition,  with  longer  times  for 

higher  and  lower  vibrational  states.  The  recorded  delay  sequence  was  explained  on 

the  basis  of  the  initially  generated  distribution  of  vibrationally  excited  CO.  The  relative 

populations  were  also  estimated  from  chemiluminescence  intensities  recorded  under 

conditions  identical  to  laser  operation  but  in  the  absence  of  mirrors;  these  are: 

Q  /..../«>,  =  0.13/0.30/0.49/0.72/1.00/1.03/0.90/0.80/0.71/0.62/0.65/0.61/0.57/ 
17  1 

0.60/0.91/1.20/2.10,  in  general  agreement  with  the  laser  delay  measurements.  While 
*  Permanent  Address:  Department  of  Pure  and  Applied  Science,  University  of  Tokyo 
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this  distribution  is  close  to  that  reported  from  other  laboratories  for  v  >  8, 

there  is  a  significant  difference  for  the  v<  7  range,  in  that  our  data  show  that 
substantial  fractions  are  also  generated  in  the  low  v  states. 
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INTRODUCTION 


3 


The  reaction  between  carbon  disulfide  and  oxygen  which,  under  suitable  conditions, 
generates  inverted  populations  of  vibrationally  excited  CO,  has  many  intriguing  aspects. 
While  this  process  has  been  investigated  assiduously  during  the  past  five  years, and  a 
concensus  has  been  reached  regarding  the  two  principal  reactions: 

CS  v  SO  (1) 

(v) 

CO  '  S  (2) 

as  yet  there  has  been  no  complete  identification  of  all  the  steps  which  do  occur,  nor  has 

any  detailed  study  been  made  of  the  lifetimes  and  concentrations  of  intermediates  in  the 

chain  which  must  be  present  in  the  flane  laser Only  a  few  quantitative  measurements 

(2-5) 

of  the  reaction  rates  have  been  published'  ’ .  For  computer  modeling  of  CS2~02  lasers 
it  is  essential  to  establish  the  partition  of  vibrational  energy  among  the  wide  range  of 
states  as  initially  produced  in  reaction  (2).  We  have  no  knowledge  of  any  proposed  3- 
center  potential  energy  surface  for  the  S/C/O  system  nor  of  trajectory  calculations  on 
such  a  surface  which  would  provide  from  a  theoretical  analysis  some  inkling  of  how  the 
exothermicity  is  distributed  in  the  product  species. 

In  this  manuscript  we  report  on  the  deduction  of  the  vibrational  population  distribution 
for  the  nascent  species  from  measurements  of  delay  times  for  lasing  (t),  at  specified 
frequencies,  subsequent  to  initiation  of  the  reaction  by  a  pulsed  electrical  discharge. 

We  have  also  measured  relative  chemiluminescence  intensities  under  conditions  identical 
to  the  above  but  in  the  absence  of  cavity  mirrors.  In  the  meantime  the  reports  by 
Hancock  and  Smith^‘6^  provided  complimentary  data  on  this  distribution,  based  on 


CS2  +  °  - * 

k 

CS  +  O  — 


! 

i 
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their  measured  chemiluminescence  intensities  in  a  rapid  flow  apparatus,  in  which  oxygen 
a 

atoms  produced  by/microwave  discharge  were  mixed  with  carbon  disulfide.  Dawson  and 
Tam'  '  suggested  corrections  to  Hancock  and  Smith’s  initial  distribution.  The  results 
of  these  investigations  are  compared  with  ours.  Also,  we  made  preliminary  measurements 
on  CS  level  populations  via  UV  absorption  spectra.  While  these  data  are  as  yet  qualitative, 
their  implication  regarding  the  mechanism  are  unexpected  and  significant. 

APPARATUS 

I.R. 

The  configuration  used  for  the/  investigations  is  shown  schematically  in 

I 

Fig.  1.  The  laser  tube  is  1  m  long,  2.54  cm  in  diameter.  Two  aluminum  ring 

I 

# 

electrodes  are  placed  at  each  end,  and  between  them  the  charge  from  a  capacitor 

1 

(0.01  pF)  is  released  through  an  ignitron  (WL  -7703),  which  is  triggered  at  a  rate  of 

the 

6  Hz.  The  current  pulse  width  depends  on  the  pressure  and  composition  of/gas,  as  well 
as  the  voltage;  shorter  pulse-widths  and  larger  peak  curients  result  from  higher  discharge 
voltages.  A  typical  combination  has  a  half -width  of  0.3  4 sec,  a  peak  current  of  550A  in 
a  mixture  consisting  of  0.05  torr  CS  ,  0.35  torr  O  ,  and  2.8  tore  He,  for  a  discharge 

Ct 

voltage  of  15  kV. 

Fine  needle  valves  (vernier  settings)  were  used  to  control  the  flow  of  each  gas.  The 

flow  rates  were  calibrated  by  measuring  the  rate  of  decline  of  pressure  in  a  vessel  of 

known  volume.  Most  of  the  experiments  were  done  with  mixtures  of  CS  0.015~  0.05  torr, 

2 

\ 

Og  0.01  ~  0.8  torr,  and  He  2  ~  3  torr,  and  a  linear  flow  velocity  of  b  m/sec. 

first, 

Two  types  of  cavities  were  tested.  In  the/two  gold  coated  spherical  Ge  mirrors  with 
a  radius  of  4  m  were  used.  The  mirrors  were  placed  1.5  m  apart,  and  one  had  an 
uncoated  hole  in  the  center  of  0.75  mm  diameter,  for  coupling  out  the  laser  emission.  The 
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total  power  generated  in  ihe  laser  was  monitored  by  <*  Au  doped  Gt  detector  through  a 

tilted  NaCjP  plate  with  rough  surfaces  used  for  attenuation,  and  an  l.R.  filter  hating  a 

flat  response  between  4. 1  ~  5. 6  \L.  The  ouiput  from  the  detector  passed  through  an  FET 

impedance  reducer,  and  was  displayed  on  an  oscilloscope  (Tekironix  535).  For  dispersed 

wavelength  studies,  an  l.R.  grating  monochromator,  Perkin- Elmer  $SG,  with  a  Bausch- 
grating 

Lomb,  150  C/mm/  blazed  at  6/t,and  a  Au  doped  Ge  detector  (Santa  Barbara  Research 
Lab,  70R)  were  used.  The  monochromator  was  calibrated  by  measuring  several  lines 
of  a  low  pressure  mercury  lamp  in  high  orders.  The  wavelength  response  of  the  mono¬ 
chromator  and  detector  system  was  calibrated  by  recording  the  emission  of  a  Nernst 
glower  lamp  whose  temperature  was  measured  with  an  optical  pyrometer,  assuming  the 
emissivity  was  constant  over  the  0.65u-5/i region.  In  order  to  measure  integrated  power 
the  vertical  signal  output  from  the  oscilloscope  was  passed  through  a  pulse  stretcher  in 
which  the  signal  was  integrated  by  a  0.47/iF  capacitor  through  a  Ge  diode,  and  lead  to  the  input  of  a  lock-in 
amplifier  (PAR  Model  120);  the  reference  signal  was  synchronized  with  the  trigger 
pulse  that  initiated  the  discharge.  The  accuracy  of  response  of  the  lock-  in  amplifier 
output  was  compared  with  the  total  power  obtained  by  integration  of  the  laser  pulse  as 
displayed  on  the  oscilloscope.  The  proportionality  between  the  two  measured  outputs 
was  good  except  for  very  weak  laser  signals. 

In  the  second  arrangement  a  wavelength  controlled  cavity  was  used.  A  grating 
(Bausch  and  Lomb,  300  i/mm,  blazed  at  3.5  ji)  was  set  at  one  end,  replacing  the  totally 
reflecting  spherical  mirror.  The  measuring  system  remained  the  same  as  for  the  con¬ 
ventional  cavity.  Since  the  laser  transitions  of  CO  are  close  to  each  other,  it  was  difficult 
to  isolate  a  single  transition  by  tuning  the  grating.  Most  of  the  tuned  laser  lines  were 


accompanied  by  2  or  3  weak  satellite  lines.  Usually,  an  effort  was  made  to  reduce  the 
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intensity  of  these  satellite  lines  to  less  than  10%  of  the  main  oscillating  lines.  Thus, 
cascading  due  to  stimulated  transitions  had  liiile  effect  on  the  measured  delay  times  of 
the  laser  pulses. 

The  fundamental  CO  emission  was  measured  from  one  end  of  the  laser  tube 
The  Brewster  angle  window  was  replaced  with  a  NaC£ 

window  set  perpendicular  to  the  axis  of  the  tube.  The  arrangement  of  the  discharge 
pulse  was  the  same  as  in  the  former  case.  The  signal  from  the  monochromator  was 

sampled  with  a  Boxcar  integrator,  PAR  lr  0,  for  a  duration  of  2.5  jusec,  using  a  constant 
preset  delay  time  after  the  discharge.  The  monochromator  was  scanned  very  slowly  over  the 

wavelength  region  that  covers  the  CO  fundamental  region,  4.5  to  5.8  fi.  Thus  the 
emission  spectrum  of  CO  at  a  given  delay  time  could  be  obtained.  The  method  used  to 
estimate  the  vibrational  population  from  the  recorded  spectrum  will  be  described  later. 

EXPERIMENTAL  RESULTS 

Laser  with  a  Conventional  Cavity 

Figures  2a  and  2b  are  oscilloscope  records  of  the  laser  pulses.  These  clearly 
show  definite  delays  after  initiation  by  the  discharge.  In  the  Og-rich  case,  the 
laser  power  peaks  very  rapidly  and  then  decreases  in  an  osci  llatory  manner .  In 
contrast,  for  low  O  content,  the  laser  pulse  delays  are  longer,  and  the 
pulses  attain  their  maximum  more  gradually.  Low  discharge  voltages 
and  excess  O.  show  pulse  shapes  very  similar  to  the  case  of  low  O 2 
content.  Thus,  two  types  of  pulse  shapes  were  found:  one  for  high  discharge 
voltage  and  O -rich  mixtures,  and  the  other  for  low  discharge  voltages  and  either 


/ 


-59- 


7 


low  or  high  O  content.  This  suggests  that  two  lasing  mechanisms  are  operative, 
possibly  differing  in  detail;  this  is  supported  by  experiments  on  the  dependence  of  the 
laser  power  on  the  discharge  voltage,  the  measured  delay  for  lasing  onset,  and  the 
laser  spectrum. 

The  oscillations  of  the  laser  pulse  may  be  related  to  radially  propagating  acoustic 

waves.  The  oscillating  period  does  not  depend  on  the  discharge  voltage.  However,  the 

addition  of  He  to  the  system,  while  maintaining  constant  CS  andO  content,  results  in  a 

2  2 

reduction  of  the  period  of  oscillation  as  well  as  its  amplitude.  The  latter  fact  can  be 

explained  by  the  higher  sound  velocity  and  the  more  homogeneous  discharge  when 

additional  He  was  present  in  the  system.  Figure  3  shows  changes  in  the  shape  of  the 

laser  output  pulse  with  increasing  He  pressures. 

The  dependence  of  the  integrated  laser  power  on  the  discharge  voltage  and  O 

2 

concentration  is  shown  in  Figure  4.  There  is  an  optimum  in  the  discharge  voltage  for 
maximizing  the  laser  power.  This  voltage  depends  on  the  O  content  of  the  system; 
the  lesser  O 2  the  higher  the  optimum  voltage.  If  it  is  assumed  that  there  is  an  optimum 
discharge  voltage  for  the  formation  of  CS  from  CS  and  that  the  concentration  of  O  atoms 
increases  with  the  discharge  voltage,  the  above  result  is  understandable. 

For  high  discharge  voltages  the  laser  power  increases  monotonically  with  CS  ,  as  does 
the  optimum  discharge  voltage.  This  is  reasonable,  because  a  larger  amount  of  energy 
is  required  to  decompose  the  CS2  to  CS.  The  overall  effect  of  added  He  is  small. 

With  increasing  He  pressure,  the  optimum  discharge  voltage  for  laser  power  shifts 
slightly  to  the  high  side,  but  the  power  does  not  change.  Apparently,  the  He  atoms  act 
as  decelerators  for  the  electrons.  The  shift  of  the  optimum  discharge  voltage  was 
also  observed  when  the  size  of  the  discharge  capacitor  was  changed.  With  0.01  jxF, 
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a  mixture  of  0.016  torr  CS  +  0.68  torr  0+1.7  torr  Ho  produces  peak  power  at  9. 5  kV: 
with  a  0.0036  pF  capacitor,  the  peak  power  for  the  same  gas  mixture  was  observed  at 
14.5  kV.  Both  cases  have  almost  the  same  electrical  energy:  0.45  and  0.38  J, 
respectively. 

The  delay  times  for  the  initiation  of  laser  pulses  are  determined  by  the  discharge 

voltage  as  well  as  the  CS^  and  concentrations.  As  shown  in  Figure  5,  the  dependence 

is  complex.  This  is  not  unexpected  since  the  nature  of  the  discharge  depends  on  the 

amount  of  CS2  and  02  present.  The  composition  factor  is  apparent  in  the  case  of  low 

voltage  discharges  (6  ~  9  kV);  a  small  increase  of  O,.  or  CS  reduces  the  laser  delay 

times  but  a  large  increase  reverses  the  trend.  While  on  one  hand  higher  pressures 

decelerate  the  electrons,  on  the  other,  they  provide  more  abundant  sources  for  CS  or 

O.  In  the  high  voltage  discharges  the  plasma  electrons  retain  sufficient  energy  to 

dissociate  O  and  CS  so  that  small  changes  in  O  or  CS„  content  do  not  produce 

significant  changes  in  the  plasma  state.  In  this  respect,  some  information  can  be 

obtained  from  measurement  of  the  current  pulses.  For  voltages  higher  than  10  kV, 

the  half-width  is  generally  constant,  at  about  0.3  jisec,  and  is  not  affected  by  small 

changes  in  CS2  or  C>2  content.  However,  for  discharge  voltages  below  10  kV  the 

half-width  is  larger,  and  an  increase  in  CS  or  O  produces  an  increase  in  the  half- 

width.  This  implies  that  the  addition  of  CS„  or  O  to  the  system  reduces  the  discharge 

peak  current,  l.e. ,  the  resistance  of  the  discharge  gas  is  larger.  Therefore,  when  the 

system  has  sufficient  discharge  energy,  the  extent  of  dissociation  is  proportional 

to  the  CS2  and  C>2  content,  and  the  lasing  process  is  accelerated.  From  Figures  , 

it  appears  that  the  inverse  of  the  delay  time  (r)  is  proportional  to  the  O  and  CS 

2  2 

concentrations  when  the  discharge  voltage  is  13  ~  15  kV: 
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i/r  *  [Cs2J  to2i 

Ninety  vibration -rotation  transitions  were  found;  these  were  assigned  to  various 
transitions  from  v  -  14—13  to  v  2—1.  These  spectra  were  taken  for  a  gas  mixture 
consisting  of  0.017  torr  CS  +  0.66  rorr  O  +1.7  torr  He;  the  discharge  voltages  were 
7,9,  and  12  kV.  In  the  7  kV  discharges,  the  recorded  laser  transitions  were  mostly 
P(16)  ~  P(ll)  of  v  -  14—13  to  v  *  4—3  vibrational  transitions.  The  spectrum  initiated 
by  12  kV  discharges  shows  P(19)~  P(14)  of  v  =  13—12  to  v  -  6—5  transitions.  The 
spectrum  produced  with  a  9  kV  discharge  exhibits  characteristics  of  both  the  low  and 
high  voltage  discharges.  Possibly  the  higher  J  transitions  are  favored  for  the  higher 
voltage  discharges  because  the  corresponding  rotational  temperatures  of  the  system 
are  higher.  In  Figure  6,  the  sums  of  intensities  of  the  various  rotational  transitions 
for  each  vibrational  transition  are  plotted  vs  (v— v-1).  Maxima  appear  at  v  -  11—10 
and  v  -  9—8  for  the  high  and  low  voltage  discharges,  respectively.  This  suggests 
that  for  the  higher  voltage  discharges  there  is  a  larger  contribution  to  the  emitted 
intensity  from  chemical  pumping  in  contrast  to  secondary  effects,  such  as  collisional 
or  stimulated  transitions. 

The  lasing  mechanism  in  this  system  is  complex  because  three  factors  strongly 
interact  to  generate  the  observed  population  distribution  among  the  lasing  levels:  the 
many  chemical  reactions  (with  some  rates  dependent,  on  the  states  of  vibrational  excitation 
of  the  reactants),  collisional  relaxation,  and  stimulated  optical  transitions.  Below, 
we  shall  describe  the  result  of  CO  spontaneous  emission  measurements,  from 
which  it  is  evident  that  collisional  relaxation  is  &  dominant  factor  in  determining  the 
intensities  of  lasing  at  the  low  vibrational  transitions.  The  pulse 
shape  for  a  specified  vibration/rotation  transition  is  essentially 
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the  same  as  the  one  for  the  total  power,  and  a  distinct  cascading  phenomenon,  such  as 

(7) 

reported  by  Patel  is  not  seen.  The  sequence  of  delay  times  for  each  vibration- 

rotation  laser  pulse  were  also  measured  (see  below)  using  the  tuned  cavity.  These  are 

significantly  different  from  observations  made  with  a  conventional  cavity. 

In  summary,  two  types  of  voltage  dependent  lasing  behavior  were  fomd.  In  high  voltage 

discharges  the  laser  pulses  are  short,  high  vibrational  and  high  J  transitions  are 

favored,  and  the  laser  power  decreases  with  increasing  discharge  voltage.  In  low 

voltage  discharges  the  delay  times  for  the  laser  pulse  are  long,  as  is  the  duration  of 

lasing,  and  the  power  increases  with  increasing  discharge  voltage.  The  optimum  laser 

power  can  be  obtained  at  an  intermediate  discharge  voltage  that  depends  on  the  CS 

2 

and  Og  content.  The  discharge  voltage  for  optimum  lasing  is  higher  with  increasing 

CS  and  decreasing  O  .  However,  the  dependence  on  CS,  is  more  sensitive  than  on  O  . 

*  ^  2  2 

When  the  O 2  content  is  fairly  high  the  optimum  laser  power  is  almost  directly 
proportional  to  the  CS  content. 

M 

Laser  with  Tuned  Cavity 

In  order  to  eliminate  the  effect  of  stimulated  transitions,  oscillations  for  specific 
rotation-vibration  transitions  were  selected  in  the  r  vity  with  ihe  aid  of  a  grating. 

These  measurements  were  carried  out  mostly  with  15  kV  discharges  because  the  former 
experiments  suggested  that  under  high  voltage  conditions  the  species  produced  by  the 
discharge  contributed  directly  to  lasing 

Typical  oscilloscope  traces  that  show  a  sequence  of  delay  times  for  lasing,  as  these 
depend  on  the  vibrational  transition,  for  the  selected  rotational  transition  P(10),  are 
reproduced  in  Fig.  7  .  Each  transition  starts  with  a  sharp  rise  in  laser  output  at  a 
distinct  delay  time  r  (J).  The  graphs  in  Fig.  8  show  the  dependence  of  r  (J) 

v,v-i  V,V-1 
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on  v— (v-1),  and  indicate  the  effect  of  added  CO  to  the  mixture.  Also  shown  are  the 
estimated  termination  times  for  the  corresponding  laser  pulses.  1*.  ;s  evident  that  the 
addition  of  a  small  amount  of  CO  has  little  effect,  on  r  .(J)  for  the  initiation  of  lasing 
for  v  *  9,  but  produces  progressively  increasing  dela>s,  up  to  about  a  factor  of  two, 
for  the  lower  v's,  down  to  v  =  5.  Depending  on  the  CO  content  the  intensities  for  the 
lower  transitions  diminish.  At  a  partial  pressure  of  CO  half  of  that  of  the  initial  CS 

2 

(0.04  torr)  the  effect  is  marked.  Added  CO  also  reduces  the  duration  of  lasing  for 
transition  from  levels  lower  than  v  =  9  (Fig.  8  )thus  accounting  for  the  decrease  in 
measured  total  integrated  intensities  at  P(9)  for  each  vibrational  transition;  the  effect 
is  marked  for  transitions  from  v$9,  as  shown  in  Fig.  9  .  The  decrement  is  monotonic 
with  the  amount  of  gas  added.  This  effect  will  be  discussed  further  below. 

At  our  operating  pressures  translation  and  rotaiion  should  be  in  equilibrium.  To 
estimate  the  rotational  temperature,  a  small  amount  of  HCj’  was  introduced  to  serve  as 
a  tracer  radiator.  Then  the  relative  emitted  intensities  of  rotational  lines  R(0)  ~  R(6) 
of  the  fundamental  vibrational  band  for  HCf  were  observed.  The  HCf  was  at  a  partial 
pressure  of  0.14  torr,  for  a  total  operating  pressure  of  3.b  torr.  The  deduced 
temperatures  were  440°  and  450°K  at  11.5  and  26.5  fisec,  respectively,  after  a  discharge 
pulse  with  15  kV.  Correction  for  self- absorption  was  made  by  the  method  of 
Ladenburg  and  Levy*01.  These  values  are  lower  limits,  because  the  added  HCf  extracted 
some  part  of  the  discharge  energy.  However,  this  temperature  was  used  for  the 
analysis  of  the  delay  data  and  the  CO  fluorescence  results.  Considering  the  accuracy  of 
data  and  the  approximations  made,  a  temperature  change  of  50°K  does  not  lead  to 
significantly  different  conclusions. 
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Measurement  of  Spontaneous  Emission 

The  emission  spectra  at  the  CO  fundamental  were  recorded 

at  selected  delay  times  after  the  discharge  pulse  to  establish  the  time 
dependence  of  the  population  distribution  for  each  vibrational  level.  The  emission  was 
too  weak  to  permit  measurement  of  individual  rotation -vibration  lines,  and  the  higher 
vibrational  transitions  overlap  appreciably.  The  monochromator  slits  were  set  so 
that  the  equivalent  triangular  slit  function  had  a  width  of  2; 9  cm"1  at  half-height.  With 
this  low  resolution  the  chemiluminescent  spectrum  appears  as  a  smooth  function 
(except  for  noise).  The  spectrum  was  scanned  slowly  (s»  1  jim/hr.). 


DATA  REDUCTION 

A  computer  program  was  written  to  calculate  the  expected  spectra  for  various 
assumed  distributions  for  comparison  with  the  spectra  recorded  at  the  specified  delay 
times.  The  intensity  of  a  vibration-rotation  line  (v1,  J?— v",  JM)  for  a  diatomic 
molecule  is  given  by 


he  to  F  S  A 

■7*T*  vM  .t"  - 


I(v\J'-v",J")  =  - J"  J  J  v’v"  _exp[-B'J'(J'-H)hcAT1 

•  a  Q  - 


CO 


(3) 


v’v" 


rot 


where  j, ^ v,rj»*  is  the  wave  number  of  the  line;  CCv,v„  is  the  wave  number  of  the  band 

origin;  Sj  is  the  rotational  transition  matrix  [ST  =  (J'+l)  for  a  P  branch  and  S  =J' 

d  J 

for  ail  R-branch  line];  F  is  the  vibration-rotation  interaction  factor;  A  is  the 

«  *  y*yM 

spontaneous  emission  coefficient  for  the  vibrational  transition  (v'^v")-  N  is  the 
population  in  the  upper  vibrational  level,  and  is  the  rotational  partition  function  . 

A  rotational  temperature  of  450°  was  assumed.  The  wave  numbers  for  the  lines 
were  calculated  from  the  molecular  constants  given  by  Mantz,  et.al.^;  values  for 
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AvV  and  Fj  were  &yen  by  Fisher (10)  and  by  Toth,  et.al. (11),  respectively.  First,  the 
contribution  from  each  vibrational  transition  to  the  observed  intensity  was  calculated 
when  the  monochromator  was  set  at  the  wave  number  using  a  triangular  slit  function 

U(v\£0.)  =  ^I(v'Jl^v»J’')f(w,w1),  (4) 

J’ 

where  N^,  =  1  and  v"  =  (v’-l).  Since  the  slit  width  is  much  wider  than  the  Doppler  line 
width,  the  spectrometer  transmission  factor  is  the  magnitude  of  the  normalized  slit  function 
at  the  line  frequency.  The  summation  was  taken  over  the  90  P-  and  R-branch  lines. 

The  observed  intensity  record  was  then  represented  by 

(5) 

V* 

The  best  set  of  relative  N^'s  that  reproduces  the  observed  Ifw^'s  [which  were  measured  at 
4  cm  1  intervals  over  the  rarge  2219-1691  cm”3}  was  sought  by  a  least  square  procedure 
for  v'  =  1  ~  17.  The  deduced  sets  of  N^'s  were  graphically  smoothed.  Figure  10  shows 
three  examples  of  the  observed  and  calculated  spectra. 

The  principal  sources  of  inaccuracy  in  the  above  procedure  are  lack  of  a  precise 
rotational  temperature,  absorption  by  water  vapor  in  the  region  of  the  higher  vibrational 
transitions,  and  overlap  of  the  CO  emission  by  that  of  OCS.  The  first  was  readily  checked 
bv  recalculation,  using  400  and  500°K  as  trial  temperatures.  The  changes  in  the  resulting 
distributions  were  less  than  10%.  Regarding  interference  by  water,  the  monochromator 
was  flushed  with  dry  Ng,  while  the  light  path  of  about  30  cm  from  the  laser  tube  to 
the  spectrometer  slit  was  lett  open  to  the  air.  Absorption  by  water  was  clearly  indicated. 
Therefore,  our  deduced  distribution  at  the  high  vibrational  states  v  =  12  ~  17  may  be 
somewhat  underestimated.  With  respect  to  OCS  absorption,  Hancock  and  Smith(5)  in  their 
study  of  the  O  +  CS2  system  found  that  emission  by  OCS  did  overlap  the  CO  fundamentals 
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However,  we  found  that  when  a  cell  10  cm  in  length  filled  with  OCS  at  1  torr  was  placed 
in  front  of  the  entrance  slit  of  the  monochromator,  no  change  in  the  spectrum  was 
observed.  Hancock  and  Smith  suggested  that  at  most  1.5%  of  the  CS  is  converted  to 
OCS  in  O  rich  mixtures.  If  the  excitation  of  OCS  is  due  to  collisions  with  CO,  10 
collisions  require  about  40  /isec.  Since  the  pressure  of  CO  cannot  exceed  the  initial 
pressure  of  CS  ,  which  in  these  experiments  is  less  than  0.05  torr,  CO  starts  at  zero  and 

Lt 

increases  exponentially  following  the  initiation  discharge  pulse,  at  a  rate  of  order 
10  sec  .  Hancock  and  Smith'  ’  showed  that  the  V-V  collision  probability  between 
CO  in  levels  of  v  =  7  ~  4  and  OCS  was  of  the  order  of  10  \  but  the  above  estimate  of 
the  collision  number  includes  collisions  with  CO  in  any  state.  Therefore,  the  contribution 
of  OCS  emission  to  the  spectrum  cannot  be  large  during  the  first  100  /zsec  after  the 
discharge.  Were  excited  OCS  generated  directly,  an  overestimate  of  CO  population  at 
v  =  5  ~  3  may  occur  depending  on  the  rate  of  OCS  formation. 

In  Fig.  11  the  deduced  relative  population  distributions  (normalized  to  unity  at  v  =  13) 
are  shown  as  a  function  of  the  delay  time  after  the  discharge  pulse.  At  a  delay  of  23  /tsec 
a  peak  in  the  population  appears  at  v  =  11  or  12  and  a  minimum  population  occurs  at  v  =  4. 

As  the  delay  time  increases ,  the  peak  shifts  to  lower  v  values  and  the  distribution  becomes 
more  Boltzmann-like.  Figure  11  also  shows  the  initial  distribution  estimated 
by  graphical  linear  extrapolation  to  t  =  0,  together  with  the  recent  results  of  Hancock, 
Morley  and  Smith. '  ' 

To  obtain  the  rate  of  formation  of  CO,  the  observed  relative  populations  in  levels 
v  =  1  to  v  =  17  were  summed;  this  is  shown  in  Fig.  12.  However,  the  population  in  the 
v  =  0  level  caunot  be  observed  by  any  emission  measurement,  and  the  relaxation  process 
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produces  a  significant  amount  of  CO  in  the  ground  slate  during  the  latter  stages. 

Therefore,  we  were  forced  to  estimate  CO(v=0)  by  a  linear  extrapolation  from  the 
higher  levels;  this  is  subject  to  considerable  uncertainty.  The  estimated  total  CO  was 
also  plotted  in  Fig.  12  .  The  carbon  monoxide  accumulates  exponentially  as  a  function 
of  time , 

[CO]  ~  (1  -  e_kt)  (6) 

4  -1 

where  k  =  2.3  x  10  sec  ,  when  the  initial  pressures  of  He,  O  ,  and  CS  are  3. 1,  0.75, 

Ct  M 

and  0.046  torr,  respectively  and  the  discharge  voltage  is  15  kV.  It  is  interesting  to  note 

that  the  rate  of  decrease  of  [CS],  as  measured  via  absorption  in  the  UV,  follows  a  decay 

(13) 

rate  that  closely  parallels  the  growth  of  CO. '  ’  Two  types  of  the  relaxation  processes 

deactivate  the  excited  CO,  V-V  and  V-T  energ>  transfers.  The  former  is  a  near-resonant 

exchange  with  a  small  amount  of  energy  (equal  to  the  anharmonicity  defect)  transferred 

to  the  translational  motion.  This  process  very  efficiently  modifies  ihe  initial  distribution 
(v) 

of  CO  generated  in  the  reaction.  However,  the  V--V  process  does  not  change  the  total 

number  of  vibrational  quanta  produced  by  the  reaction.  In  contrast,  when  V-T  energy 

transfers  occur,  the  total  number  of  quanta  are  reduced.  The  average  number  of 

vibrational  quanta  for  CO  is  estimated  to  be  6.0,  which  corresponds  to  42%  of  the  exo- 

(v) 

thermicity  of  the  pumping  reaction:  CS  +  O  —  CO  +  S.  This  average  number  of 

quanta  is  retained  during  the  first  60  usee  after  the  discharge,  and  then  gradually 

decreases.  Thus,  prior  to  60  u sec,  V-T  relaxation  is  not  as  important  as  V-V  relaxation, 

which  dominates.  The  decrement  rate  for  the  number  of  vibrational  quanta  after  60  u sec 
3  -1 

is  about  4  x  10  sec  ,  which  is  much  larger  than  the  spontaneous  emission  rate  for 
2  -1 

v  =  6—5  (1.61  x  10  sec  ).  V-T  relaxation  of  CO  in  collisions  with  He  or  O^  is  very 
slow  at  room  temperature,  and  therefore,  the  CO  is  efficiently  deactivated  either  at  the 
walls  of  the  tube  or  by  collision  with  the  active  species  such  as  O,  S  or  CS. 
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Interpretation  of  t  ’ s 

As  shown  above,  the  total  CO  concentration  increases  exponentially  to  an  equilibrium 
value.  Before  relaxation  processes  become  dominant  the  population  of  level  v  may  be 
expressed  by 

Ny  -  ),  (7) 

where  N  is  the  final  concentration  of  CO,  and p  =k  /£k  ,  the  fraction  of  production 

v  V  v  V 

channeled  to  that  level.  The  distribution  of  CO  changes  rapidly  due  to  V  -V  relaxation 

as  shown  in  Fig.  11.  For  these  experiments  the  partial  pressure  of  CS  was  0.05  torr, 

and  this  limits  the  amount  of  CO  that  can  be  produced.  An  upper  bound  for  the  collision 

5  -1 

frequency  between  CO  molecules  is  4  x  10  sec  ;  this  is  an  overestimate  since  CO 
increases  gradually  during  the  initial  period  of  the  reaction  and  reaches  its  maximum 
value  only  when  all  the  CS  is  converted.  We  recorded  delay  times  of  less  than  20  jusec, 

a 

except  for  the  lower  vibrational  transitions.  The  latter  are  involved  in  V-V  relaxation, 
because  the  addition  of  a  small  amount  of  CO  affects  significantly  their  delay  times. 

Hence  laser  pulses  from  14—13  to  7—6  can  be  interpreted  in  terms  of  the  initially 
generated  vibrational  distribution.  For  these, collisional  or  radiative  deactivation 
processes  other  than  the  V-V  relaxation  cannot  make  an  appreciable  contribution  since 
the  average  vibration  quanta  of  CO  are  conserved  in  the  system  during  the  first  60  psec 
after  the  discharge. 

To  deduce  population  ratios  for  adjacent  levels  from  laser  delay  data  an  additional 
assumption  was  made,  that  the  loss  in  the  laser  cavity  is  independent  of  a  small 
change  in  the  grating  angle,  which  in  the  range  of  7°.  The  theoretical  gain  for  a  P- 
branch  transition,  (v,J-l)—  (v-l,J)  is: 


/ 
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„  3  4 

8ff  c  A  , 

_ v.v-1 

i 

3kT(2ffkT  /m)2 

x  J[N  B  exp  f-F  (J-l)  —  j  -  N  exp  {-F  ,(J)  ~  }  ]  (8) 

1  v  v  v  kT  v-1  v- 1  1  v-1  kT  '  ' 


where  A  =  transition  matrix  element  for  v  —  (v-1) 
v,  v-1 

s*  molecule  density  of  CO  in  level  v 
m  =  mass  of  CO 

T  =  translational  --  rotational  temperature 

2  2 

F  (J)  =  rotational  term  value  for  vibrational  level  v  =  B  J(J+i)-D  J  (J+l) 

v  v  v 

For  a  sequence  of  transitions  (v,  J-l)  —  (v-l,"J),  with  J  fixed,  the  time  dependent  relative 

gain  after  the  discharge  pulse,  may  be  written  in  terms  Eq.  (7), 


_lr*  llC  llC 

k  (J)  =  const* A  (1  -e  )  .  [<p  B  exp{-F  (J-l)— ;}  -<p  B  exp[-F  ,(J)— }]. 
v,v-l  v, v-1'  1  v  v  1  v  kT  v-1  v-1  1  v-1  kT 


In  the  limiting  case,  as  t  —  «, 


00  K(J  ,  Kn 

a  (J)“const- A  ,[?)  B  exp[-F  (J-l)-~ j  -  ,B  exp[-F  ,(J)  r^}]»  (10) 

v,v-l  v,v-l  v  v  v  kT  v-1  v-1  v-1  kT 

This  is  the  maximum  gain  were  the  chemically  produced  excited  species  not  to  make 

transitions  to  lower  levels.  Equation  (9)  simplifies  to 


«V,V-1(J)  * av,v-l(J)  (1  "  e"kt)*  (9,) 

Laser  oscillations  start  when  the  gain  exceeds  the  cavity  loss.  Therefore,  the  delay 
time  r  ,  (J)  may  be  defined  by 

Tv,v-l(J)  55  "(lA)  fn  11  "  (LAv,v-l(J))1  *  (11) 

where  L  is  the  loss  in  the  cavity  (equal  to  a  at  threshold).  When  a  value  for  k  is 
available  it  is  more  convenient  to  work  with  the  equation. 
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T  ,(J)  «  L ,/oT  .(J) 

V,V*“1  V,V“1 

where  T  „(J)  s  l  -  expf-kr  ,(J)1.  When  relative  values  for  adjacent 's  are 
v,v-l  v,v-l  *v 

available  one  may  compare  T  „  (J)  with  that  of  T  „  (J): 

J  v,v-l  v-1, v-2 

he  he 

T  (J)  A  (<p  Ap  „  )B  expf-F  (J  l)“*}-B  , expf-F  „(J)r^} 

v-1. v-2'  ;  v.v-l  lVyv-r  v  v'  ;kTJ  v-1  vl'  7kTJ _ 

T  /t\  A  he  he 

v,v-l'  ;  v-1, v-2  B  , expf-F  ,  (J-l)r~}  -  (to  n/ip  ,)B  „ expf-F 

’  ’  v-1  1  v-1'  kT;  '*v-2  v-l  v-2  1  v-2  'kTJ 


18 

(12) 


(13) 


Thus,  once  a  population  ratio  for  one  pair  of  adjacent  vibrational  levels  is  specified, 

the  population  ratios  for  other  pairs  may  be  derived. 

Additional  information  on  population  ratios  may  be  derived  from  the  dependence  of 

T  -  (J)  on  J.  From  Eqs.  (10)  and  (12),  the  following  relation  is  obtained  for  a  fixed 
v,v-i 

vibrational  transition  at  different  J’s  ; 


O  ,  B 

1/T  ,(J)«J[1- - -  — — 

v,v-l  to  B 

^v  v 


he 


e*t',Pv-l,J)'Fv(J'1,,i?rj1 


(14) 


Relative  values  of  l/T  ?  (J)  for  P(3)  to  P(22)  of  the  12—11  transition  for  a  range  of 

v,\  -J. 

population  ratios  /<p  at  T  =  450°K  are  plotted  vs  J  m  Fig.  13 .  Superposed  on  these 

X  X  Xu 

4  X 

are  l/T  from  observed  fs  for  k  2 . 06  x  10  sec  .  Although  there  is  considerable 
v,  v-i 

scatter,  on  comparing  the  experimental  points  with  the  calculated  curves  best  agree¬ 
ment  is  indicated  for  (<p  Ap  j  in  the  range  of  0.8  ~  0.9.  This  range  is  consistent  with 
the  value  of  0.88  derived  from  the  CO  spontaneous  emission  measurements. 

To  estimate  population  ratios  for  this  range  of  (<p  /<p  )'s,  Eq.  (13)  was  applied, 

XI  Xu 

assuming  0.80,  0.85,  and  0.90.  The  results  are  shown  in  Table  I.  They  are 
compared  with  the  relative  population  distributions  deduced  from  the  spontaneous 
emission  measurements  obtained  in  this  laboratory  and  those  reported  by  Hancock, 


i 
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Morley,  and  Smith  .  The  assumption  =  0.8  gives  a  population  ratio  for 

which  is  negative;  this  is  physically  unacceptable.  However,  for 
(<Pl/<P12)  =  °-85~  °-9  the  relative  populations  derived  from  the  r  values  are  in  general 
agreement  with  those  from  the  spontaneous  emission  measurements,  except  for  the 
low  v's,  7  ~  5.  As  shown  in  Fig.  8  ,  for  transitions  lower  than  8—7  the 

relatively  long  delay  is  significantly  affected  by  added  CO.  Hence  it  is  not  valid  to 
assume  that  the  population  of  levels  lower  than  v  *  8  are  predominantly  determined  by 
the  chemical  production  step  at  a  time  when  the  laser  begins  to  oscillate  at  these 
transitions. 

We  may  now  compare  the  gain  calculated  from  the  distribution  based  on  the 
spontaneous  emission  measurements  with  the  laser  delay  data  shown  in  Fig.  8  .  Though 
fluorescence  distributions  have  been  recorded  for  only  5  different  times  after  the 
discharge,  and  these  are  not  sufficient  to  permit  unravelling  of  the  entire  history  of  the 
gain  for  each  rotation-vibration  transition,  reasonable  estimates  of  the  gain-time 
relation  can  be  made  using  eq.  8,  since  Nv  is  now  a  known  function  of  time.  The  results 
are  plotted  in  Fig.  14a, b.  The  P(9)  transitions  of  13  -  12  to  5  —  4  have  positive  gains 
which  initially  increase  with  time,  attain  maxima  and  then  decrease  to  negative  values. 
The  gain  for  the  8  —  7  transition  remains  positive  for  the  longest  time,  while  that  for 
9  —  8  is  largest.  The  gains  for  the  4  —  3  and  3  —  2  transitions  are  negative  initially; 
these  increase  to  a  positive  value  before  decreasing  again  to  become  negative.  From 
Figs.  8  and  9,  the  maximum  power  appears  at  10  -  9  and  the  laser  duration  is  longest 
for  8  -  7.  Also,  a  short  lasing  duration  was  found  for  the  higher  transitions  along  with 
short  delay  times,  while  the  4  —  3  and  3  —  2  transitions  have  long  delay  times. 
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Reaction  Mechanism 

The  data  summarized  in  the  preceding  section  suggest  that  the  mechanisms  for  the 
high  and  low  voltage  discharges  differ  to  some  extent.  When  the  discharge  voltage  is 
high  the  electrical  energy  is  sufficient  to  dissociate  most  of  the  CS  and  the  O  i  ,  so  that 

M  Z 

an  inverted  CO  population  and  hence  lasing  occur  directly,  i.e. , 

tit  CS  J  ® 

' '  2  ICS  +  S  -  102  kcal 

(ii)  O  -  O  +  O-  119  kcal 

(iii)  CS  +  O  -  C01V)  +  S  +  85  kcal 

The  recently  redetermined  value  for  the  heal  of  dissociation  of  CS  ,  as  measured  by 

z 

(14) 

Okabe  ,  was  used.  In  contrast,  during  low  voltage  discharges  an  inverted  CO 
population  occurs  subsequent  to  the  initiation  of  a  chain  reaction.  The  principal  steps 
in  the  sequence  are: 

(ii)  02  ^  O  +  O 


(iii) 

CS  +  0 

— 

co(v)  +  s 

+  85  kcal 

(iv) 

o  +  cs2 

— 

CS  +  so 

+  21.0  kcal 

(V) 

s  +  o2 

- 

so  +  0  + 

5.6  kcal 

(Vi) 

so  +  o 

z 

— 

so2  +  0 

+  12 . 5  kcal 

(vii) 

o  +  cs2 

- 

ocs  +  s 

+  54  kcal 

(viii) 

s  +  cs2 

— 

CS  +  S' 
2 

-  2.0  kcal 

(ix) 

CS  +  o 
z 

- 

ocs  +  o 

+  39.0  kcal 

(x)  SO  +  0  +  M  —  SO  +  M  +  131.9  kcal 

z 

of  these  (vi),  (viii) ,  (ix)  and  (x)  are  of  lesser  importance;  the  chain  is  prorogated  via  (v). 
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The  argument  for  two  types  of  initiation  processes  is  supported  by  our  observations  of  CS 

absorption  in  the  ultraviolet.  Subsequent  to  a  high  voltage  discharge  the  concentration  of 

CS  decays  exponentially  approaching  a  constant  value,  while  after  a  low  voltage  discharge 

(13) 

the  concentration  of  CS  remains  almost  constant  for  more  than  one  msec. '  The 
latter  indicates  that  CS  is  a  reaction  intermediate  which  attains  a  stationary  stale  con¬ 
centration,  whereas  in  the  former  case,  the  decay  rate  of  CS  agrees  with  the  observed 
Increasing  rate  of  production  of  CO.  Also,  an  optimum  discharge  voltage  to  generate 
CS  in  [CS  <•  He]  mixtures  was  found,  and  this  corresponds  to  the  observed  optimum 
discharge  voltage  for  laser  power  production.  Apparently,  when  discharge  voltages 
greater  than  the  optimum  are  used  CS  decomposes  to  give  species  other  than  CS.  The 
experimental  results  on  CS  monitoring  in  the  UV  and  computer  calculations  of  the  reaction 

kinetics  will  be  presented  in  a  later  report. 

(15) 

Recently,  Rosenwaks  and  Yatsiv  suggested  that  the  reaction 

(xi)  CS*  +  O  -  COiV)  +  SO  +  91  kcal 

M 

contributes  substantially  to  population  inversion.  Our  observations  of  CS  indicate  that 
a  rapid  reaction  between  CS  and  O  does  not  occur  since  the  concentration  of  CS  (v-0. .  .4) 

b 

remains  constant  long  after  the  end  of  the  CO  formation  period,  even  with  O  present 

2 

in  excess.  That  the  life-time  of  CS  produced  by  the  flash  photolysis  was  not  affected 

by  the  presence  of  O >  was  also  reported  by  Callear^6\  The  low  probability  for 

•  (5) 

reaction  (xi)  was  also  discussed  by  Hancock  and  Smith  . 

Comparison  with  Other  Experiments 

The  only  published  data  on  the  delay  sequence  for  CO  lasing  are  those  of  Gregg  and 

(17) 

Thomas  .  They  observed  laser  pulses  in  flash  photolyzed  CS  -0  mixtures  using  a 
cavity  with  a  grating  and  a  rotating  mirror  and  found  two  minima  in  delay  times,  those 
associated  with  13—12  and  10—9  transitions.  This  lead  them  to  the  conclusion  that  the 
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(v) 

CO  was  produced  chemically  as  well  as  by  collisions  with  electronically  excited  SO  . 

2 

The  latter  part  of  their  conclusion  is  in  question.  Hancock  and  Sinith(1S)  found  only  a 

weak  CO  emission  associated  with  the  slow  recombination  reaction:  SO  +  O  —  SO* 

2 

Also,  our  observation  of  the  correlation  between  the  CS  and  CO  concentrations  shows 
that  the  main  source  of  CO  is  the  reaction  between  CS  and  O.  However,  Gregg  and  Thomas' 
T's  are  roughly  in  agreement  with  our  results,  except  for  our  finding  one  minimum  in  the 
delay  time  sequence  at  12—11. 

Literature  reports  on  the  effect  of  added  CO  on  laser  power  suggest  that  the  manner 

of  excitation  may  be  an  important  factor.  Pollack*, 19)  who  initialed  lasing  in  [CS  40  +He] 

2  2 

mixtures  by  flash  photolysis,  and  Arnold  and  Kimbell?0’  who  used  an  axial  spark  discharge, 
found  that  the  addition  of  CO  to  the  mixture  decreased  the  output  power.  On  the  other 
hand,  Suart,  Arnold  and  Kimbell(21)  observed  enhanced  lasing  when  they  added  cold  CO 

to  a  CS^/He  flow,  and  mixed  it  with  discharged  oxygen;  the  same  result  was  reported  by 

Foster (21a).  They  attributed  this  effect  to  an  increase  in  the  population  inversion 

density  through  selective  depopulation  by  V-V  exchange  of  the  lower  vibrational  states 

produced  via  chemical  excitation,  by  the  ground  state  carbon  monoxide.  Similarly, 

Jeffers  reported  a  factor  of  two  increase  in  power  in  a  transverse  flow  configuration 

when  he  added  CO  to  the  CS2  flow  and  subsequently  mixed  these  with  discharged 

[02  +  He],  A  positive  effect  due  to  carbon  monoxide  was  also  reported  by  Searle3  and 
(23) 

Djeu  in  a  CSg/Og  flame  laser.  Our  experimental  configuration  is  similar  to  that 
used  by  Arnold  and  Kimbell  and  we  found  that  the  addition  of  CO  to  the  mixture  resulted 
in  a  decrease  of  laser  power  both  in  grating  tuned  and  in  conventional  cavities. 

Figures  9  and  15  show  the  results  for  the  former  and  the  latter  cases,  respectively.  In 
these  experiments  the  added  CO  might  be  excited  by  the  pulsed  discharge,  but  the 
population  in  the  higher  levels  would  be  smaller  than  in  the  ground  state  by  several 
orders  of  magnitude.  This  is  suggested  by  experiment  and  by  theoretical  predictions*24, 2 5) 
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Our  observations  on  the  effect  of  added  CO  may  be  qualitatively  correlated  to  the  sub¬ 
stantial  production  of  CO(V^7)  via  reaction  (iii).  For  an  initial  distribution  such  as 
shown  in  Figure  11  the  chemically  unexcited  diluent  merely  raises  the  relative  popu¬ 
lation  in  the  low  v  states  to  produce  a  distribution  with  a  shallow  minimum,  one  that 
resembles  a  somewhat  relaxed  system.  This  accelerates  relaxation  to  an  equilibrium 
state  and  accounts  for  the  observed  shortening  of  laser  duration  for  each  transition, 
as  shown  in  Figure  8. 

/g\ 

Recent  measurements  of  relative  rates  for  (iii)  by  Hancock,  Morley  and  Smith  and 
by  Foster*21^  provide  no  information  for  vibrational  levels  v<  6;  where  their  data 
overlap  ours  the  agreement  is  good,  particularly  with  the  latter,  'n  view  of  the  assigned 
error  limits.  H.M.&S.  made  two  types  of  measurements.  In  one  they  estimated  the 
stationary  distribution  of  CO(V)  in  a  flowing  [CS2  +  O]  mixture.  The  recorded  light 
intensities  were  interpreted  in  terms  of  specific  chemical  excitation  rates,  using  the 
relative  spontaneous  emission  intensities  and  collisional  relaxation  rates.  They 
indicated  the  possible  presence  of  large  errors  in  the  rates  for  v  <  7 ,  because  the  latter 
were  given  by  the  difference  between  deactivating  rates  of  CO  from  v  to  lower  levels, 
and  from  higher  levels  to  v.  In  their  second  method  they  directly  monitored  the  popu¬ 
lation  of  excited  CO(V)  in  flash-photolyzed  [CS2  +  02  +  Ar]  mixtures,  using  the  CO  cw 
laser  lines  20  —  19  to  5  —  4  transitions  as  a  light  source.  Though  the  concentration  of 
CO  is  of  the  order  of  10~3  torr,  the  observed  distribution  at  50  jisec  after  the  flash  is 
determined  by  both  the  chemical  production  rate  and  V-V  relaxation.  Extrapolation  to 
zero  time  shows  <p .  and  <p_  w  0.  In  our  measurements  of  CO  emission  the  scanning  time 
of  the  monochromator,  which  required  2  hrs.  to  cover  the  CO  fundamental  region, 
limits  the  integration  time  of  the  signal,  and  the  signal-to-noise  ratio  is  not  as  good  as 
desired  for  the  weak  bands.  The  lower  transitions,  with  small  radiative  transition 
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probabilities,  do  not  overlap  each  other  as  much  as  do  the  higher  transitions.  In  any  case 
the  observed  spectrum  is  quite  different  from  the  one  calculated  using  the  distribution 
given  by  Hancock,  et.al. ,  as  seen  from  Fig.  10(c).  In  this  respect,  an  interesting  analysis 
was  presented  by  Dawson  and  Tarn'  '  who  calculated  the  stationary  vibrational  state 
distribution  of  CO,  including  chemical  formation  of  CO  in  each  level,  V-V  and  V-T 
energy  transfers,  spontaneous  radiation  loss,  and  removal  of  CO  by  pumping.  They  could 
explain  the  two  sets  of  experimental  CO  distributions  reported  by  Hancock  and  Smith 
by  using  the  initial  formation  rates  proposed  by  Hancock,  Morley  and  Smith*6*.  However, 
the  results  of  Arnold,  et.al.  which  requiied  a  distribution  with  a  large  rate  of  CO 
formation  was  difficult  to  simulate  without  assuming  CO  formation  in  the  low  levels. 

Since  the  distribution  with  a  high  formation  rate  of  CO  should  be  more  sensitive  to  the 
form  of  the  initial  distribution  of  CO,  this  conclusion  supports  the  initial  formation  rates 
of  CO  given  by  our  experiments.  Note  that  the  analysis  by  Dawson  and  Tam  was 
made  on  the  basis  of  calculated  V-V  rates  according  to  the  modified  SSH  theory,  which 
leads  to  rates  for  the  low  v  levels  that  are  smaller  by  a  factor  of  about  6,  than 
those  of  Hancock  and  Smith,  or  the  ones  estimated  from  the  Sharma-Brau  theory.  There¬ 
fore,  the  conclusion  derived  by  Dawson  and  Tam  requires  further  analysis. 

When  the  chemical  reaction  is  initiated  by  a  discharge  the  initial  conditions  may 
differ  significantly  from  those  iu  (O  +  CS^)  flow  systems.  We  assumed  that  the  plasma 
state  disappeared  within  about  1  fisec,  since  the  half-width  of  the  current  pulse  was 
about  0.3  fisec.  However,  ionic  or  electronically  excited  species  with  lifetimes  longer 
than  0.3  /isec  might  contribute  to  the  formation  of  CO  by  reactions,  such  as 
CS+  +  O  —  CO  +  S+;  CS+  +  O  —  CO+  +  S,  or  CS*  +0  —  CO  +  SO,  etc. 

Vibrational  states  of  CO  in  up  to  v  =  17  were  seen  in  the  present  system.  The 
relative  populations  at  v  =  13  ~  17  agree  well  with  those  based  on  the  flash-photolyzed 
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/M 

(CS  +  O  )  mixtures  as  reported  by  Hancock,  et.al. '  '  Using  the  new  value  for  the 
heat  of  formation  of  CS^14\  the  exothermicity  for  (iii)  is  84.6  kcal/mole;  this  permits 
excitation  to  v  =  15.  However,  as  suggested  by  Hancock,  et.al. ,  CS  is  also  vibrationally 
excited  since  CS  in  levels  v  =  0  ~  4  were  detected  in  our  U.  V.  absorption  experiments 
with  a  vibrational  temperature  of  about  12000°K*13\  If  the  vibrational  energy  of  the  CS 
can  be  converted  to  that  of  CO,  CS(v=4)  can  supply  14.6  kcal  extra  and  the  v  =  16  level 
for  CO  can  be  reached. 

SUMMARY 


(1)  In  the  pulse  initiated  [CS  +  O  +  He]  laser  system  there  is  an  optimum  discharge 

Ct  c* 

voltage  for  maximizing  laser  power.  The  optimum  voltage  is  higher  for  mixtures  lean 

in  O  and  rich  in  CS  .  The  maximum  power  of  the  laser  is  proportional  to  the  CS  content 
2  2  ^ 

In  O  -rich  mixtures. 

2 

(2)  When  the  discharge  voltage  is  higher  than  the  optimum,  CS  and  O  are  generated 
by  dissociation  of  the  reactants  in  the  discharge,  and  these  are  directly  involved  in  the 
lasing  reaction  [O  +  CS  —  CO^  +  S],  When  the  discharge  voltage  is  low,  a  chain 
reaction  is  initiated. 

(3)  The  delay  sequence  of  laser  pulses  for  each  vibrational  transition  was  recorded 
using  a  grating  cavity.  The  minimum  delay  occurs  for  the  12  —  11  transition,  with 
longer  delays  for  both  the  higher  and  lower  states. 

(4)  The  time-resolved  emission  spectra  of  CO[v  —  v-1]  fundamentals  were  observed 

(v) 

in  order  to  deduce  the  vibrational  distribution  of  chemically  generated  CO  .  The 
distribution  at  23  fisec  subsequent  to  the  discharge  pulse  shows  a  maximum  at  v  =  12, 
a  minimum  at  v  =  4,  and  a  population  that  increases  toward  lower  v. 

ft 

(5)  The  relative  population  for  each  vibrational  level  in  v  =  14  ~  5  was  calculated 
from  the  laser  delay  data;  the  results  were  found  to  be  consistent  with  those  obtained 
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from  the  CO  spontaneous  emission  measurements.  Also,  from  the  time-history  of 

the  vibrational  distribution  one  can  predict  qualitatively  the  power  and  the  duration  of 

the  lasing  pulse  for  each  vibrational  transition.  These  agree  with  observations. 

(V) 

(6)  The  shape  of  the  normalized  nascent  CO  population  distribution  [Fig.  11] 
suggests  the  presence  of  two  mechanisms  for  the  generation  of  vibrationally  excited 
CO,  but  it  is  not  known  which  reactions  other  than  [O  +  CS  —  CO  +  S]  contribute  to 
the  populations  in  these  lower  levels  for  pulse-discharged  [CS  +  O  +  He]  mixtures. 

a  a 
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TABLE  I 

Relative  Populations  of  CO  Vibrational  Levels 


Laser  Delay  Measurements 

]  1R  Emission  Measurements 

Assi 

a)  / 
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this  work 

Hancoek.  et.al. 
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0.13 
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0.18 
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0.30 
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0.32 

15 

0.49 

~0.2 

0.58 

14 

1.07 

0.93 

1. 06 

0.94 

1.02 

0.98 

0.72 

0.64 

0.90 

13 

0.93 

1.00 

0.90 

1.00 

0.87 

1.00 

1.00 

1.00 

1.00 

12 

0.90 

0.93 

0.85 

0.90 

0.80 

0.87 

1.03 

0.87 

0.85 

11 

0.88 

0.84 

0.83 

0.77 

0.75 

0.70 

0.90 

0.80 

0.65 

10 

0.89 

0.75 

0.81 

0.64 

0.69 

0.52 

0.80 

0.66 

0.55 

9 

0.88 

0.66 

0.78 

0.51 

0.55 

0.36 

0.71 

0.61 

0.41 

8 

0.87 

0.58 

0.73 

0.40 

0.18 

0.20 

0.62 

0.27 

0.32 

7 

0.90 

0.50 

0.72 

0.29 

-2.9 

0.036 

0.65 

~0.06 

0.17 

6 

0.89 

0.45 

0.60 

0.17 

0.61 

0.05 

5 

0.87 

0.39 

0.29 

0.05 

0.57 

4 

0.35 

0.015 

0.60 

3 

0.91 

2 

1.20 

1 

2.10 

(a)  Derived  from  O  +  CS  system 

2 

(t  Derived  from  flash-photolyzed  CS  +  O 

2  2 
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LEGENDS  FOR  FIGURES 


Figure  1. 
Figure  2. 


Figure  3. 
Figure  4. 
Figure  5. 
Figure  6. 

Figure  7. 


Figure  8. 


Figure  9. 


Figure  10. 


Figure  11. 


Schematic  for  tuned  laser  intensity  measurements. 

Signal  profiles  as  a  function  of  discharge  voltage  and  O  content. 

2 

Time  progresses  from  left  to  right. 

(a)  Abscissa  200  jusec/cm;  ordinate  0.2  V/cm. 

(b)  Abscissa  100  ^sec/cm;  ordinate  0.2  V/cm. 

Signal  profiles  as  a  function  of  He  content. 

Integrated  power  dependence  on  O  partial  pressure. 

Dependence  of  delay  times  on  operating  conditions. 

Integrated  laser  power  output,  as  a  function  of  vibrational  level, 
normalized  to  peak  value;  CS2  0.017  torr  +  C>2  0.68  torr  +  He  1.7  torr. 
Change  of  time  delay  for  initiation  of  P(19)  laser  pulse,  as  a  function 
of  vibrational  level  (abscissa:  20  jzsec/div;  ordinate  for  v  =  13  —  12 
through  v  =  6  —  5  transitions:  0.01,  0.05,  0.05,  0.02  V/div.) 

Time  history  of  P(9)  lasing  line,  for  no  added  CO  (a);  0.016  torr  CO 
added  (b);  0.032  torr  CO  added  (c). 

Dependence  of  integrated  power  for  P(9)  on  vibrational  transition  and 
the  amount  of  added  CO.  [Conditions  same  as  in  Figure  8.] 

Observed  (O)  vs  calculated  (T)  [via  eq.  5]  superposed  spectral  intensities 
of  CO  spontaneous  emissions:  0.045  torr  CS  +  0.76  torr  O  +3,0  torr  He 
15  kV  discharge;  (a)  23  ^sec,  (b)  148  jjsec  after  the  discharge  pulse, 

(c)  calculated  spectrum  on  the  basis  of  Hancock,  Morley  and  Smith 
distribution. 

Normalized  vibrational  population  at  selected  delay  times. 
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Figure  12. 


Figure  13. 


Figure  14. 


Figure  15. 


Total  population,  total  excited  state  population,  total  vibrational  quanta 

(v) 

and  average  quanta  per  CO  ,  as  a  function  of  time.  The  experimental 
data  are  from  Figure  11. 

Dependence  of  l/T  ^^(J)  onP  ratio*  ’I^ie  observed  (O)  were  normalized 
to  the  calculated  curves  at  P(ll). 

Relative  gain  vs  time  after  initiation,  measured  at  P(9)  for  various 
v—  (v-1).  Derived  from  experimental  data,  Fig.  11. 

Dependence  of  total  integrated  power  on  the  amount  of  CO  added, 
for  various  discharge  voltages. 


(For  total  power  measurement) 


FIGURE  1. 


Schematic  for  Tuned  Laser  Intensity  Measurements 
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FIGURE  2a 
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FIGURE  2b 


FIGURE  3 
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Partial  Pressure  of  02°r  CS2  (torr) 


Dependence  of  delay  times  on  operating  conditions 


Intensity  =  £  Pv  v-i^^e*a*'ve) 


Integrated  Laser  power -output  as  a  function  of 
vibrational  level,  normalized  to  peak  value  ,  CS 2 
O.OI7  torr  +  O2  0.68  torr  +  He  1.7  torr. 
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FIGURE  7.  Change  of  time  delay  for  initiation  of  P(10) 
laser  pulse,  as  a  function  of  vibrational 
level  (abscissa:  204s/div;  ordinate  for 
v  =  13— 12  through  v  -  6— 5  transitions: 
0.01,  0.05  ,  0.05  ,  0.02  v'div). 
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FIGURE  8.  Time  history  of  P(9)  lasing  line,  for  no  added  CO(a);  0.016  torr  CO 
added  (b)-,  0.032  torr  CO  added  (c). 
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FIGURE  9.  Dependence  of  integrated  power  for  P(9)  on  vibrational  transition 
and  the  amount  of  added  CO.  [Conditions  same  as  in  Figure  12], 


Normalized  Population 


Fig.  11  Normalized  vibrational 
population  ot  selected 
deloy  times.  [Experimental 
conditions  ore  the  same 
os  in  Fig.  15]  . 
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Fig.  13  Dependence  of  l/T,2—  M(0)  on  4>  ratio.  The  observed  points  (o)  were 
normalized  to  the  calculated  curves  ot  P  (II). 
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C.  Studies  of  the  0  -  C2H2  Reaction  (R.  A.  McFarlane) 

The  reaction  of  atomic  oxygen  with  acetylene  can  lead  to  the  production 
of  vibrationaliy  excited  carbon  monoxide  via  either  or  both  of  the  following 
reactions. 

0  +  C2H2  -*•  CH2  +  CO*  +  51  kcal/mole  (1) 

0  +  CH2  CO*  +  2H  +  71  kcal/mole  (2) 

It  has  been  the  put  pose  of  our  studies  during  the  past  period  to  determine 
whether  as  a  result  of  the  reactions  a  population  inversion  exists  between 
vibrationaliy  excited  levels  of  the  CO  ground  electronic  state.  A  chemical 
laser  based  on  this  reaction  and  operating  in  the  5  u  spectral  region  would 
have  significant  advantages  over  present  systems  in  regard  to  fuel  costs 
and  greatly  simplified  handling  of  both  fuel  and  reaction  products. 

Our  report  for  the  previous  period  described  the  atomic  oxygen  source 
using  a  1  kw,  c.w.  magnetron  and  described  our  measurements  on  atom 
concentration  available  for  our  studies  of  the  above  reactions.  The 
accomplishments  of  the  present  period  are  discussed  in  detail  in  the 
accompanying  papers  which  have  been  published  or  submitted  for  publica¬ 
tion. 

Measurements  of  Vibrational  Population  Distribution 

Following  the  injection  of  acetylene  into  the  fast  flowing  stream 
of  atomic  oxygen  (typically  10  ^  mole/sec)  observation  was  made  of 
the  first  overtone  emission  from  the  chemically  formed  carbon  monoxide 
in  the  spectral  range  2.3  -  2.8  microns.  Scans  were  made  using  a  one-half 
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meter  monochrometer  equipped  with  a  room  temperature  ?bS  cell  and  the 
system  sensitivity  was  measured  using  a  600°K  blackbody  source. 

A  computer  program  was  set  up  to  produce  a  synthetic  spectrum  using 
as  Input  parameters  the  vibrational  populations  and  the  rotational 
temperature.  The  program  Included  the  spectroscopic  constants  of  the  CO 
molecule,  the  system  sensitivity  and  the  spectroscopic  slit  function.  The 
input  parameters  were  adjusted  to  give  a  synthetic  spectrum  which 
matched  the  observed  scan  and  by  this  means  it  was  possible  to  determine 
the  relativa  populations  of  levels  v  ■  2  to  v  ■  14  under  a  wide  variety 
of  reaction  conditions .  With  the  large  amount  of  Helium  diluent  used  the 
rotational  temperature  was  typically  though  not  always  very  close  to 
300 °K. 

At  low  fuel  levels  corresponding  to  low  concentrations  of  chemically 
formed  CO,  a  near  Boltzman  distribution  of  population  was  observed  correspond¬ 
ing  to  a  vibrational  temperature  of  approximately  8000°K.  As  the  fuel  level  was 
incraased  the  effect  of  the  Traanor  pumping  mechanism  became  apparent 
wheraby  low  vibrational  lavels  ware  found  to  have  a  local  vibrational 

temperature  of  3500°K  and  the  higher  vibrational  lavels,  v  ■  5  to  v  ■  14 
a  local  vibrational  temperature  as  high  as  14000°K.  This  lattar  value 
is  typical  of  that  found  in  alectrically  excited  CO  lasers  and  is  adequate 
to  provide  a  partial  inversion  on  P  branch  transitions  of  J  >  20. 

Studies  were  made  of  the  affect  on  the  vibrational  population  of 
adding  cold  CO  premixed  with  the  C2H2  fuel.  The  process 


CO(v-O)  +  CO (v)  -*-C0(l)  +  CO(v-l) 


(3) 
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can  be  expected  to  dramatically  modify  the  distribution  where  the  probability 
for  this  one  quantum  exchange  is  large.  At  a  temperature  of  300°K  this 
probability  is  largest  for  small  v  and  peaks  near  v  ■  3  -  4.  The  addition 
of  0.1  torr  cold  CO  resulted  in  a  very  high  but  positive  local  vibrational 
temperature  for  levels  v  ■  5  and  v  ■  6  while  the  addition  of  0.2  torr  cold 
CO  resulted  in  a  total  vibrational  inversion  for  levels  v  ■  5  to  v  ■  7. 

There  are  therefore  two  different  regimes  in  which  the  oxygen-acetylene 
reaction  could  be  employed  for  laser  purposes. 

1.  High  fuel  condition  -  taking  advantage  of  Treanor  pumping 
to  populate  high  v  levels. 

2.  Preferential  V-V  relaxation  of  lower  vibrational  levels  in  a 
time  short  compared  with  V-T  processes. 
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In  Ira  red  chemiluminescence  Irom  Hie  reaction  of  atomic  oxygen  with  acetylene  was  studied  lor  CO  laser  applica¬ 
tion.  A  totally  inverted  vibrational  population  between  a  =  5  and  v  =  7  was  obtained  by  preferentially  relaxing  the 
lower  vibrational  levels  with  added  cold  CO.  Treanor  pumping  at  high  fuel  flows  was  also  observed. 


In  this  letter  we  report  the  first  observation  of  a  to¬ 
tally  inverted  vibrational  population  distribution  in 
carbon  monoxide  formed  in  the  chemical  reaction  of 
oxygen  with  acetylene.  It  is  apparent  from  the  results 
that  this  reaction  is  potentially  valuable  for  use  in  CO 
laser  applications. 

The  primary  reactions  are  [  1 1 

O  +  C2H2  -*  CO  +  CH2  +  5 1  kcal/mole , 

O  +  CH2  -  CO  +  2H  +  7 1  kcal/tnole  (1) 

for  low  fuel  flows.  These  reactioi,:  are  sufficiently  en¬ 
ergetic  to  permit  the  excitation  of  vibrational  levels  in 
CO  up  to  p=  14.  Previous  studies  [2)  have  found  that 
as  a  consequence  of  the  anharmonic  nature  of  the  mo¬ 
lecular  potential  the  relative  populations  in  higher  vi¬ 
brational  levels  can  be  increased  through  intramolecular 
V-V  energy  transfer,  producing  a  highly  non-Boltzmann 
distribution.  In  the  present  work,  advantage  is  taken 
of  the  V-V  energy-transfer  process 

CO(v)  +  CO(v=0)  -*■  CO(e- 1 )  +  CO(v=  1)  +  Ap  . 
Calculations  [2,3|  of  the  probability  of  this  process 

f  Research  supported  by  the  Advanced  Research  Projects 
Agency  and  monitored  by  ONR  under  Contract  N00014-67- 
A-0077-0006,  and  also  through  the  Materials  Science  Center, 
Cornell  University. 


show  that  it  peaks  between  t>=4  and  v=6  for  tempera¬ 
tures  between  300  and  700°K,  decreasing  only  slight¬ 
ly  as  resonance  is  approached  with  decreasing  v  and 
falling  off  much  more  rapidly  at  higher  values  of  t>.  It 
has  been  possible  therefore  with  the  addition  of  cold 
CO  to  the  reaction  to  significantly  depress  the  popula¬ 
tions  of  levels  v  <  5  and  provide  a  total  vibrational  in¬ 
version  between  v=S  and  v=7. 

Using  a  fast  flow  system,  1R  chemiluminescence  of 
CO*  was  measured  in  the  fundamental  and  the  first 
overtone  emission  regions.  Atomic  oxygen  at  flow  rates 
in  excess  of  10~4  mole/sec  was  produced  in  a  micro- 
wave  discharge,  the  concentration  being  established  by 
titration  with  N02.  Acetylene  was  injected  radially  in¬ 
to  a  teflon  reaction  tube  3/8  inch  i.d.  and  18  inches 
long.  The  system  was  pumped  with  a  100  c.f.m.  blow¬ 
er  backed  by  a  model  1397B  Welch  mechanical  pump. 
A  linear  flow  speed  of  greater  than  80  m/sec  in  the 
mixing  zone  was  measured  under  typical  conditions. 

Emission  from  the  reaction  products  was  monitored 
through  NaCl  windows  in  the  side  wall  of  the  reaction 
tube.  A  silicon  lens  focused  the  light  from  a  small  re¬ 
gion  in  the  reaction  zone,  through  a  chopper,  onto  the 
entrance  slit  of  a  1/2  meter  Jarrell-Ash  Ebert  mono¬ 
chromator,  equipped  with  a  105  1/mm  grating  blazed 
at  8.6  microns  in  first  order.  A  room-temperature  PbS 
photoconductive  detector  was  used  for  the  overtone 
measurements  and  a  liquid-nitrogen  cooled  Au  -Ge  de¬ 
tector  for  observations  in  the  fundamental  region.  The 
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detecior  ouiput  signal  was  amplified  by  a  Keithley 
Model  1 03  low-noise  amplifier  and  a  Princeton  Applied 
Research  Model  JB-5  lock-in  amplifier.  Specira  were 
recorded  on  a  linear  strip  chart  recorder.  The  entire 
optical  path  was  Hushed  with  dry  nitrogen  to  minimize 
absorption  by  atmospheric  carbon  dioxide  and  water 
vapour. 

A  broad-band  emission  of  substantial  intensity 
throughout  the  region  of  the  CO  fundamental  made  it 
difficult  to  determine  the  vibrational  populations  from 
the  fundamental  emission.  The  first  overtone  emission 
spectrum  between  2.3  and  2.8  microns  was  distinct 
however  and  permitted  the  determination  of  the  rela¬ 
tive  population  distribution  by  a  method  described  by 
Karl  et  al.  [4|.  A  computer  is  used  to  generate  a  simu¬ 
lated  spectrum  to  compare  with  the  experimental 
scans.  It  accounts  for  the  instrumental  sensitivity  vari¬ 
ation  as  a  function  of  wavelength,  the  spectral  slit 
width,  the  Einstein  A  coefficient  of  each  upper  level 
and  the  molecular  spectroscopic  constants.  The  rota¬ 
tional  temperature  and  the  relative  population  of  each 
vibrational  level  are  adjusted  to  obtain  a  match  be¬ 
tween  the  synthetic  spectrum  and  the  experimental 
scan. 

Fig.  I  shows  the  vibrational  population  distribution 
of  CO*  at  a  time  approximately  100  microseconds 
after  injection  of  the  acetylene.  The  distribution  did 
not  alter  appreciably  for  observations  made  several 
centimeters  downstream.  An  addition  of  cold  CO  [5] 
to  the  fuel  dramatically  modifies  the  vibrational  dis¬ 
tribution  leading  to  a  total  vibrational  inversion  with 
0.2  torr  CO  added  to  the  acetylene.  The  data  have  been 
plotted  taking  account  of  relative  signal  levels  under 
different  conditions  of  added  cold  CO.  The  data  is 
clearly  consistent  with  the  peak  in  the  V-V  exchange 
probability  predicted  to  be  near  i>=5  although  no  at¬ 
tempt  has  yet  been  made  to  set  up  a  comprehensive 
model  which  includes  reactant  mixing,  the  detailed 
chemistry,  V  V  and  V-T  relaxation  and  radiative 
decay.  Our  observations  at  fuel  Hows  six  times  that  of 
fig.  I  but  without  the  addition  of  cold  CO  show  clear¬ 
ly  the  effects  of  V-V  pumping  predicted  by  Treanor 
et  al.  [2|.  The  vibrational  temperature  of  approximate¬ 
ly  8000°K  which  characterizes  the  Boltzmann  distribu¬ 
tion  of  population  at  very  low  fuel  levels  becomes  mod¬ 
ified  at  high  fuel  levels.  A  local  vibrational  tempera¬ 
ture  of  approximately  350()°K  was  determined  for 
levels  v=2,  to  v=4,  while  for  levels  v=6  to  v=  14  popu¬ 
lations  were  characteristic  of  a  !4000°K  vibrational 
temperature.  This  latter  value  is  typical  of  electrically 
560 
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l;ig.  I.  Relative  vibrational  population  distribution  of  CO*. 
Pressures:  C2  H2  0.05  lorr,  02  0.68  torr.  He  1.67 

torr.  Temperature  rrot  =  350°K.  |C2*H2 1  =  0.25  X  It)-4 

mole/sec,  |0)  =  1.7  X  10 -4  mole/ sec.  o  Without  addi¬ 
tional  told  CO.  ii  0.05  lorr  additional  cold  CO.  A  0.1  ton 
additional  cold  CO,  e  0.2  torr  additional  told  CO. 

excited  CO  lasers  which  operate  under  conditions  of 
partial  inversion. 

It  is  apparent  that  both  modes  of  operation  are  po¬ 
tentially  useful  for  CO  laser  application.  Since  the  re¬ 
action  rate  of  oxygen  with  acetylene  is  greater  than 
that  for  the  reaction  of  oxygen  with  carbon  monoxide 
the  above  results  also  suggest  a  free  burning  flame 
laser  in  which  a  mixture  of  acetylene  and  carbon  mon¬ 
oxide  is  burned  in  an  oxygen  atmosphere. 


References 

1 1 1  D.M.  Creek,  C.M  Melliar-Smilh  and  N.  Jonathan,  J.  Chem. 
Soc.  A  (1970)646. 

|2|  C.E.  Treanor.  J.W.  Rich  and  R.G.  Reh  n,  J.  Chem.  Ph)s. 
48(1968)  1728. 

13|  W.Q.  Jeffers  and  J.D.  Kelley,  J.  Chem.  Phys.  55  ( 1971 ) 
4433. 

|4 1  G.  Karl,  P.  Kraus  and  J.C.  Polanyi,  J.  Chem.  Phys.  46  ( 1 967) 
224. 

151  G.  Hancock  and  l.W.M.  Smith,  Appl.  Opt.  10  (1971)  1827. 


-105- 


Measurement  of  V-V  Energy  Transfer  Probabilities  In  CO-CO  Collisions 

The  population  data  obtained  upon  the  addition  of  cold  CO  to  the 
system  were  used  to  determine  the  V-V  energy  exchange  probabilities  of 
process  (3)  for  excited  vibrational  levels  v  ■  3  to  v  ■  9.  Under  the 
conditions  of  the  experiment  the  volumetric  pumping  rate  greatly  exceeded 
all  other  loss  processes  for  a  given  vibrational  level  population.  This 
permits  the  series  of  coupled  steady  state  equations  describing  the 
population  densities  to  be  uncoupled  and  the  collisional  deactivation 
rate  due  to  added  cold  CO  to  be  determined  separately  for  each  vibrational 
level.  The  detailed  analysis  is  discussed  in  an  appended  paper. 

Reoionably  good  agreement  is  obtained  with  data  reported  by  Hancock, 
and  Smith  in  studies  of  the  CS^-0^  reaction  and  under  very  different 
flow  conditions.  The  "Fast  Flow  Approximation"  used  in  the  present  study 
substantially  simplifies  the  analysis  of  the  data  and  avoids  uncontrolled 
wall  effects  and  others  which  would  complicate  the  interpretation. 

Calculations  were  made,  following  the  work  of  Caledonia  and  Ceater, 
of  the  V-V  exchange  probabilities  for  single  quantum  exchange  in  collision. 

Both  long  and  short  range  interactions  are  included  in  the  calculation  and 
data  is  now  available  over  a  sufficiently  wide  range  of  vibrational  levels 
to  indicate  that  both  contributions  are  required  to  correctly  describe  the 
collisional  process.  Good  agreement  was  obtained  with  our  experimental 
measurements  and  those  of  Hancock  and  Smith. 

A  Time-Dependent  Solution  of  the  Kinetic  Equations  for  the  Vibrational  Populations 

For  the  system  under  study  it  is  neceisary  to  provide  some  redistribution 
of  vibrational  energy  in  the  carbon  monoxide  before  conditions  of  total  or 
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partial  inversion  maintain.  This  process  is  a  dynamic  one  and  for  laser 
applications  it  is  necessary  to  know  at  what  time  optimum  population 
conditions  occur.  A  computer  solution  was  obtained  for  the  coupled 
equations  describing  the  population  densities  for  each  vibrational  level 
as  a  function  of  time.  This  is  discussed  in  detail  in  an  appended  paper. 

The  mixing  and  chemistry  are  assumed  to  be  very  rapid  and  the  model 
includes  all  single  quantum  V-V  energy  exchanges,  V-T  energy  exchange 
and  spontaneous  radiative  decay. 

The  results  of  calculations  demonstrate  the  preferential  deactivation 
of  low  v  levels  discussed  above  as  well  as  the  Treanor  pumping  of  high  v 
levels.  In  general  it  is  found  that  the  populations  relax  through  a  series 
of  non  Boltzmann  distributions  and  at  long  times  the  populations  do  indeed 
approach  a  Boltzmann  distribution. 

Using  as  initial  conditions  those  for  the  experiments  on  the  effect 
of  adding  cold  CO  to  the  fuel,  it  was  possible  to  reproduce  the  observations 
that  0.1  torr  added  cold  CO  gives  an  increased  but  positive  local  vibrational 
temperature  while  0.2  torr  leads  to  a  total  vibrational  inversion. 

The  computation  was  used  to  study  the  time  development  of  CO  population 
using  results  for  the  initial  population  distribution  from  the  CS2-02  reaction 
determined  by  Hancock  and  Smith  and  by  Tsuchiya,  Nielsen  and  Bauer.  Only 
preliminary  calculations  have  been  made  to  date  but  it  is  clear  from 
comparing  them  with  the  experimental  time  development  data  reported  by 
the  Cornell  group,  that  substantial  populations  in  v  ■  0  and  v  ■  1  are 
needed  to  reasonably  simulate  the  time  development  found.  The  modeling  for 
this  system  described  in  the  appended  paper  shows  an  overall  time  development 
which  is  somewhat  too  slow  compared  to  experiment  and  appears  to  underestimate 
the  magnitude  of  Treanor  pumping  which  is  operative. 
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To  accurately  describe  any  of  the  systems  where  the  computation  might  be 
applied  it  is  clearly  desirable  to  include  chemical  reactions  as  time 
dependent  contributors  to  the  populations.  Using  known  or  estimated 
reaction  rates  this  represents  a  trivial  extension  of  the  present  program. 
Recent  calculations  of  multiple  quantum  V-V  processes  by  Dillon  and 
Stephenson  will  be  Included  in  future  computations . 

Continuing  Activity 

In  addition  to  the  elaboration  of  the  computation  discussed  above 
work  is  in  progress  to  set  up  a  laser  l  ased  on  the  O-CjHj  reaction  and  to 
carry  out  parametric  studies  of  its  performance.  To  this  end  we  have: 

a)  Begun  the  installation  of  a  2000  cfm.  blower  system  to  provide 
an  order  of  magnitude  increase  in  pumping  capacity  over  that 
now  being  used.  This  will  be  used  to  permit  operation  of  a 
transverse  flow  laser  system. 

b)  Constructed  a  c.w. ,  grating  tuned  CO  laser  for  diagnostic 
studies  of  both  the  present  linear  flow  system  and  the 
planned  transverse  flow  system. 

c)  Constructed  an  electronic  system  for  differential  gain  and 
loss  measurements  on  any  chemical  laser  system  of  interest. 

TVo  pyroelectric  detectors  are  employed  and  normalization 
circuitry  is  provided  to  eliminate  errors  resulting  from 
probe  laser  power  fluctuations. 
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Meajurement  of  V-V  Energy  Transfer  Probabilities 
in  CO-CO  Collisions  Following  the  Reaction 
of  Oxygen  with  Acetylene. 

Y.  S.  Liu,  R.  A.  McFarlane  and  G.  J.  Wolga 


Abstract 

Infrared  chemiluminescence  has  been  studied  from  vibrationally 
excited  CO  formed  in  the  reaction  of  oxygen  with  acetylene.  By 
monitoring  the  first  overtone  emission  at  various  levels  of  added 
cold  CO  in  a  fast-flow  system,  we  were  able  to  determine  the  V-V 
energy  transfer  probabilities  ^  for  3  £  v  £  9.  The 

analysis  is  greatly  simplified  by  using  a  fast-flow  approximation 
in  which  the  volumetric  pumping  rate  greatly  exceeds  the  removal 
rate  for  other  loss  processes.  The  results  are  compared  with 
previously  reported  data  and  with  the  probabilities  determined 
theoretically  where  both  long  and  short  range  interaction  potentials 
are  Included. 
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INTR0DUCTI0N 


In  recent  years  studies  of  inter-  and  intra-  molecular  energy 
transfer  in  diatomic  molecules  have  been  motivated  by  interest  in 
chemical  lasers.  Various  experimental  methods  have  been  employed  to 
measure  the  rates  of  vibrational  energy  transfer  which  occur  upon 
collision.  Among  them  I.R.  chemiluminescence  has  been  exploited 
most  extensively.  The  method  has  been  successfully  applied  by  Polanyi 
and  his  co-workers ^  to  determine  the  relative  rate  constants  for  the 
production  of  specific  energy  levels  of  the  product  molecule  in  hydrogen 

halide  reactions  and  many  others.  More  recently,  the  quenching  of  I.R. 

(2) 

chemiluminescence  has  been  used  by  Hancock  and  Smith'  '  to  measure  the 
rates  of  de-excitation  of  vibrationally  excited  CO  by  ground  state  CO 
and  other  molecules.  They  observed  the  I.R.  chemiluminescence  from 
vibrationally  excited  CO  formed  in  the  reaction  of  oxygen  and  carbon 
disulfide  and  studied  the  quenching  of  the  CO  first  overtone  spectrum 
upon  the  addition  of  various  gases. 

There  are  two  times  that  have  to  be  considered  in  I.R.  chemiluminescence 


experiments:  the  spontaneous  radiative  lifetime  x  ^  and  the  col- 

spont . 

lisional  de-activation  time  x  , .  of  the  molecule.  The  former  is 


—  bxmw  *  11,  ,  w*  wilt-  lUUACUUlCI  iUC  lULUlbl  0.0 

collision 

characterized  by  the  spontaneous  transition  probabilities  Av,  ;  the 
latter  depends  on  the  collisional  de-activation  probability  and  on  the 
density.  For  the  vibrationally  excited  CO  molecule,  the  spontaneous 

_3 

radiative  lifetime  is  on  the  order  of  10  seconds.  A  pressure  of  a  few 
Torr  is  required  for  the  collisional  de-activation  time  to  be  of  the 


same  order. 


Chemiluminescence  experiments  are  frequently  carried  out  using  a 
low  flow  condition  eet  by  a  relatively  small  capacity  pump.  The  low 
pressure  required  is  established  by  moving  only  a  small  amount  of 
material.  In  the  low  flow  system,  however,  undesirable  and  unavoidable 
experimental  problems  such  as  the  de-excitation  due  to  surface  effects, 
the  cascading  process  from  the  vibrationally  excited  states  and  the 
spatial  inhomogeneity  of  the  species  have  to  be  taken  into  account. 
These  can  greatly  affect  the  accuracy  of  rate  determinations  based  on 
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I.R.  measurements  and  can  lengthen  the  analysis  of  the  experimental 
results. 

In  view  of  these  difficulties  encountered  in  the  conventional 
I.R.  chemiluminescence  experiments  using  the  low  flow  condition,  it 
appears  more  appropriate  to  operate  I.R.  chemiluminescence  experiments 
at  a  high  flow  rate.  The  pressure  can  be  reduced  by  increasing  the 
volumetric  pumping  rate.  The  "Fast  Flow"  condition  is  satisfied  by 
carefully  choosing  the  experimental  parameters  such  that  the  volumetric 
pumping  rate  is  made  much  faster  than  the  spontaneous  radiative  rate 
and  collislonal  deactivation  rates.  Under  these  conditions,  using 
the  appropriate  approximations,  the  master  rate  equations  can  be 
de-coupled  and  solved  in  a  very  simply  way.  The  approach  is  hence  called 
the  "Fast  Flow  Approximation". 

The  vibrational  energy  exchange  probabilities  measured  in  the 

(2) 

present  work  are  compared  with  those  obtained  by  Hancock  and  Smith' 
in  an  I.R.  quenching  chemiluminescence  experiment  and  with  theoretical 
values  calculated  by  including  both  long-range  dipole-dipole  interactions 
and  the  short  range  repulsive  Interactions.  Reasonably  good  agreement 
is  obtained  between  the  present  values  and  those  reported  by  Hancock 
and  Smith.  This  is  particularly  satisfying  in  view  of  the  very  great 
experimental  differences  in  the  two  determinations.  Both  sets  of 
data  are  also  in  good  agreement  with  the  calculated  V-V  energy  exchange 
probabilities . 

THE  EXPERIMENT 

Vibrationally  excited  CO  is  produced  in  the  primary  reactions:'  * 

0  +  C2H2  CO*  +  CH2  +  51  Real /mole 

(1) 

0  +  CH2  - ^  CO*  +  2H  +  71  Kcal/mole  . 
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The  exothermicity  of  the  reactions  is  sufficient  to  produce  vibrationally 
excited  CO*  up  to  v  *  14.  In  the  present  work  chemiluminescence  of  CO* 
was  measured  in  the  first  overtone  emission  bands  between  2.3  and  2.8 
microns  as  shown  in  Figure  1(a).  The  spectra  were  monitored  as  different 
amounts  of  cold  CO  were  added  to  the  acetylene  and  injected  into  the 
mixing  zone.  The  addition  of  cold  CO  affects  not  only  the  overall 
intensity  of  the  spectrum  but  changes  the  intensity  distribution  as 
well.  This  latter  effect  is  caused  by  the  preferential  deactivation 
of  the  lower  vibrational  levels  of  the  vibrationally  excited  CO 
molecules.  Figures  1(b)  and  1(c)  show  two  spectra  recorded  at  different 
levels  of  added  cold  CO. 

The  experimental  apparatus  is  shown  in  Figure  2.  Atomic  oxygen  at 
a  flow  rate  of  10-^  mole/sec  was  produced  in  a  microwave  discharge.  The 
oxygen  concentration  was  determined  by  titration  with  NO^*  Acetylene 
was  injected  radially  through  small  orifices  into  a  teflon  reaction  tube 
3/8  inch  i.d.  and  18  inches  long.  Precautions  were  taken  in  choosing 
the  dimensions  to  assure  that  the  low  pressure  and  fast  flow  conditions 
were  satisfield.  The  system  was  pumped  with  a  100  cfm  blower  backed  by 
a  moder  1397B  Welch  pump.  A  linear  flow  speed  of  83  m/sec  in  the  mixing 
zone  was  measured  under  normal  operating  condition. 

The  infrared  emission  was  observed  in  a  direction  transverse  to  the 
flow.  The  observations  were  made  through  a  side  window  at  various 
positions  along  the  flow  downstream  from  the  mixing  point.  The  infrared 
emission  was  focussed  with  a  silicon  lens  through  a  chopper  onto  the 
entrance  slit  of  a  Jarrell-Ash  1/2  meter  Ebert  monochromator.  A  room 
temperature  PbS  detector  was  used.  The  signal  was  amplified  by  a 
Keithley  Model  103  low  noise  amplifier  and  processed  through  a  Princeton 
Applied  Research  Model  JB-5  lock-in-amplifier  phase  locked  at  a  reference 
frequency  set  by  the  light  chopper.  The  spectra  were  recorded  on  a 
Hewlett-Packard  dual-pen  strip  chart  recorder.  The  spectrometer  employed 
a  grating  with  105  lines/mm  blazed  at  8.6  microns  in  first  order.  The 
optical  path  internal  and  external  to  the  spectrometer  was  flushed  with 
dry  nitrogen  to  avoid  the  atmospheric  water  vapor  absorption  which  overlaps 
strongly  with  the  CO  fundamental  and  somewhat  less  with  the  overtone. 
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The  vibrational  population  distributions  were  determined  by  a 
computer  simulation  method  introduced  by  Karl,  Kruss  and  Polanyi^. 

A  spectrum  can  be  synthesized  to  match  the  experimental  scan  with 
vibrational  populations  and  the  rotational  temperature  as  the  input 
parameters  and  using  the  molecular  and  spectroscopic  constants,  the 
frequencies  of  transition  and  the  instrumental  characteristics.  A 
preliminary  discussion  of  the  determination  of  the  vibrational 
populations  for  the  present  system  has  been  reported  In  the 

present  work  we  will  describe  the  determination  of  the  V-V  energy 
exchange  probability  for  the  process 

C0(0)  +  CO(v)  -*■  C0(1)  +  CO(v  -  1)  +  AE  (2) 

using  the  vibrational  population  data  of  reference  (5). 

A  STEADY  STATE  SOLUTION  FOR  VIBRATIONAL  POPULATIONS  USING  THE  FAST-FLOW 
APPROXIMATION 


The  steady  state  master  equation  for  the  population  densities  N 

\ 

of  the  vibrational  level  v  can  be  expressed  as 


0  - 


dt 


N  -  R  + 
v  v 


max 


n*l 


(A  +  K  [X]  }N 

v  t*  n | v  v  +  n *  v  4*  n 


(3) 


U, 


v,v  -  n 


+  K 


v,v 


n**l 


„[X)  ]  +  K*  +  Kp>Nv 


The  positive  term  represents  the  total  production  rate  of  the 
population  density  into  vibrational  level  v.  The  negative  term 
describes  he  rate  of  removal  of  molecules  from  the  vibrational  level  v. 
The  notation  is  as  follows: 

is  the  chemical  production  rate  that  produces  population  in 
the  vibrational  level  v. 

Aviy  is  the  spontaneous  radiative  rate  of  the  transition  from 

v*  to  v. 


i 
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Kv>v  is  the  colllslonal  deactivr ,iOn  rate  at  which  a  molecule 

is  moved  from  the  state  v*  to  the  state  v. 

[X]  is  the  concentration  of  molecule  X  that  causes  deactivation. 

Only  CO  is  included  in  the  following  analysis, 
is  the  total  rate  of  deactivation  by  surface  recombination 
and  diffusion. 

Kp  is  the  volumetric  pumping  rate. 

The  master  equation  can  be  simplified  with  the  following 
approximations. 

(1)  The  Fast  Flow  Approximation 


Kv  +  n,v[X] '  Av  +  n,v*  Kw  <  <  Kp 


(3)  V  K  .  -  K  .  , 

L.  v  +  n,v  v  +  l,v 

n"l 

(4)  Rv  •  constant 

The  first  approximation  is  valid  if  the  volumetric  pumping  rate  is  much 
faster  than  the  spontaneous  radiative  rate,  the  colllslonal  deactivation 
rate  and  diffusion  rate.  The  second  and  the  third  assumptions  imply  that 
the  spontaneous  radiation  is  confined  in  the  fundamental  and  the  first 
overtone  emissions  and  colllslonal  energy  transfer  is  a  single  quantum 
exchange  process.  Further,  we  assume  that  the  reaction  rate  is  not 
affected  by  the  addition  of  another  gas. 

With  these  approximations,  the  steady  state  master  equation 
without  and  with  cold  CO  can  be  written  as. 


5 


-115- 


0  "  97  Nv°  “  Rv  +  *Av  +  l,v  +  Kv  +  l,v  M}Nv  +  1  +  Av  +  2  vNv  +  2 


"  {Av,v  -  1  +  Av,v  -  2  +  Kv,v  -  1  ^  +  Kp  +  K„>  N° 


¥t  Nv  “  Rv  +  {Av  •»-  l,v  +  Kv  +  l.v  ^  }Nv  +  1  +  Av  +  2  vNv  +  2 


r,v  -  1  +  Av,v  -  2  +  Kv,v  -  lM  +  *p  +  *w> 


where 


[q]  is  the  concentration  of  CO  formed  in  the  reaction. 

[Q]  is  the  concentration  of  CO  added.  [Q]  »  [q]  . 

N°  is  the  concentration  of  vibrational  population  in  level  v 

without  additional  cold  CO. 

Nv  i8  the  concentration  of  vibrational  population  in  level  v  with 
the  additional  cold  CO. 

If  we  subtract  and  rearrange  the  two  equations  and  keep  the  significant 
terms,  we  obtain 

jr  *  x*  * 1  ^ <7> 

If  the  left-hand  side  is  plotted  against  the  magnitude  of  the  bracket  on  the 

right-hand  side  using  the  measured  values  of  N  and  N°  and  the  values 

v  v 

[Q]  of  the  additional  cold  CO,  Ky  +  can  be  determined.  The 

concentration  of  CO  formed  in  the  reaction  [q  ]  ,  is  determined  using  the 
known  reaction  rate  constant  and  the  known  concentrations  of  oxygen  and 
acetylene.  It  is  estimated  as  [q]  -  5.4  x  10“10  mole/cm3.  Together  with 
the  values  of  Kp  and  A v  +  the  values  of  the  brackets  in  Equation  (7) 

can  be  determined  for  each  vibrational  level  v.  The  results  are  shown  in 
Figure  3  for  the  vibrational  levels  v  -  2,4,  and  8.  The  deactivation 
rat®  constant  Kv  +  ^  can  be  determined  from  the  slope. 
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The  V-V  energy  transfer  probabilities  of  the  process  (I)  can  be 
determined  from  the  measured  values  of  the  de-activation  rate  constants. 

pO.l  +  Kv.v  -  1 
v-v  -  1  zco  -  CO  • 

14  -1  -3 

The  CO-CO  collision  frequency  is  calculated  to  be  1.86  x  1C  mole  cm 
sec  ^  at  a  temperature  of  300°K. 

Our  determinations  of  K  -  and  P^’^  ,  are  summarized  in 

V.V  -  1  v,v  -  1  (2) 

Table  1.  Also  shown  are  the  data  reported  by  Hancock  and  Smith'  1  from 
their  measurements  on  the  quenching  of  CO  chemiluminescence  resulting  from 
the  reaction  of  oxygen  and  carbon  disulfide. 

CALCULATION  OF  THE  VIBRATIONAL  ENERGY  TRANSFER  PROBABILITIES 

The  V-V  energy  transfer  probabilities  for  the  process 

C0(0)  +  CO(v)  -►  C0(1)  +  CO(v  -  1)  +  AE 

must  be  calculated  including  both  long  range  and  short  range  Interaction 
potentials.  In  general  the  V-V  energy  transfer  in  CO-CO  collisions  can  be 
resonant  or  highly  non-resonant  depending  on  the  energy  levels  involved.  No 
single  theory  can  adequately  explain  the  energy  transfer  mechanism.  The 
short  range  interaction  from  the  molecular  repulsive  potential  has  been 
treated  by  Rapp  et  al. ^  and  the  long  range  interaction  due  to  molecular 
multipole  moments  has  been  treated  by  Sharma  and  Brau^.  Detailed 

calculations  Including  both  contributions  to  the  transition  probability  have 

(8)  (9) 

been  reported  by  Jeffers  and  Kelly  and  by  Caledonia  and  Center v 

Figure  4  shows  the  measured  V-V  exchange  probabilities  compared  to  our 

calculated  values  based  on  the  formulation  of  Caledonia  and  Center.  The 

experimental  data  supports  the  notion  that  both  long  and  short  range 

interactions  must  be  part  of  an  accurate  description  of  the  V-V  energy 

exchange  process. 

DISCUSSION  AND  CONCLUSION 

The  "Fast-Flow"  approximation  used  in  the  determination  of  V-V 
energy  exchange  probabilities  results  in  a  substantial  simplification  of 
the  system  of  simultaneous  equations  which  describe  the  steady  state 
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vibrational  population  distribution.  Diffusion  and  wall  deactivation 
processes  can  be  made  insignificant  compared  to  the  volumetric  pumping 
rate  as  can  colllslonal  processes  other  than  those  of  interest.  For 
example  under  the  conditions  of  the  experiment  deactivation  of  vibrationally 
excited  CO  by  oxygen  or  helium  is  estimated  to  be  10*  times  smaller 
than  by  the  added  cold  CO.  Deactivation  by  water  formed  during  the 
combustion  of  acetylene  is  difficult  to  estimate.  Since  the  formation 
of  water  probably  occurs  in  a  reaction  step  subsequent  to  the  primary 
steps  leading  to  CO*  it  appears  likely  that  V-V  and  to  r.ome  extent  V-T 
equilibration  will  take  place  before  the  appearance  of  a  substantial 
amount  of  water. 


Knowledge  of  the  time  development  of  the  vibrational  populations, 
particularly  a  population  inversion,  is  critical  for  the  application  of  this 
reaction  in  a  chemical  laser.  A  time-dependent  solution  of  the  kinetic 
equations  of  the  vibrational  populations  has  been  carried  out  Hind  will  be 
discussed  in  another  paper ^10^.  The  solution  includes  all  single 
vibrational  quantum  exchange  collisions  between  levels  v  *  0  to  v  ■  14, 

V-T  and  spontaneous  radiative  pro 'esses.  The  computation  agrees  with  the 
data  of  reference  (5)  that  under  the  conditions  of  the  experiment  0.2  Torr 
added  cold  CO  can  give  rise  to  a  total  vibrational  population  inversion  while 
0.1  Torr  added  cold  CO  can  increase  local  vibrational  temperatures  but 
fails  to  produce  an  inversion.  The  computation  also  indicates  that  for 
the  CSj-O  system  a  substantial  amount  of  CO  in  the  v  »  0  level  is  required 
in  addition  to  the  initial  distributions  reported^  in  order  to 

approximate  the  observed  time  dependence  of  the  distributions  for  this 

_  (11) 

system 
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FIGURE  CAPTIONS 


FIGURE  1 

FIGURE  2 

TABLE  1. 

FIGURE  3 

FIGURE  4. 


CO  first  overtone  spectrum  showing  the  effect  of  adding 
cold  CO. 

A  schematic  diagram  of  the  experimental  setup. 

V-V  energy  transfer  probabilities  ..  and  the  values 

V  |V"1 

of  the  de-activation  rate  constants  K  ,  of  CO-CO  collisions. 

v,v-l 

Determination  of  the  de-activation  rate  constant  K  . . 

V  |  V**J. 

'Jollisional  V-V  energy  transfer  probabilities  of  CO.  The 
measured  data  is  compared  with  a  theoretical  calculation 
which  includes  both  long  range  and  short  range  Interactions 
at  T  -  300°K- 


-120- 


^  ° 
cr  ° 

o3 

to  ° 

2  CD 

o  ? 

CO  Q 

co  Q 
==  < 


L X 

O 


LU 

llJ 
2  h- 
O  O 
H  Ixl 


q: 

iJ 

> 

o 

I— 

CO 

ct: 


lx 

Ix 

LxJ 

IxJ 

X 

h- 


CO 

£ 
o 

O  X 
O  CO 


i 

•H 


,7- 


/ 


-122- 


T  C\  L 1  £  | 

V-V  ENERGY  TRANSFER  PROBABILITIES  Py°^  AND  VALUES  OF 
THE  DE-ACTIVATION  RATE  CONSTANTS  K„  „  ,  OF 
CO  (O)  +  CO  (v)  —  CO  (1)  +  CO(v-l) 

RESULTS 


Vibrational 

^EVEL _ Present  Results  Hancock  &  Smith's1 


V 

kv,v-i 

K  P  °>1 

V,V-jl  "  V  ,  V— 1 

3 

2,il3(-12)  8.5(-3) 

4 

1.8AC-12)  6.2C-3) 

2,0(-12)*  6.7C-3) 

5 

. 

1.19 (-12)  4.0C-3)  j 

1 . 3 (-12)  A,3(-3) 

6 

8.67C-13)  2.9C-3) 

6. 3 (-13)  2.1(-3) 

7 

5.2(-13)  1.7A(-3) 

3.'l(-13)  1.15(-3) 

8 

3 .62C-13)  1.21 (-3) 

1 .77 (-13)  5.9(-A) 

9 

2.38C-13)  8,06(-N) 

1,13 (-13)  3.8(-A) 

1  G,  Hancock  &  I.W.M,  Smith,  Appl,  Opt,,  10,  1827,  1971 
2,0(-12)  =  2,0  x  10  12  cm3  -  molecule”1  -  sec”1 
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Pie j^re  \  COLLISIONAL  V-V  ENERGY  TRANSFER  PROBABILITIES  OF  CO; 

THE  SOLID  LINE  IS  THE  THEORETICAL  CALCULATION  INCLUDING 
BOTH  LONG  RANGE  AND  SHORT  RANGE  INTERACTIONS  AT  T  =  300 
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ABSTRACT 

The  dynamic  behavior  of  a  relaxing  CO-He  mixture  has  been 
studied.  The  time  evolution  of  the  non  -  equilibrium  vibrational 
population  distribution  was  obtain  id  by  numerically  integrating  the 
master  equations  for  the  vibrational  population  including  Vibrational- 
Vibrational,  (V-V),  Vibrational -Translational,  (V-T),  and  the 
spontaneous  decay  processes.  The  result  shows  the  characteristic 
relaxation  through  the  V-V  energy  transfer  of  an  anharmonic  molecule. 
The  results  show  good  agreement  with  the  experimental  data  from 
chemiluminescence  experiments  in  the  following  two  reactions 
(i)  O  +  C2H2  -♦  CO*  +  CH2  +  51  Kcal/mole;  O  +  CH2  -»  CO*  +  2H  + 

71  Kcal/mole;  and  (ii)  O  +  CS2  -♦  CO  +  CS,  O  +  CS  -»  CO  +  S  + 

75  Kcal/mole. 
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INTRODUCTION : 

Studies  of  the  vibrational  energy  distribution  of  an  anharmonic 
molecule  in  non- equilibrium  states  have  received  increasing  attention 
in  recent  years.  This  is  true  not  only  because  of  the  relative  simplicity 
of  the  system  for  the  study  of  the  relaxation  mechanisms  in  a  diatomic 
molecule,  but  also  because  of  its  importance  of  the  understanding  of  the 
molecular  kinetics  in  a  chemical  or  molecular  laser  system,  and  for 
the  improvement  of  the  performance  of  such  laser  systems. 

Recently,  the  I.  R.  chemiluminescence  method  has  been  successfully 

applied  to  study  the  vibational  population  distribution  of  the  product 

(12  3) 

molecules  in  exothermic  reactions.  '  *  ’  'It  has  been  observed  that  the 
products  are  frequently  in  highly  vibrationally  excited  states  following 
the  reaction,  and  then  relax  through  a  rapid  V -V  energy  transfer  process. 
Studies  of  I.  R.  chemiluminescence  provide  information  on  the  vibrational 
energy  distribution  at  various  stages  of  the  relaxation  during  a  kinetic 
process.  However,  due  to  the  complexity  of  the  kinetics  involved,  there 
has  not  been  a  comprehensive  study  of  the  dynamic  behavior  of  the 
vibrational  relaxation  mechanism  in  a  diatomic  molecular  system. 

MODELING  OF  THE  RELAXATION  PROCESS 

In  this  work,  we  report  a  successful  modeling  scheme  for  study  of 
the  time-dependent  behavior  of  the  vibrational  energy  transfer  in  a  relaxing 
CO-He  mixture.  The  set  of  master  equations  for  the  vibrational  popu¬ 
lation  has  been  numerically  integrated  by  the  Runge-Kutta-Treanor 

(4) 

method.  The  time  evolution  of  the  vibrational  population  thus 
obtained  shows  good  agreement  with  our  data  measured  experimentally 
from  the  I.  R.  chemiluminescence  of  the  excited  CO  produced  in  the 
exothermic  reaction  of  oxygen  with  acetylene. '  '  We  have  been  able 
both  experimentally  and  computationally  to  demonstrate  that  a  total 
inversion  of  the  vibrational  population  occurs  as  a  consequence  of  the 
preferential  de-activation  of  low  vibrational  levels  due  to  the  presence 
of  an  excess  amount  of  ground  state  CO.  We  have  also  applied  this 
modeling  method  to  examine  the  vibrational  population  data  recently 

1. 
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reported  by  Tsuchiya  et  al.  ^in  a  chemiluminescence  experiment  on 
a  pulse -discharged  CS2/Oz/He  system.  The  results  show  a  reasonably 
satisfactory  agreement  although  some  discrepancy  exifcts  as  to  the  overall 
relaxation  rate.  The  time -dependent  solution  of  the  vibrational  popu¬ 
lation  was  obtained  numerically  by  solving  the  set  of  master  equations 
for  the  vibrational  population  including  V -V %  V-T  and  the  spontaneous 
radiative  process  with  a  given  set  of  the  initial  populations  experimentally 
determined  from  chemiluminescence  data. 

The  master  equation  we  have  used  for  the  vibrational  population  of 

the  level  V  is: 


N..(t)  -  Z „  _  __  J  y  |pK-l,,K  ^  I»j  n  1 

at  K  CO-CO  t-ofj +1,  j  NK-lNj+l  J ,  j  +1  jJ 

+  j?1[P^-V.jNK+INj-!  -Pj*.r-INKNj]  +  NCO(t>  [YK-1,KN1 


-  UK,K-1NK  +UK+1,K  NK+1  •  UK,K+1Nk]  +  AK+2,KNKt2+AK+l,KNK+l 

'  AK,K-1  NK  "  K-2NK  ,  (1) 

where  N„,  N.  are  the  population  densities  for  levels  K  and  j  at  time  t. 

K  J 

pK,  K-l  y_y  energy  transfer  probability  in  the  CO-CO  collisional 

j»j+1 

process; 

CO(j)  +  CO(K)  ->  CO(j+l)  +  CO(K-l)  +  AE. 

2. 
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(7) 

This  probability  was  calculated  by  combining  the  Sharma-Brau  theory' 
and  that  of  Rapp^  to  include  both  long-range  and  short-range  interactions. 
The  measured  probability  has  recently  been  reported^* ^  by  several 
workers  and  the  results  show  good  agreement  with  the  calculated  values 
thus  obtained  for  the  room  temperature  case.  *  ^ and  UK 

are  the  V-T  energy  transfer  probabilities  in  CO-CO  and  CO-He  collisions 
at  room  temperature,  and  the  values  employed  were  calculated  from  SSH 

theory.  Zcq  qo  and  ZCO-He  are  co^sional  frequencies  for 

CO-CO  and  CO-He  collisions.  N^Q<t)  is  the  concentration  of  CO  at  time  t. 

NHe  *8  of  He‘  K  1  an<*  K-2  are  t*ie  sPontaneous  radiative 

rates  of  the  level  K  to  levels  K-l  and  K-2  respectively.  In  the  above 

equation,  only  a  single  quantum  exchange  is  included  in  the  V-V  exchange 

process.  The  set  of  equations  (1)  for  the  vibrational  levels  V  =  0  to  V  =  14 

(4) 

has  been  integrated  numerically  by  the  Runge-Kutta-Treanor  method. 

This  numerical  integration  technique  has  been  widely  applied  for  solving 
coupled  non-linear  differential  equations  such  as  are  found  in  modeling 
coupled  chemical  reactions. 

Figure  (1)  shows  the  vibrational  population  distribution  of  CO* 

measured  experimentally  from  the  chemiluminescence  data  from  excited 

CO  formed  in  the  oxygen-acetylene  reaction.  ^5»10^  The  vibrational 

population  was  determined  from  the  first  overtone  chemiluminescence 

spectrum  between  2.3  and  2.8  microns^bv  using  a  computer  simulation 

method  first  described  by  Karl,  et  al.  A  numerical  method  is  used 

to  generate  a  simulated  spectrum  to  compare  with  the  experimental  scans. 

It  takes  into  account  the  instrumental  sensitivity  variation  as  a  function 

of  wavelength,  the  spectral  slit  width,  the  Einstein  A  coefficient  of  each 

upper  level  and  the  molecular  and  spectroscopic  constants.  The  rotational 

temperature  and  the  relative  population  of  each  vibrational  level  are 

adjusted  to  obtain  a  match  between  the  synthetic  spectrum  and  the 

experimental  scan.  More  detailed  discussion  of  this  procedure  was 

(5) 

reported  previously.'  '  The  figure  shows  the  non -equilibrium  relative 
vibrational  population  distribution  of  CO*  at  a  time  approximately 
100  microseconds  after  the  injection  of  the  acetylene  into  the  atomic 
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oxygen  stream.  The  addition  of  cold  CO  to  the  fuel  dramatically 
modified  the  vibrational  distribution.  Preferential  de -activation  of  the 
lower  vibrational  levels  leads  to  a  total  inversion  of  the  vibrational  popu¬ 
lation  with  0.2  Torr  CO  added  to  the  acetylene.  The  data  have  been 
plotted  taking  account  of  the  relative  signal  levels  under  different 
conditions  of  added  cold  CO.  It  is  seen  from  the  curve  that  as  a 
consequence  of  the  anharmonic  nature  of  the  molecular  potential,  the 
process 

CO(0  )  +  CO(V)  -*CO(l)  +  CO(V-l)  +  &E 

has  a  maximum  transfer  probability  between  V  =  4  and  V  =  6  and  decreases 
rapidly  at  higher  values  of  V.  With  the  addition  of  cold  CO  to  the  reaction, 
it  has  therefore  been  possible  to  significantly  depress  the  populations 
of  levels  V  <  5  and  produce  a  vibrational  population  inversion. 

Figure  2(a)  shows  the  calculated  time -evolution  of  the  vibrational 
population  model  for  the  conditions  shown  in  Figure  1  in  the  presence  of 
0.2  Torr  cold  CO.  The  initial  vibrational  population  was  taken  from  the 
experimental  data^'  ^  and  was  a  nearly- Boltzmann  distribution  similar 
to  the  top  curve  shown  in  Figure  1.  The  system  was  assumed  to  be  isothermal 
at  a  temperature  300°K  and  no  attempt  was  made  to  model  the  mixing 
mechanism  in  the  flow.  The  populations  are  seen  to  be  changing  more 
rapidly  for  vibrational  levels  V  <  8  and  a  total  inversion  occurs  at  a 
time  about  150  microseconds  after  the  collisional  relaxation  begins. 

The  populations  relax  through  a  series  of  non-Boltzmann  distributions 
with  the  high  vibrational  levels  over -populated.  As  the  relaxation  time 
reaches  1  millisecond,  a  Boltzmann  distribution  appears  for  level  V  <  8. 

This  is  shown  in  Figure  2(b). 

Figure  3(a)  shows  the  time -evolution  of  the  vibrational  population 
with  the  same  conditions  except  the  added  cold  CO  was  changed  to  0.1 
Torr.  It  is  seen  from  the  curves  that  the  vibrational  population  changes 
moderately  during  the  relaxation  process  but  a  total  inversion  does  not 
appear  in  this  case.  The  population  curve  corresponding  to  162.  9  (xsec 
after  the  start  of  relaxation  shows  very  good  agreement  with  the  curve  (A) 
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in  Figure  1.  Figure  3(b)  shows  the  long-time  behavior  for  this  case 

Figure  4  shows  the  experimental  results  of  a  study  of  the  relative 
populations  of  vibrationally  excited  CO  produced  in  a  pulse-discharged 
CS2  +  02  +  He  mixture  as  reported  recently  by  Tsuchiya  et  al.  ^  They 
deduced  the  relative  populations  from  the  chemiluminescence  of  CO*  at 
selected  delay  times  after  the  discharge  pulse  which  initiated  the  reaction  and 
established  the  time  dependence  of  the  population  for  each  level.  A 
modeled  time -evolution  of  the  vibrational  populations  for  this  case 
is  shown  in  Figure  5.  The  two  figures  show  a  similarity  of  the  time- 
dependent  behavior  of  the  relative  vibrational  populations.  However,  some 
discrepancies  exist.  First,  the  computed  relaxation  process  appears 
slower  than  that  determined  experimentally.  Since  in  our  model,  only 
V-V  exchange  by  CO-CO  collisions,  V-T  exchange  by  CO-CO,  and  CO-He 
collisions  and  the  spontaneous  radiative  decay  were  considered,  any  other 
de-activation  processes  will  speed  the  relaxation.  However,  our  modeling 
scheme  does  provide  a  method  to  verify  some  experimental  conditions 
assumed.  Second,  in  our  model,  initial  populations  for  the  ground  level 
and  the  first  level  must  be  determined  by  extrapolating  from  the 
populations  in  the  adjacent  levels  assuming  that  a  nearly  Boltzmann 
distribution  exists  for  these  low  levels.  Since  it  is  not  possible  to 
determine  the  population  for  the  lowest  level  by  a  chemiluminescence 
experiment,  the  question  of  how  these  lowest  levels  are  populated  in  the 
nascent  species  in  these  two  types  of  reaction  has  not  yet  been  solved. 

Failure  to  include  a  substantial  V  =  0  population  from  the  CS  -  O  reaction 
in  the  calculation  resulted  in  a  time  development  very  different  from  that 
observed  experimentally. 

DISCUSSION: 

Our  results  for  the  time -dependent  behavior  of  the  vibrational 
populations  have  shown  a  rather  satisfactory  agreement  with  the 
experimental  data  presently  available  from  chemiluminescence  experiments. 
The  model  includes  V-V  and  V-T  processes  and  radiative  decay  in  a 
relaxing  CO-He  mixture.  It  is  a  rather  simple  one,  and  the  numerical 
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integration  technique  we  employed  is  straightforward.  The  agreement 
with  the  experimental  data  suggests  that  this  modeling  scheme  is  widely 
applicable  to  general  kinetic  problems  in  other  diatomic  molecular  laser 
systems.  It  can  be  elaborated  with  more  complicated  relaxation 
mechanisms,  namely,  multiquantum  vibrational  exchange,  chemical 
reactions,  fluid  mixing,  and  stimulated  emission  to  model  the  kinetics 
of  a  transverse -flow  chemical  laser  system  which  is  of  particular  current 
interest. 
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FIGURE  CAPTIONS: 

FIGURE  1  Relative  vibrational  population  distribution  of  CO’'. 

Pressures:  C2H2  0. 05  Torr,  02  0.  68  Torr.  He*  1.67 

Torr.  Temperature  Trot=  3 00°K.  [C2H2  ]  =  0.  25  x  10"4 

mole/sec.  [0]=  1.7  x  10""4  mole/sec.  o  without 
additional  cold  CO.  a  0.  05  Torr  additional  cold  CO. 

6  0.1  Torr  additional  cold  CO.  ®  0.2  Torr  additional 
cold  CO. 

FIGURE  2(a)  The  time  evolution  of  the  vibrational  population  of  a 

relaxing  CO-He  mixture  in  the  presence  of  0.2  Torr  cold 
CO.  The  partial  pressure  of  He  is  1.  70  Torr:  that  of 
CO  is  0.25  Torr.  T  =  300°K.  The  calculation  includes 
V-V,  V-T  and  radiative  decay  processes. 

FIGURE  2(b)  The  long-time  behavior  of  2(a). 

FIGURE  3(a)  The  time  evoluation  of  the  vibrational  population  of  a 
relaxing  CO-He  mixture  in  the  presence  of  0.  1  Torr 
cold  CO.  The  partial  pressure  of  He  is  1.  70  Torr;  that 
of  CO  is  0.  15  Torr.  T  =  300°K.  The  calculation  includes 
V-V,  V-T  and  radiative  decay  processes. 

FIGURE  (b)  The  long-time  behavior  of  3(a). 

FIGURE  4  Vibration  population  distribution  of  CO*  produced  in  pulsed - 

discharged  CS2  +  02  +  He  mixtures.  The  time- 

dependence  of  the  population  is  determined  by  the  delay  time. 

The  partial  pressure  of  CS2  is  0.045  Torr;  that  of  °2  is 

0.  76  Torr  and  He  is  3.0  Torr.  15  KV  discharge.  Data  taken 
from  Tsuchiya,  et  al.  (&) 

FIGURE  5.  Calculated  time -evolution  of  the  vibrational  populations 
for  the  experiment  of  Figure  4. 
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V.  Molecular  Relaxation  via  Atomic  Collisions 

Three  studies  are  currently  in  progress  to  measure  atcm  collision 
Induced  vibrational  relaxation  in  Important  molecular  laser  systems • 

The  collisional  processes  being  studied  are: 

HF(v  -  1)  +  X  ■+  HF(v  -  0)  +  X  +  AE  ;  X  -  H,F 

C02(0n0)  +  X  -*■  C02(000)  +  X  +  AE  ;  X  -  O.N.H.F 

C02(00°1)  +  X  -*>  C02(nm  )  +  X  +  AE  ;  X  »  H,F,0,N 

A.  Atomic  Relaxation  of  HF(v  -  1)  (G.  J.  Wolgs) 

At  the  time  of  submission  of  our  Semi-Annual  Report  for  1972  we  discussed 
the  status  of  this  experiment.  In  this  report  we  shall  complete  the  descrip¬ 
tion  of  the  experiment  which  is  now  in  progress. 

The  capacitance  diaphragm  pressure  gauge  has  proven  to  be  very  useful 
and  reliable.  In  the  following  figures  we  show  a  schematic  drawing  of  the 
gauge  and  typical  calibration  runs.  The  gauge  is  constructed  from  materials 
(Telfon,  PVC,  Monel)  that  are  corrosion  resistant  and  no  difficulty  is 
anticipated  in  applying  it  to  our  work  with  HF,  F2,  etc. 

The  HF,  v  •  1  -*■  0  pump  laser  was  redesigned  and  significantly  improved. 

We  now  obtain  approximately  20-30  milli  joules/pulse  on  the  v  -  1  -*■  0 
transition.  This  will  be  more  than  adequate  for  our  experimental  needs. 

The  fluorescence  cell  was  designed  and  constructed.  A  schematic  drawing 
of  the  cell  is  attached  as  a  figure.  The  cell  design  stressed  use  rf  construc¬ 
tion  materials  that  could  withstand  the  corrosive  gases  and  maintenance  of 
a  good  filling  factor  for  the  pump  radiation  and  good  optical  coupling  for 
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detection  of  fluorescence. 

We  have  made  studies  of  the  atom  wall  recombination  rate  for  H  atoms 
in  our  flow  system.  Knowledge  of  this  loss  mode  for  H  is  essential  to  the 
quantitative  studies  we  have  undertaken, 

At  the  time  of  preparation  of  this  report,  experiments  studying  the 
process  HF(v  *  1)  +  H  -*■  HF(v  ■  0)  +  H  +  AE  were  in  progress. 

B.  Atomic  Relaxation  of  CC^fOnO)  (S.  H.  Bauer) 

Construction  of  a  high  power  TEA  CO2  laser  is  nearly  completed.  We 
plan  to  use  this  unit  for  simultaneous  irradiation  of  CC^-M  mixtures  with 
4.2  pm  (HBr  laser)  and  10.6  or  9.6  ym  (CO2  laser)  pulses  in  order  to  obtain 
(concurrently)  substantial  populations  of  C02  molecules  in  the  (00°1)  and 
(02°0)  or  (10°0)  states  and  to  measure  the  fluorescence  lifetimes  at 
4.2  pm  and  15.5-16.5  ym.  These  decay  times  will  permit  us  to  obtain 
transition  probabilities  for  the  various  collisional  vibrational  deexcitation 
steps,  at  selected  cell  temperatures,  with  specified  colliders  (M) .  Indeed, 
we  anticipate  eventually  to  go  over  to  flow  cells  for  measuring  the  relative 
efficiencies  of  atomic  relsxers. 

C.  Atomic  Relaxation  of  CC^fOO0!)  (G.  J.  Wolga) 

A  parallel  experiment  to  that  of  part  A  has  been  set  up  to  study  the 
process  C02(00°1)  +  X  +  CC^Cnm  )  +  X  +  AE  ;  X  -  H,F,0,N  .  A  second 
EPR  apparatus  for  measuring  atom  concentrations  is  in  operation.  The  pump 
laser  for  the  laser  induced  fluorescence  studies  will  be  shared  with  the 
HF  deactivation  experiment  except  that  it  will  be  operated  as  an  HBr  chemical 
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